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The probability function of Menke has been used as a 
basis for a statistical treatment of the intermolecular po- 
tentials in solutions of symmetrical molecules for which the 
entropy of mixing is the same as for an ideal solution of the 


same composition. The resulting expressions enable calcula- 
tions to be made of the departure of iodine solutions from 
the ideal solution laws which agree satisfactorily with the 
observed departures. 





HE senior author has called attention! to 

the existence of families of solubility curves 

for a single solute in a number of solvents where 
both components are comparatively nonpolar 
and the chemical effects of solvation, association 
and compound formation are absent. Such solu- 
tions have been called regular, and it has been 
suggested that the regularity is due to a com- 
pletely random distribution of the component 
molecules resulting from thermal agitation, such 
as exists in ideal solutions. Such randomness 
cannot, of course, persist to indefinitely low tem- 
peratures in the presence of orienting or associat- 
ing forces or unequal molecular fields. On lower- 
ing the temperature to the critical solution 
temperature, for example, the two species begin to 
unmix, giving a very different statistical distri- 
bution, and submicroscopic aggregates certainly 
form above this point; however, the range of 
temperature and concentration in which this ag- 
gregation occurs to any appreciable extent is 
certainly very small, as shown by theoretical and 
experimental studies of opalescence in this re- 
gion.? Where this randomness exists, it suggests, 





1 J. H. Hildebrand, J. Am. Chem. Soc. 51, 66 (1929); also 
Proc. Nat. Acad. Sci. 13, 267 (1927). 

* Cf. M. v. Schmoluchowski, Ann. d. Physik [4] 25, 205 
(1908); A. Einstein, ibid. 33, 1275 (1910); Wo. Ostwald, 
ibid. 36, 848 (1911). 
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because of the well-known connection between 
entropy and probability, that the entropy of 
mixing for a regular solution should be the same 
as for an ideal solution of the same mole fraction. 
Writing the partial molal entropies of component 
X,in a regular and in an ideal solution as 5,’, and 
Si‘, respectively, the above conclusion can be ex- 
pressed by the equation, 


$;"-—S,'= . (1) 


Since the corresponding partial molal free ener- 
gies are related to the entropies by the equation, 
0(F;’ —F,')/8T = —(S:"—S:'), we can integrate, 
and write 


Fi’ —F,'=f(N), (2) 


where N denotes mole fraction. The form of f(N) 
is the subject of this paper. It will be approached 
through the consideration of the energy of mix- 
ing. This is readily possible for regular solutions, 
since, when As=0, AF=AuH, and since, further, 
An = AE when A(PV) is negligible. 

In order for a solution to be ideal, the potential 
of a given molecule must not alter when one 
species among the surrounding molecules is re- 
placed by another. A solution will be regular even 
if this molecular potential varies, provided that 
the randomness of arrangement is still ‘‘ideal,”’ 
i.e., complete. This can hardly be the case, in 
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general, unless the molecules are symmetrical and 
nonpolar. (Isolated cases may be cited where a 
change of species, even with polar molecules, 
would probably not materially alter the arrange- 
ment, as with two isomeric esters.) We will as- 
sume equal molecular volumes, which allows a 
more satisfactory statistical treatment, before 
taking up the more difficult case of unequal 
volumes. 

The form of f(N) has been deduced by van 
Laar,? from van der Waals theory, and by 
Heitler* by aid of the assumption that the liquid 
has a cubic lattice structure. The former treat- 
ment shares the unconvincing character of the 
van der Waals equation, and the latter depends 
upon a picture which cannot be true. Heitler 
failed, moreover, as shown by his chosen ex- 
amples, to recognize the essential limitation to 
systems that we have called regular. The follow- 
ing derivation is, nevertheless, to be regarded as 
a refinement of his method. 

The nature of liquid structure has been re- 
vealed by the studies of x-ray scattering made 
recently by H. Menke‘ in Debye’s laboratory. 

The structure of a liquid is expressed by a 
probability function, W (designed W, by Menke; 
the subscript is here omitted for the sake of sim- 
plicity) which is defined as follows. The probabil- 
ity that an element dV of a volume of liquid, V, 
will contain the center of a given molecule (e.g., 
the nucleus of a central atom) is dV/V. The 
probability that a given pair of molecules will 
occupy two such elements is (dV/V)*, provided 
they are sufficiently far apart that they do not 
interfere with each other. When, however, the 
two elements of volume are very small, and their 
distance apart, 7, is of the same order of magni- 
tude as the molecular radius, the volume of the 
molecules themselves will help to determine 
their positions, and the probability must be ex- 
pressed by W(dV/V)?, where W is a function of r. 
Fig. 1A gives the plot of this function for mercury 
as determined by Menke from x-ray scattering. 

It is possible to use W to calculate the number 


3 J. J. van Laar, Zeits. f. physik. Chemie A137, 421 (1928). 

4 Heitler, Ann. d. Physik [4] 80, 630 (1928). Cf. also J. H. 
Hildebrand and E. J. Salstrom, J. Am. Chem. Soc. 54, 4257 
(1932). 

5 (a) H. Menke, Phys. Zeits. 33, 593 (1932); cf. further, 
(b) O. Kratky, ibid. 34, 482 (1933). 
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Fic. 1. A, Menke probability function for mercury; B, 
number of atoms within distance r of a central atom, for 
liquid and solid; C, potential function, area under which 
indicates small contribution to heat of vaporization of 
molecules beyond the “‘first layer.” 


of molecules (atoms in the case of mercury) sur- 
rounding a central one at various distances. The 
number of nuclei per cc in the liquid as a whole is 
N/v, where N is the Avogadro number and v the 
molal volume. The number of molecules in a 
large sphere of liquid would of course be 
(N/v)(4/3)7r. In a small sphere, however, the 
number surrounding a central one is 


(w/v) f 4rr?Wadr. 
0 


Fig. 1B shows this number for various values of 
r calculated by graphic integration. The distribu- 
tion for solid mercury is given for comparison 
(without correcting for the difference in density) 
and shows very strikingly how far the structure 
of the liquid differs from that of the solid. 

The number in a spherical shell of radius r and 
thickness dr is (4r%N/v)Wr’dr. The number of 
molecules that can be regarded, in succession, as 
“central’’ to a shell is NV, hence the total number 
of pairs at distance r from each other is (27 N?/V) 
Wr'dr (the numerical coefficient is changed from 
4 to 2 to avoid counting each pair twice). If the 
potential of each pair is ¢ (a function of r) the 
potential of all the pairs at distance 7 is 
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(27N?/v) ¢Wr'dr, and the potential energy of the 
whole liquid, per mole, is 


E” = (2x N2/v) f yWradr. (3) 
0 


If we assume with London,* and others, that 
the attraction potential between molecules is 
given by —k/r®, and neglect the repulsive poten- 
tial, the potential energy of a mole of liquid 
would be 


EE” = (—2en*k/v) [ (W/r)ar. (4) 


Curve C, Fig. 1, shows the curve of W/r‘ 
plotted against 7, the area under which represents 
the above integral. It is here given to illustrate 
how little is contributed by molecules beyond 
those immediately adjacent. 

We will next consider a solution of +72 
moles of species X; and X2, whose molal volumes 
are identical, v, and which mix without volume 
change. The number of molecular centers of both 
species per cc is the same as for the pure sub- 
stance, i.e., N/v; the probability function should 
be the same for both species, hence the total 
number of molecules of both kinds at distance r 
from a central molecule is (44N/v)Wr'dr, and 
since there are (%,;+72)N possible central mole- 
cules, the total number of pairs at distance r from 
each other is (27.N?2/v)(1-+2) Wr’dr. The chances 
that these pairs are of the respective species 
X1Xi, XoXo and X1Xe are n;"/(m+N2)?, n2"/ 
(m+n)? and 2nin2/(m1+n2)?, hence the number 
of pairs of each kind are, respectively 

















2rN? n;? 2xrN? ne? 
waimcomieds -Wrdr, . -Wr'dr 
Vo m+ne Vo m+n. 
and 
4nN? NN 
-Wr'dr. 
Vosm+Ne 


The potentials of single pairs are functions of 7 
of the same form, designated by ¢u, ¢22 and ¢i2, 
where the subscripts denote the species forming 
the pair. The potential of all the X1X1 pairs 
at distance r is then (2rN?/v)-[m?/(m-+m2) ] 
‘giWr'dr, and the potential of all Nuk pairs 





°F. London, Zeits. f. Physik 63, 245 (1930). 


throughout the liquid is (27 N?/v) - [m:?/(m+72) ] 
- So? euWrdr, and similarly for the X2X»2 and 
X,X» pairs. Adding these gives the potential 
energy of the entire liquid of m:+2 moles, desig- 
nated by EF’. The partial molal energy of X, in 
this solution is E,;= (0E’/dm),,, which gives 


2rN? 1 
Vv (11+)? 


xf euWerdr—nd [ v2 Wrdr 
0 0 


+n f eaWrar], (5) 
0 





E,’ = 


(m+ 2n ne) 


Subtracting §£,°’ for a mole of the pure Xj, 
setting E;/—E," = E,—E,°, the difference in poten- 
tial energy equal to the difference in total 
energy, and collecting yields 


E, —E,°= Qnwne/v)|2f/ g2Wrdr 
0 


-f euWrdr— eaWr'dr] (6) 
0 0 


If we set gu=—h1/r*, etc., we obtain an ap- 
proximation formula 
4 27 N?N2? oo W 
E, — E)° =—_———_ (Ri + Roe — 212) —dr. (7) 
V . ¥ 


If, as has been frequently assumed, ie 
= (Rikee)?, then 


ku +t+-Roe — 2ki2= [ (R11)! — (Ree) *)?. (8) 


By setting up the equations in the form of Eq. 
(3) for both pure components it can be seen that 
Eq. (7) can be transformed into 


E, — E,° = (E1°)* — (E2°)* Py’. (9) 


This, unlike 6, can hardly be a rigid equation on 
account of the additional assumptions involved, 
especially Eq. (8), and it is not our purpose to 
subject it to a searching experimental test at this 
juncture. One illustration only will be given. 
Iodine and carbon disulfide have nearly identical 
molal volumes, 59.0 cc and 60.7 cc, respectively, 
at 25°C. Both have zero dipole moment. The 
solubility of iodine in carbon disulfide agrees 
closely! with the equation In N,‘/N:=)N2?/RT7, 
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where 0 is a constant. This corresponds to the re- 
lation F;’—Fi;'=bn.*. The value of 5 calculated 
from solubility is 870 cal. The energies of vapor- 
ization at 25° are, respectively, 11,300 cal. 
and 5970 cal., from which we calculate 
[ (Ey°) #— (E_°)* P= 840 cal., a sufficiently good 
agreement. 

Eq. (6), at any rate, should correspond well for 
systems conforming to the conditions postulated. 
The function, W, should depend chiefly on the 
molecular radius, for symmetrical molecules, and 
be relatively independent of the species of mole- 
cule. It should flatten out as the temperature in- 
creases and the departure from hexagonal close 
packing®» becomes greater. 

Let us now consider a solution of two species 
whose molecular volumes are different, and 
where a single common probability function can 
not suffice. Consider a spherical shell of thickness 
dr, at a distance 7, from the center of a central 
molecule. The volume of this shell is dV =4mr*dr. 
The number of X; molecules whose centers are 
in 1 cc of liquid is m,:N/V, where V is the total 
volume of the solution. The shell under considera- 
tion would contain ,;NdV/V molecular centers 
of X, if there were no interference due to the size 
of the central molecule, but since such inter- 
ference exists a probability function will have to 
be introduced for each combination of species. 
This will be designated by subscripts, Wi for the 
X,X, combination, etc. 

The total number of molecules in the solution, 
each able to play the réle of ‘‘central’’ molecule, is 
(m,+n2)N. The number of central molecules of 
X, is m,N. The number of X; molecules in the 
shell surrounding each is mNWydV/V. The 
number of X,X; ‘‘pairs’”’ throughout the solution 
at distance r is 4n",N?WudV/V. (The 3 is intro- 
duced to avoid counting each pair twice.) The 
potential of each X1X, pair is gu, hence the po- 
tential of all such pairs at distance r is 
(2rnN?/V)+guWur?, and the potential of all 
X,X;, pairs throughout the liquid at all distances, 
is (2rN°n?2/V) So® euWur'dr. We will abbreviate 
this expression by writing for the integral the 
symbol /f(11), giving (27N?n,*/V) f(11). For 
the X2X-2 we get the corresponding expression, 
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(2N2n2?/V) f (22) and for the X:X2 and X2X, 
pairs, which are identical, (4r7N?nn2/V) f (12). 
We will pause here to consider the probability 
functions, Wi, Woe, Wie and W2. It would seem 
that the position of the first maximum, as shown 
in Fig. 1A, should be the same for the solution as 
for the pure substance, since it depends chiefly on 
the sum of the two molecular radii, increased 
somewhat by thermal agitation. The position of 
the second maximum would, however, be differ- 
ent, for, if X» is thought of as smaller than X,, 
two molecules of X; might be separated by an X2 
instead of another X;. The frequency of this 
separation would depend upon the proportion of 
X:2 molecules, hence W;; for the solution would 
be a function of the composition beyond the first 
maximum. A quantitative formulation of this 
dependence seems very difficult and we shall, for 
the present, take refuge in the fact that the poten- 
tial is but little affected by the course of W much 
beyond the first maximum, as illustrated by the 
area under the curve in Fig. 1C, permitting us to 
set Wy for solution and pure X; equal to each 
other as a first approximation. We are undertak- 
ing an experimental study of these functions. 
We should note, further, that although Wi. and 
Wa might appear to be different, the number of 
X,X_2 and X2X, pairs is the same, hence / (12) 
= f{(21). We are now prepared to write for .the 
total potential energy of 1+ moles of solution, 


B!= (2aN*/V)] mt f (11) +me f (2) 


+2nm f (12) (10) 


The partial molal energy of Xi, i.e., E;’, is ob- 
tained, as before by differentiating with respect 
to m, at constant m2, remembering that V=mvVi 
+n2V2, approximately. We should write, more 
strictly, V=mVi+meVe2, where V; and Vz are the 
partial molal volumes, but this could hardly be 
worth the added labor in view of the other ap- 
proximations involved. Carrying out this differ- 
entiation and subtracting the value of E,” from 
Eq. (3) gives 


ks N2Ve2 2 2 1 1 ; 
f—2'=2n07(—~_) vf — fan-— fan-— fea]. (11) 
m1VitMNeV2 ViV2 vi? V2" 
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We may regard this equation as valid for any 
system of the sort here considered. To apply it 
rigidly will, however, require more data than we 
now have concerning probability functions and 
intermolecular potentials. Indeed, one of its prin- 
cipal uses is likely to be the throwing of light 
upon the latter question. In the meantime we 
may derive from it an approximation formula 
which may be tested by the aid of solubility data 
now at hand. 

First, we will neglect the repulsive potential 
and write for the attractive potential simply 
= —k/r®, This may make an appreciable error in 
dealing with mercury at ordinary temperatures,’ 
which is not much expanded over its volume at 
0°K. Fig. 1B shows the average number of 
molecules in the liquid which are closer together 
than in the solid. For substances at temperatures 
nearer their boiling points the effect of repulsion 
would be less. 

Next, we will assume, by analogy with Eq. (8), 


that 
faa=f (11) x f (2). (12) 


The physical significance of this is as follows. 
Since the form of the function, W, is the result of 
a thermally disturbed close packed arrangement, 
it may be assumed to be, for all spherical mole- 
cules, approximately the same function of the 
molecular domain. This could hardly be strictly 
true, for even the curve for a pure liquid will 
surely flatten with increasing thermal agitation, 
hence curves for different liquids must be, in 
general, slightly different. However, if we neglect 
differences in flattening as an effect of higher 
order, we may express the different W’s as the 
same function of the position of the first maxi- 
mum, r, for each, i.e., W=f(r/r). The curves ob- 
tained by Menke for mercury and gallium, when 
plotted against r/r, agree very closely, probably 
within the limit of error. We will therefore as- 
sume that the integral k fo° (W/r‘)dr, if trans- 


formed into 
k 7° sryi (7 r 
“J, )G)G)- 
redo \r r r 


would be the same for all substances so far as the 





7 J. H. Hildebrand, Phys. Rev. 34, 984 (1929). 


integrand is concerned. Introducing this, with 
subscripts appropriate to each combination, into 
Eq. (12), and cancelling out the integrand, leaves 


ky2/Ti2° = (Riskeo/Ti*t20*)}. (13) 


If we write fi2=43(fu+fe2), which should be very 
nearly true, Eq. (13) becomes 


(ru+f22)* 
kig =—————- (Rurkne). (14) 


8r11°Fo0" 


The coefficient of ki:k22 does not differ much from 
unity, even when ry and fee differ considerably; 
for example, it is only 1.06 when ry and fee 
differ by 50 percent. Eq. (14) is therefore nearly 
identical with Ey. (8), which is generally regarded 
as approximately true. We are therefore en- 
couraged to try the effect of combining Eq. (12) 
with Eq. (11), and then combine the result with 
Eq. (3) applied to each pure component. The 
final result is 


E —E;,° 


Ne2Ve 2 E;° 4 E2” oe 
ont a 
N1VitNeVe2 Vi Ve 
The equation derived by van Laar* with the aid 
of the van der Waals equation reduces to this 
under the conditions herein postulated. The same 
equation has been given by Scatchard? by the aid 
of some of the same assumptions as here em- 
ployed, plus the arbitrary one that the cohesive 


energy of a mole of mixture is given by the ex- 
pression 





AnVP ny? + 2a2V1V2N1N2+A22Ve"nN2” 





(m1Vi+M2V2)” 

where the a’s are regarded as constants (equiva- 
lent to the van der Waals constants.) 

Experimental data now being obtained in this 
laboratory upon tetrahalides and their solutions 
will enable us to make a more searching test of 
the laws of regular solutions than is now possible. 
In the meantime we may cite certain iodine solu- 
tions as most nearly conforming to the require- 
ments. We have data for the solubility of iodine® 


8 Geo. Scatchard, Chem. Rev. 8, 321 (1931). 

9 J. H. Hildebrand and C. A. Jenks, J. Am. Chem. Soc. 
42, 2080 (1920); G. R. Negishi, L. H. Donnally and J. H. 
Hildebrand, ibid. Nov. (1933); cf. also reference 1. 

























in the tetrachlorides of carbon, silicon and titan- 
ium, all very symmetrical molecules with zero 
electric moment. To these we shall add the solu- 
tions in carbon disulfide. Table I gives values of 
F.’—F,' for iodine, calculated from the experi- 
mental data by the expression RT In (Ne'/Ng), in 
comparison with values of E2”—E2°, calculated 
from Eq. (12). The values of vi, Ve, E:° and E2° 
used are indicated in Table I. They are taken 
from the usual sources which need not here be 
elaborated, since the purpose is merely illustra- 
tive. The values are given for 25° only, since the 
uncertainty regarding the change of the heat of 
vaporization with temperature hardly makes it 
worth while to attempt to check the small 
temperature trend. 

The agreement between the figures in the last 
two columns may be regarded as very satis- 
factory, indicating that the solubility of iodine 
could be calculated by the aid of Eq. (12) with- 
out serious error. 
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TABLE I, Test of Eq. (12) by means of solubility data for 











iodine at 25°C. 
Sol- Vcc E° Molepercent Fe—F,* E2.—E2° 
vent per mole cal. I, in soln. cal. cal. 
TiCl, 110.5 8370 2.150 1365 1535 
SiCl, 115.4 6610 0.499 2220 2320 
CCl, 97.1 7088 1.594 1725 1635 
CS. 60.7 6070 5.76 770 815 


Iodine 59.0 11330 








The solubility data for sulfur and for stannic 
iodide might also be used, but uncertainties re- 
garding heats of vaporization at the tempera- 
tures in question diminish the significance of the 
comparisons, and they are here omitted. 

One of the most important uses for the equa- 
tions here derived will doubtless be to throw 
light upon the important question of the poten- 
tial functions for pairs of like and unlike mole- 
cules. We expect to be able to report upon this in 
the near future. 
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Technique for the Electrolytic Production of H*H’O 


HuaGu S. TAYLorR, HENRY EYRING AND ARTHUR A. Frost, Frick Chemical Laboratory, Princeton, N. J. 
(Received October 2, 1933) 


A detailed description of a method of preparing heavy water in quantity is given and yields obtained over a 


period of months reported. 





UMEROUS requests for a detailed de- 

scription of our methods for the production 
of water rich in heavy hydrogen isotope, first 
outlined elsewhere,! prompt the following more 
detailed discussion of our present practice in the 
light of experience attained during the past 
summer in the treatment of 600 gallons of used 
electrolyte with the production of some 80 cc of 
water containing upwards of 90 percent heavy 
hydrogen. 


ELECTROLYSIS WITHOUT RECOVERY 


Electrolyte from commercial electrolytic cells 
is distilled in order to remove the large amounts 
of carbonate and hydroxide present. Sodium 
hydroxide is added to make the distillate 
approximately 0.5 molar. This solution is electro- 
lyzed in a battery of 210 cells, made from 
hydrometer jars, 4 cm diameter and 25 cm high, 
each of approximately 200 cc capacity, water- 
cooled by immersion in a tank containing flowing 
water and carrying a current of six to seven 
amperes. The electrodes are made of nickel 
plate, a strip 60 cm by 3 cm being bent twice at 
right angles so as to form anode and cathode in 
neighboring cells. The battery is subdivided into 
units of 30 cells each, the potential of the direct 
current supplied from a motor generator being 
adjusted so as to yield the amperage employed. 
After three days, when the electrolyte has 
diminished to 1/6 or 1/7th its original volume, 
it is replaced by a fresh batch of electrolyte. 
This process is so conducted that, while one cell 
is being emptied with the aid of suction, a 
neighboring cell is being filled with fresh electro- 
lyte by a siphon, the whole process requiring 





‘H. S. Taylor and H. Eyring, Proc. Am. Phil. Soc. 72, 
255 (1933). 
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about one hour. The concentrated electrolyte 
from the first electrolysis is partially neutralized 
by bubbling carbon dioxide through it so as to 
form the carbonate, with fuchsin as an indicator, 
and then distilled. This water is now added to 
another group of cells containing water of the 
same grade but still containing all of the sodium 
hydroxide originally added with the water at 
the start of the first electrolysis. This procedure 
obviates setting aside a portion of the concen- 


‘trated liquor before distillation which would 


later be added to the distillate to give the 0.5 
molar solution for the second electrolysis. In 
this way three successive electrolyses are carried 
out resulting in water of density D7" =1.001, 
approximately, and which therefore contains 
about 2.5 percent of the heavy hydrogen isotope. 


ELECTROLYSIS WITH RECOVERY 


From this stage onwards the hydrogen evolved 
during electrolysis contains an appreciable 
amount of the heavy isotope which should be 
recovered. This has been effected in the following 
way. Water containing sodium hydroxide is 
electrolyzed in cells of 200 cc capacity carrying 
a current of 10 amperes each and cooled in a 
trough of running water. The mixed electrolytic 
gas passes successively, as shown in Fig. 1, 





()<— recovered water 


Fic. 1. Electrolytic concentration apparatus. 
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through a spray-trap, a mercury bubbler, 
through a T-tube immersed in mercury to form 
a safety valve for releasing excessive pressure, 
and finally flows out through a capillary nozzle 
of Pyrex glass where it burns as a flame. The 
water formed by the combustion is condensed in 
an inclined quartz tube surrounded at the lower 
end by a water jacket. A red-hot platinum 
filament is placed just beyond the flame in order 
to re-ignite the gas in case the flame should go 
out, which may happen every few minutes if the 
size of the capillary nozzle is not perfectly 
adjusted for the rate of gas flow from the cells. 
For a given rate of gas flow, if the capillary is 
too fine it is easily clogged up, presumably by 
alkali from the cells; on the other hand, if the 
capillary is too wide, the flame strikes back 
through the tubing but is prevented from reach- 
ing the cells through the mercury obstruction 
(about 1 cm) in the mercury bubbler. If the 
condenser tube in the neighborhood of the flame 
becomes too hot, the flame jumps to the platinum 
wire where it burns incompletely. This is pre-’ 
vented by cooling the tube slightly at this point, 
a fact that indicates the effect is due to an 
excessive pressure of water vapor which is 
known to hinder the explosive reaction. For two 
cells each carrying 10 amperes the optimum size 
for the capillary has been found to be about 
0.3 mm. Each nozzle has a lifetime of about one 
day but is easily replaced by a new nozzle. 
Practically a quantitative recovery of the water 
can be obtained by this method. 


YIELD 


In electrolyzing a dilute solution of the heavy 
water down to about 1/6th of the original 
volume, it is found that the concentration of the 
heavy hydrogen in the heavier 1/6th has in- 
creased by a factor of about four while the 
concentration in the lighter 5/6ths is about 
2/5ths that of the original. The recovered water 
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TABLE I, 
Water obtained from Density % hydrogen 
electrolysis No. D° which is heavy 
I 0.998 — 
II 0.999 0.5 
Ill 1.001 ae 
IV 1.007 8 
V 1.031 30 
Vi 1.098 93 
VII 1.104 99 








is therefore added to the water one stage back in 
the process. With the present arrangements 
seven successive stages suffice to produce water 
which contains about 99 percent of the hydrogen 
in the heavy form. Table I indicates the quality 
of the water now obtained at each stage of the 
process. 

Table II shows the amounts of water put 
through each stage of the process during the 
period May 9th to September 27th, 1933. 


TABLE II. 








Amt. electrolyzed 


Electrolysis No. May 9-Sept. 27, 1933 





I 610 gal. (of commercial electrolyte) 
II 90 gal. 

III 521 

IV 10.15 | 

V 2.00 1 

VI 420 cc 

VII 82 cc 








It is found that about 1 cc of water, 95 percent 
heavy, is obtained for every seven gallons of 
starting material. Assuming that the commercial 
electrolyte liquor contains one part in 3000 heavy 
hydrogen, it is easily calculated that approxi- 
mately 10 percent of the heavy isotope of 
hydrogen present in the original starting ma- 
terial is obtained in the heavy water. Most of 
the loss occurs in the first three electrolyses 
where no attempt is made to recover the heavy 
hydrogen in the gas evolved during electrolysis. 
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The Ground State of the Hydrogen Molecule’ 


HuBert M. JAMES AND ALBERT SPRAGUE COOLIDGE, Harvard University 
(Received October 3, 1933) 


The method used by Hylleraas in treating the He atom 
has been extended to the Hz molecule. The method consists 
of setting up a wave function as a series in the five variables 
required, electronic separation being introduced explicitly 
as one of the variables. The coefficients are then deter- 
mined so as to produce the lowest energy. The energy found 
is within 0.03 v.e. of the most probable experimental value, 
while the form and location of the potential energy curve 


for various internuclear distances agree with those deduced 
from spectra to within similar limits. The value of the func- 
tion is computed for several configurations of the electrons, 
and compared with other approximations. Application of 
the method to other problems is discussed. A method is 
given for the numerical solution of secular equations of 
high degree. 





INTRODUCTION 


N constructing an approximate wave function 
for the hydrogen molecule in its normal state 
(regarding the nuclei as fixed), one may follow 
several lines of attack. Heitler and London make 
use of the unchanged functions of the individual 
atoms, regarding the interatomic forces as small 
perturbations. In a molecule as close-coupled as 
He this implied persistence of the identity of the 
individual atoms is quite unjustified, and the 
resulting wave function is a very poor approxi- 
mation. In particular, it is inadequate in that it 
is a function of only four electronic coordinates, 
instead of the required five. Aside from con- 
venience, the only advantage of the H—L scheme 
seems to be that it yields a result in which the 
energy of the molecule is represented as a sum 
of terms, some of which are just the energy of 
the separated normal atoms, so that these may 
be cancelled out and the remaining terms called 
‘binding energy.” 

Several attempts have been made to improve 
the H-L method by grafting upon it the variation 
principle; the original atomic wave functions are 
modified by the introduction of one or two 
arbitrary parameters, which are then determined 
so as to give the lowest energy. Such schemes do 
give somewhat better results, but do not escape 
the essential shortcomings of the original method 





‘A preliminary notice appeared in Phys. Rev. 43, 588 
(1933). Some minor numerical inaccuracies in that notice 
have been corrected in the present paper. 
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—the implication of atomic individuality, and 


the omission of an essential coordinate. Indeed, 
they actually sacrifice the peculiar advantage of 
the pure perturbation method, since the un- 


changed atomic energy no longer appears in the 


result. 


In view of the theoretical weakness and 


practical inadequacy of these methods, it seemed 
to us worth while to explore rather carefully a 
third possibility. Abandoning all concern with 
individual atoms, we have tried to build up from 
the ground a suitable molecular wave function, 
containing the full number of required coordi- 
nates. While guided by certain considerations of 
limiting forms of this function, we have placed 
our main reliance on the variation principle, 
introducing a large number of parameters to be 
determined by its use. 


DiscUssION OF METHOD 


Concerning the true function, we know the 


following: The ground state being a = state, 
the function must have rotational symmetry 
about the nuclear axis; that is, it may contain 
only the difference between the azimuthal angles 
of the two electrons, and not either angle alone. 


It must be symmetrical in the two electrons and 
also in the nuclei. In those parts of phase-space 
corresponding to one electron at a large distance 
from the nuclei, the wave function should 
approximate the product of a H-like function for 
the distant electron and an H,*-like function for 
the close electron. The last two conditions 








suggest the use of elliptical coordinates and the 
insertion of an exponential factor reducing to e” 
for either electron at large distances. The neces- 
sary flexibility can best be secured by expressing 
the rest of the function as a power series in the 
five variables, with coefficients to be determined 
by variation to produce the minimum energy. 
Such coefficients can be handled much more 
satisfactorily than parameters appearing in other 
ways (for example, in exponents). Finally, 
Hylleraas,?, working on another two-electron 
problem, the helium atom, has shown the ad- 
vantage of taking the interelectronic distance 
itself, rather than the difference in azimuthal 
angle, as one of the independent variables. 

Let the distances between the several particles 
involved, expressed in terms of the first Bohr 
radius, carry the usual symbols R, fia, 710, fea, 
Yop, 12. Then our coordinates are 


M _ (ria tri) /R, Ae = (Teatro) /R, 


Ma = (‘ia—"w)/R, He = (fea— Yan) /R, 


p= 2ri2/R, 
and our trial function is 
Y= LmnjkpCmnikpLmnjkp ], 
where [mnjkp] stands for 
(1/2) e—8O14¥2) (Ay Ao" Uy!" p? + Ar" Ao” U1 * Uo’ p”). 


The summation is to extend over positive or 
zero values of the indices, subject to the re- 
striction required by nuclear symmetry that 
j+k must be even, and taking as many terms as 
shall prove necessary to give an acceptable ap- 
proximation for the energy. This approximation 
will be found by computing the minimum value of 
SS V¥HvdVdV2 which can be obtained by suit- 
able choice of the coefficients C, subject always 
to the normalizing condition { fy*¥d Vid V.= 1, 
and must always lie above the true energy.* The 
same choice of C’s will then ensure that the trial 
function is as nearly as possible a solution of the 
wave equation Hy= Ey. The errors in the func- 
tion itself will greatly exceed that in the energy; 


2E. A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 
3 See, for instance, J. K. L. MacDonald, Phys. Rev. 43, 
830 (1933). 


4C, Eckart, Phys. Rev. 36, 878 (1930). 





826 H. M. JAMES AND A. S. COOLIDGE 


Eckart* has shown that if the errors are small, 
the former will be of the order of the square root 
of the latter. 

In accordance with what has been said about 
the limiting form of the function, the exponent 
5 should have the value R/2, but additional 
flexibility can be secured without much addi- 
tional calculation by regarding it as another 
arbitrary parameter to be varied. It should be 
noted that the effect of such variation upon the 
wave function will be slight (at least in the 
region where the function is important), provided 
that the coefficients are allowed to re-adjust 
themselves. For a small change in 6 is equivalent 
to multiplication by a rapidly converging power- 
series in \; and do, the effect of which can be 
almost completely absorbed by appropriate 
changes in the coefficients of the smaller powers. 
Now, the whole computation must be done from 
the ground up for each new value of 6, while 
different values of R can be introduced with little 
trouble. Let us limit ourselves to a definite 
number, s, of terms in the series for y. Then 
for each value of 6 there will be one R for which 
the best energy will be lower than it would have 
been for the same R if we had used any other 6. 
To find this, we must plot a series of curves 
(similar in appearance to Morse curves), each 
giving the energy computed with a given 6 for 
various R’s. The envelope of this family of 
curves will then show the lowest energy ob- 
tainable for any R with the corresponding most 
appropriate 6, and each 6 used in computing 
will be the best possible for that particular R at 
which its curve osculates the envelope. (It turns 
out that 6 should somewhat exceed R/2.) 

For a given 6 and R, it can be readily*® shown 
that the values of the C’s which minimize the 
computed energy are those which satisfy the 
system of equations 


(7, — S31) Ci + (2 — ASi2) Co+ ee 
+ (Hi,— Sis) C, =0 


(Hy»— ASi2) Cit (He- Soe) Co+ ae 
+ (Hes — S2.)C,=0 


(Hyp — S14) Cy + (Hos — XSo5) Co + + « 
+ (Iss— S55) Cs = 0. 


































Here \ is a Lagrangian multiplier, and the 
condition of compatibility is the secular equation 


Hy—d\Sy Ty — dSi2 Ths — Siz 
The —_ AS2 To —_ Soe Fo; me So, =(), 
TIis— Si,  $EHes— Ses E55 — S55 








The J/;, and the S;, are the matrix components 
of the Hamiltonian energy operator and of unity, 
respectively, between the fth and gth terms in 
the series, considered as numbered consecutively 
from 1 to s. Details of their computation, of the 
solution of the secular equation, and of the 
determination of the C’s will be reserved for 
the Appendix. As usual, the lowest of the s values 
of X turns out to be the required minimum 
energy, and can be found without determining 
the C’s. 

By way of finding out how many and which 
terms in the series would have to be used in 
order to get a good approximation for the wave 
function, we confined ourselves in the beginning 
to the equilibrium distance, R=1.4, and to 
6=0.75. It soon became apparent that the first 
few terms alone were capable of giving a far 
better energy value than any previously reported, 
and that incorporation of additional terms pro- 
duced only rapidly diminishing improvements. 
In some cases, even terms involving low powers 
of the variables made such slight changes as to 
warrant leaving them out. Of course, the in- 
corporation of each new term required the 
computation of its matrix elements with all the 
terms already in use and the solution of a new 
secular equation of higher degree. Strictly, a 
term should not be neglected until it had been 
tested in combination with all other terms. But 
it soon became clear that the improvement 
obtainable by any given term became progres- 
sively less important as the number of other 
terms present increased, so that we were justified 
in rejecting any term which at any stage in the 
building up of our function was found to produce 
a negligible improvement in the energy. 


RESULTS 


An idea of the results obtained is given in 
Table I, in which the energy at several stages is 
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TABLE I, 

Total Binding 

energy energy Internuclear 

atomic electron- distance, 

Function units volts Bohr radii 

One term — 2.189 — 2.56 1.40 
5 terms — 2.33290 — 4.507 1.40 
11 terms — 2.34580 — 4,682 1.40 
13 terms — 2.34693 — 4.697 1.40 
Without 72 — 2.3154 —4,27 1.40 
Heitler-London —2.21 —2.9 1.40 
Heitler-London —2.24 —3.2 1.51 
Wang —2.278 — 3.76 1.42 
Rosen — 2.297 —4.02 1.416 
Observed —4.73+40.04 1.40 








compared with that given by other methods, and 
with the observed value. The last is obtained by 
adding 0.27 v.e., the zero-point energy of the 
normal molecule, to the heat of dissociation, 
which is 4.46+0.04 according to Richardson and 
Davidson,® while Mulliken® gives 4.44. Much of 
the table is copied from Rosen,’ who, however, 
seems to have neglected the zero-point energy 
in giving the observed value. The values in the 
last column are those for which the various 
functions (other than the series functions here 
presented) give their best results. They are 
sufficiently near R=1.4 for purposes of com- 
parison, except the Heitler-London function, 
whose minimum energy lies at considerably 
greater R, so that it seems worth while to 
include also the energy given by this function at 
the true distance. 











TABLE II. 

Terms Coefficients in normalized functions 

00000] —- 1.69609 2.23779 2.29326  —- 2.22350 
| 00020 | 0.80483 1.19526 1.19279 
[00110 } —0.27997 —0.49921 —0.45805 
| 10000 } —0.60985  —0.86693 —0.82767 
[ 10200 | —0.13656 —0.17134 
F10020} —0.07214 —0.12101 
[10110 ] 0.14330 0.12394 
| 20000 } 0.06621 0.08323 
[00001 } 0.19917 0.33977 0.35076 
r00021 4 0.07090 
[00111 } —0.03143 —0.01143 
| 10001 | — 0.03987 
[00002 | —0.02456 —0.01197 











5 QO. W. Richardson and P. M. Davidson, Proc. Roy. Soc. 
A123, 466 (1929). 

®R.S. Mulliken, Rev. Mod. Phys. 4, 78 (1932). 

7N. Rosen, Phys. Rev. 38, 2099 (1931). 
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The terms used in the first four functions of 
Table I, together with their normalized coeffi- 
cients, are given in Table II. It is interesting to 
note that the simple exponential wave function, 


e~ 8 Qut Ao) =e r1aq/2 «K ew 1b/2  e7T2al2 XK e7 726 2 


already gives a binding energy comparable to 
that of the H—L function; indeed, by a suitable 
choice of 6 the result can be made better, when 
for both functions R=1.4. In the energy thus 
calculated, there is nothing resembling the ‘‘ex- 
change integrals’ of the H—-L treatment; this 
raises the question whether the importance of 
the ‘exchange terms,” frequently assumed to 
represent the essential nature and magnitude of 
chemical binding, may not have been over- 
emphasized. The five-term function is offered as 
a practical compromise between simplicity and 
accuracy. The reason for including both 11- and 
13-term functions is really accidental. Due to a 
mistake, we at first believed that the contribu- 
tions from the terms [00021 ] and [10001] could 
be neglected, and we therefore took the 11 terms 
as a basis for investigating the effects of varying 
5 and R (discussed later). The results showed 
that in fact we had hit upon the best possible 
values of both. In going over the work, we 
discovered the mistake, and found that a slight 
improvement could be obtained by including the 
two given terms, and we therefore offer the 13- 
terms function as the best known approximation 
to the true wave function of He. It did not seem 
worth while to repeat the variations of 6 and R 
with this slightly different form of the function. 

The energy given by the 13-term function is 
within the range fixed by experiment, but higher 
than the most probable value. It is beyond 
doubt that if still more terms were included, the 
result could be still further depressed. We have 
made a rather careful estimate of what could be 
gained in this way, reasoning by analogy from 
the contributions actually found from series of 
related terms; it seems safe to say that we have 
reached the limit of convergence within from 
0.01 to 0.05 volt. The theoretical energy of the 
hydrogen bond is, then, 4.73+0.02 volt, as 
compared with the experimental 4.73+0.04 volt. 

We found it possible to reach the value — 4.27 
with a combination of terms with p=0, and 
therefore not including ri. in the wave function. 
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(The coefficients were not determined in this 
case.) This is substantially the same as the value 
reported by Hylleraas* as the best which he 
could obtain by an involved computation which 
he does not reproduce, and is probably close to 
the limit obtainable with wave functions in four 
coordinates. 

It may be of interest to compare the magni- 
tudes of some of the functions themselves, 
computed for various positions of the electrons, 
and to see how nearly they satisfy Schrédinger’s 
equation. Table III shows values of y and of 


TABLE III. 








Coordinates 





Yia 0.35 0.35 0.35 0.35 1.4 1.4 

Yib 1.4 1.4 1.4 1.4 1.4 1.4 

rea 0.35 0.7 0.7 0.7 1.4 1.4 

r2b 1.75 0.7 1.4 1.4 1.4 1.4 

ri2 0.525 0.7425 0.3556 1.0334 0 2.4249 

¥1-¢ — - 0° 180° 6° 180° 
Values of ¥ (normalized) 

H-L 0.0532 0.0849 0.0533 0.0533 0.0219 0.0219 
Rosen 0.0702 0.1215 0.0710 0.0710 0.0249 0.0249 
5-term 0.0798 0.1172 0.0698 0.0777 0.0162 0.0271 
11-term 0.0783 0.1180 0.0667 0.0777 0.0130 0.0271 
13-term 0.0769 0.1162 0.0656 0.0772 0.0142 0.0269 
Values of (H-E)y, in 27.08 volts 
H-L —0.0425 —0.0481 0.0460 —0.0522 0.0313/p —0.0052 
Rosen —0.0310 —0.0015 0.0879 —0.0429 0.0356/ p —0.0038 
5-term 0.0148 0.0093 0.0653 —0.0097 0.0102/ p —0.0002 
11-term 0.0117 0.0109 0.0356 —0.0060 —0.0034/p —0.0005 
13-term 0.0055 0.0060 0.0269 —0.0021 0.0028/ p —0.0009 








UI—E)y, calculated in each case for R=1.4; 
the unit of energy is here (as in all our compu- 
tations) 27.08 v.e., twice that adopted for 
purposes of comparison in Table I. 

It is natural to regard the discrepancies 
between the various functions as a rough indica- 
tion of their probable errors. The differences 
shown in Table III are representative of a large 
number of points which have been calculated. 
Evidently the H—L function is much too diffuse, 
giving ¥” only about half its proper value in the 
region of importance, and a corresponding ex- 
cessive value in more remote regions in order to 
preserve normalization. Between the other func- 
tions no systematic differences are conspicuous 
except when we compare points in phase-space 
for which the elliptical coordinates of the two 
electrons are respectively the same, but the 
differences in azimuth, and hence the inter- 
electronic separations, are different. Columns 3 


8 E, A. Hylleraas, Zeits. f. Physik 71, 739 (1931). 
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and 4, 5 and 6, offer such comparisons. Here the 
present functions show their characteristic ad- 
vantage of being able to assume smaller values 
for the smaller electronic separations. 

The discrepancies are in general in good 
agreement with Eckart’s' criterion, which can 
be thus stated: the root-mean-square error in y 
is equal to (e/AE)', where e is the error in the 
calculated energy, and AE represents a kind of 
average of the separation between the lowest 
energy level and the other levels of the same 
symmetry. According to Mulliken,® the second 
lowest ',+ state has an energy about 20 volts 
above the ground state (for R=1.4). The state 
of complete ionization of the fixed-center mole- 
cule would evidently have the positive absolute 
energy 19.3 volts (due to nuclear repulsion), or 
51 volts above the ground state. If we take AE 
as 35 volts, e as 1.8, 0.71, 0.22, 0.05, and 0.03 
volts for the five functions, respectively, the 
mean errors should be 23 percent, 14 percent, 
8 percent, 4 percent, and 3 percent. The errors 
in ¥’, and presumably in quantities calculated 
by its aid, will be double these figures. 

Corroboratory evidence is offered by the values 
of JI—E)y. This quantity, which vanishes for 
the correct function, seems on the basis of five 
computed points to tend to amount to something 
like 80 percent, 40 percent, 12 percent, 9 percent, 
6 percent, of y. Reasoning similar to Eckart’s 
(see Appendix) leads to the conclusion that the 
root-mean-square average of these ratios should 
be somewhat greater than (eX AE)?, or 29 per- 
cent, 19 percent, 10 percent, 5 percent, 4 percent, 
respectively. The magnitude of (H—E)y/y 
varies greatly from point to point, but at each 
point there is a clear tendency for the values to 
run proportional to the square-roots of the 
assumed energy errors, thus indicating that those 
errors are indeed of the right order of magnitude, 
and that the correct function would give the 
experimental value. 

It is instructive to see what happens when the 
electrons coincide and r2.=0. Since the term 
1/re in the potential energy becomes infinite, 
Schrédinger’s equation requires either that y 
vanish or that a cancelling infinite term arise 
from the Laplacian or kinetic energy part of the 
operator H. The second alternative is evidently 
the correct one. In the fifth column of Table III, 
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the entries under (J7—F)yW show the coefficient 
of the 1/p term in Jy, which outweighs all 
other terms as p becomes small. In the H—L 
and Rosen functions, this term comes solely 
from the potential energy; in the series functions, 
especially the 13-term function, cancellation by 
a term from the Laplacian is approximately 
complete. The function itself shows no sign of 
approaching zero as it is made more accurate. 
In order to see whether still lower energy 
values could be obtained by varying 6 and/or R, 
we computed a number of energies for different 
values of these quantities, using the same eleven 
terms in the series. As a basis of comparison, we 
took as “‘observed’’ values those given by a 
Morse curve constructed as recommended by 


Mulliken, 
U(r) =D.[1—e-*e) P; a=(27?ucw?/hD,.)}. 


taking D,=4.73 e.v., wo. =4375cm™, r.= 1.40 ay. 
u=0.8309 XK 10-*4 g. 








TABLE IV. 

R 1.2 1.3 1.4 1.5 1.6 1.7 
Binding (6 =0.875 —4.67 -4.63 -—4.52 —4.35 
energy,}6=0.75 —4.41 -—-4.62 -4.68 -—4.63 —4.49 
electron-} Morse —4.49 -—-4.67 -—4.73 —4.68 -—4.57 —4.40 

volts L 








The results are summarized in Table IV. The 
significance of these results will best be seen 
from the graph, Fig. 1. It will be noticed, first, 
that the choice already discussed, 6=0.75, R 
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Fic. 1. Energy of He in electron-volts, against separation 
of the nuclei, in Bohr radii. (Energy of dissociated molecule 
taken as zero.) Lowest curve is “experimental’’ Morse 
curve. Circle indicates point computed with 6=0.75, 11- 
term function. Cross indicates point computed with 
5=0.875, 11-term function. Double circle indicates point 
computed with 6=0.75, 13-term function. 
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TABLE V. 
Coefficients 

PR ae Radiant secede eects sy okt cdr Pah aca set th eG tea ho ee ates 
Terms R=1.2 R=1.3 R=1.4 R=1.5 R=1.4 R=1.5 R=1.6 R=1.7 
[00000 } 2.64085 2.47282 2.29326 2.10831 2.55009 2.42055 2.28454 2.14343 
[00020 | 1.05889 1.14862 1.19526 1.20876 1.32769 1.46843 1.56080 1.61252 
[00110 | — 0.35865 — 0.44227 — 0.49921 — 0.53385 — 0.48136 — 0.59679 — 0.68500 — 0.74964 
| 10000 | — 0.75010 — 0.82703 — 0.86693 — 0.87515 — 0.76638 — 0.84500 — (0.89956 — 0.92919 
[10200 } — 0.07001 — 0.10926 — 0.13656 — 0.15443 — 0.06026 — 0.12001 — 0.16569 — 0.19941 
| 10020 } — 0.00474 — 0.04501 — 0.07214 — 0.08931 +0.03746 — 0.02644 — 0.07409 — 0.10817 
10110 | 0.05723 0.10936 0.14330 0.16388 0.03528 0.11560 0.17625 0.22060 
| 20000 | 0.06082 0.06100 0.06621 0.07168 0.07933 0.07208 0.07421 0.07970 
100001 ] 0.46092 0.39769 0.33977 0.28909 0.54697 0.48577 0.42713 0.37287 
100111 | — 0.01097 — 0.02408 — 0.03143 — 0.03508 — 0.00316 — 0.02137 — 0.03420 — 0.04269 
[00002 } — 0.02616 — 0.02596 — 0.02456 — 0.02264 — 0.03099 — 0.03198 — 0.03135 — 0.02974 











= 1.4, gives the lowest energy of any calculated. 
Second, that the curve for 6=0.75 lies parallel 
to the Morse curve and 0.05 above it at R=1.4, 
while the same is true of the curve for 6=0.875 
for R=1.6. In view of the small effect of changes 
in 6, it seems safe to conclude that, within the 
region studied, and to an accuracy of 0.01 v.e., 
the envelope of a family of curves having various 
6 values would be given by calculating the 
energy for each R with 26=1.08R, and that this 
envelope would lie 0.05 v.e. above the observed 
Morse curve. It follows that the calculated 
equilibrium distance and fundamental vibration 
frequency agree sensibly with the observed. It 
must be admitted, however, that the cogency of 
this reasoning would be much re-enforced by an 
actual computation with 6=0.625, say, in order 
to exclude the conceivable possibility that for 
some R near 1.4 a still lower result might 
come out. 

The coefficients of each term in the series 
were computed for all cases. They vary in a 
regular manner as R and 6 are changed, showing 
no such sharp maxima or minima as the param- 
eters in Rosen’s paper. They are shown in 


Table V. 


OTHER PossIBLE APPLICATIONS 
The excited = states of He which, like the 
normal state, are singlet, with wave functions 
symmetrical in the nuclei, can be attacked by 
the same general method. Hylleraas and Und- 
heim® have shown that the second lowest root 
of the secular equation is an upper limit for the 


®E. A. Hylleraas and B. Undheim, Zeits. f. Physik 65, 
759 (1930). 


energy of the second lowest state of the given 
symmetry, and so on. Thus, by investigating 
the values of the higher roots for various choices 
of R and 6 we may hope to gain approximations 
to the potential curves for the excited states. 
Of course, this may demand the use of more 
terms than are needed for the ground state, or 
even the use of the exponential e~*:~®2% in 
place of e~®™:+%), The various singlet = states 
antisymmetric in the nuclei can be investigated in 
the same way, by using those terms in the series 
with 7+ odd, and the triplet 2 states by using 
the series in which differences instead of sums 
occur: 


y = 2m oe njkpe 002) (Ay™ de "Uy jug* p? 


—Ai"Ao™ U1" uo! p?). 


Once the computation for the singlet states has 
been carried out, the triplet states can be 
treated with very little labor. Some work on 
excited He has been done in this laboratory, and 
will be communicated later. So far, the results 
are disappointing. 

In applying a similar method to molecules 
with but two electrons outside closed shells, 
one would attempt to represent the inner shells 
by the usual atomic functions, the valence 
electrons by a series function. When only terms 
with p=0 are used, the work is straightforward 
and not too difficult. The introduction of the 
terms with 72 involves the evaluation of some 
exceedingly formidable integrals, and it seems 
that some modification of these terms will be 
needed. This is now being investigated in relation 
to the normal states of LiH and Liz. 
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MATHEMATICAL APPENDIX 


Evaluation of matrix components 

Consider the matrix components of unity and of energy arising between any two terms in the 
series from which the wave function is constructed, say the terms [matajakapa] and [mynojokopo ]. 
Each term contains two sub-terms, which we may designate Ya, Ya’, and Ys, ¥s’, respectively; they 
differ by interchange of indices between electrons. The matrix component of unity will have the form 


a f f dV id Vabedo + f f dVid Vb! + ‘ f d Vid Vabe!o+ | f dVid Vai’. 


It will be sufficient to consider the first of these integrals, which we may call s,». Four similarly 
related integrals will form the matrix component of the energy; the first of these will be 


hab = { d Vid Vowallpr 


-{ dVid Vowal. —3Vr- 2V2+(1/R) [14+2/p—4r1/ (Ax? — wn”) —4)o/ (Ao? — ps”) } Wo. 


Owing to the fact that the complete terms ya+ya’ and ¥,+y»’ contain the coordinates of the 
two electrons symmetrically, we may use instead the simpler integrals 


Yom f f d Vid Vebol — Vat + (1/R) {1+2/p—8Y0/ (ds? — nn?) } Wo. 


Then 
Hay=h' aa th'a'y +h’ av’ +h’ o's. 


In performing these integrations, we were not able to utilize the device adopted by Hylleraas for 
He, using p directly as one of the variables of integration. We found it necessary to transform the 
powers of p into elliptical coordinates, by means of the equations 


p=AP +e? wi? + Me” — 2 — 2Ardopipe — 2C (Ar — 1) (Ae? — 1) (1 — mi?) (1 — ys") } cos (¢1— ¢2), 


p=> peP.()or()P. (m1) P,” (ue) cos v(gi— ¢2), (notation of Zener and Guillemin"’) 
T =0 2 1 


=0 v 


D£=2r+1; D,=(—1)’X2(2r+1)[(7—v)!/(7 +)! for v>0. 


Now, the required quantities may all be found in terms of integrals defined thus: 


1 
X(m, n, 7]; k, p) ~ Ag? Lf ff (Aq? — py?) 78D) yo” Uy o* p?AA ded Md pod gid vo. 
TT 


The indices are all zero or positive, except that p may have the value —1. 
Remembering that d Vid V2= (1/64) R®( Aq? — py”) (Ao? — po”) d Ard Aod wd pod gid geo, we find at once 


Sab = (1/64) R°LX (mat my, na tnot+2, jatjo, Ratko, Pat po) 


—X(mMatmy, Nato, jatjo, Ratkot+2, pat ps) ]. 
We shall abbreviate this to 


Sav = (1/64) R°LX (02000) —.X (00020) J. 


0 C, Zener and V. Guillemin, Jr., Phys. Rev. 34, 999 (1929). 
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In hav’, the only troublesome part is that arising from the Laplacian. Insofar as the coordinates 
of electron 2 are concerned, we may write y as constant X L5(Az2) M/5(u2)Ps(p). Then 
V2"Vo/Wo = V2 Le/Lot+ V2? Me/Mot+V2Ps/Pot2VeLe:VePs/LoPot+2V2Ms:-VeP1/MsP2; 
(since the term in V2l,-V2M,z always vanishes, the coordinates being orthogonal). Now, in general, 
V°f(p) =L4/R’p? ] X L2pdf/dp+ p*d*f/dp*]; 
VF(A, mM) = L4/R2(? — w?) |X [2d Of/OX+ (M?— 1) d2f/ON? — 2p Of/du+ (1 —w?)0*f/du"]; 
VA(A, H) VA, Hy ¢) = L4/R2(M?— yw?) |X L(A? — 1) (Af/Ad) (0g /AA) + (1 — wv?) (Of/Ou) (Og/du) J. 


Putting 
Ly = eho", My = ue", Py=p, 
we find 
[ (Ae? — 1) /Ae JOP s/Od2 = 5 poP opp? — Av>+ do? — wr? — po? + 2Arwre/Ae 
C(1 — me”) /u2 JOP /Op2= — > poP op? p? — dr? — do? — wr? + mo? + 2Ardous / Ue ] 
and finally 


WVaV2Wo = (4WaWo/R? (Ao? — wo”) |X {2(4 —Sd2) + (1 — 1/do?) (124(4 — 1) — 2642+ 62d?) 
— 2kot(1/pe?— 1) (Ro(Ro—1)) + (po/p”) L (bo +1) (Ao? — me”) 
+ (ny— 5X2) (p? — +e? — pr? — po? + 2AiMibe, ne) — kop? — Ax? — do? — wr? + oe? + 2r1AoH1/ M2) 1 } . 


This can evidently be integrated in terms of the functions X previously defined. But a somewhat 
shorter expression for the integral can be found. By Greene’s theorem, for properly continuous 
functions which vanish canonically at infinity, 


f dV aV2bo= 4 f dV 2(PaV bot WoV2Wa— 2V 2a" Var). 
We note 
Voba Vado / Pals = Vala Vale/LaLo+V2Ma-V2Mo/MaMot+V2P a: V2Po/PaPs+V2la:V2Po/LaPs 
+V2Ma:VoPs/MaPot+Vols:V2Pa/LePatVeM oe: V2Po/ MPa. 
When this is evaluated, some cancellation occurs, and the final formula is 
h’ 6 = (1/64) R® {| X (02000) — X (00020) + 2X (0200 — 1) — 2X (0002 — 1) —8X (01000) } 
— (1/64) R4{ 0 (ta— ms)? — (Ra— ko)? + (atmo) — (Ratko) + (Pa po) (Ma— Ns — ka tke) |X (00000) 
—465X (01000) — [(a— 1s)? — (nats) |X (0— 2000) + [(ka— ke)? — (Ra the) |X (000 — 20) 
+[(pa— po)? + pat pot (pa— po) (Na—nvo+ka— ke) JLX (0200 —2) —X (0002 — 2) ] 
— (pa— po) (Ma — Nv — (Ra — ke) ) LX (2000 — 2) +X (0020 — 2) ] 
+2(pa— pv) (Ma—ms)X (1—111—2) —2(pa— po) (Ra — kv) X (111—1—2)}. 


(The same abbreviation as before has been made in the argument of the X’s.) 
In discussing the computation of these integrals, it will be expedient to introduce the abbreviation 
a= 26. Let us define the more general function 


X’(m, n, J, k, p) =Z’(m+2, n, j, k, p) —Z"(m, n, j +2, k, p); 


Z’(m, n, j, k, p) =(1/4") LS Sf [ecomorereuinator Me cost (1— ex)ddudadud ad end e' 


M? = (Ay? — 1) (Ae? — 1) (1 — mn?) (1 — we?) 
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This evidently reduces to the previously defined X(m, n, j, k, p) when v=0. 
For p=0, we have 
Z(m, n, j, k, 0) =4Am(a)An(a)/(G+1)(R+1), 
Z'(m, n, j, k, 0) =0, 
Z"(m, n, j, k, 0) =8[Am42(@) — Am(a) LA ngo(a) — A nla) J/(G-+1)(G+3) (R41) (k+3), 


when j and k are both even; otherwise the integrals vanish. (Here, as usual, A ,(a)=f"X"e~™ dX. 
For details of computation and tables, see Rosen’s paper.) 


For p= —1, the Neumann expansion of 1/p yields 
Z(m,n,j,k, —1)= me (27+1)R,(7)R-(k)H,(m, n, a), 
7=0 


Z'(m, n, j,k, —1)=-¥ [1/r(7+1)JR,’(a) Re’ (k) X[(27 +1) /7(7 +1) JH’ (m, n, a), 
t=1 


Z" (m,n, j,k, —1) =4 3 [1/(27+3)(27—1)(r+2)(7+1)7(7—1) JR,(9) Re’ (R) 
 [ (27 +3) (27 +1) (27 —1) /(7 +2) (7 +1)7(7 —1) JH,’ (m, n, a) 


4255 (274+1)(R,(j) —R-(j+2)) (Re(k) —R-(k+2)) 


T=0 


< L7,(m+2, n+2, a) —H,(m+2, n, a) —H,(m, n+2, a) +H,(m, n, a) ]. 


Here, R,’(j) = Stidu(1 —p?)”*P,”(u)u?. Rosen has discussed and tabulated this integral for v=0. 
The other cases are readily obtained from the elementary theory of associated harmonics. They 
exist only for even or for odd 7, according as 7+ is even or odd, and vanish for all values of 7 greater 
than j+¥, so that the summation over 7 contains at most a few terms. 

The other integral 


ao co Ay de 
IT,’(m, n, a) = anf Adel (Ar? — 1) (Ao? — 1 Pe ( Jor( Jarmrare-svrns 
neal | 1 Ne Ai 


is familiar in molecular problems, at least for »=0. Rosen gives a method for computing F/o(m, n, a) 
(his formula A19). With the aid of his function S(m, n, a) (formulas A13, A14), and the following 
recurrence relations, other values can be found readily: 


IT,(m, n, a) =Ho(m+1, n+1, a) —S(m, n+1, a) —S(n, m+1, a). 

IT,(m, n, a) = (1/7?) [(27 —1)?27,_1(m +1, n +1, a) +(7—1)*H,_2(m, n, a) 
— (2r—1) (27 —3) (H,_2(m4+2, n, a) +H,-2(m, n+2, a))+2(27—1)(27 —5) H,-3(m+1, n+1, a) 
— (27—1)(27—7) (H,_4(m42, n, aw) +H ,_4(m, n+2, w))+2(27—1)(27-—9) H,_5(m+1, n +1, a) 
—-++] _— either 

(for even r) —(27—1)[Ho(m4+2, n, a) +Ho(m, n+2, a) —S(m+1, n, a) —S(n+1, m, a) ] 

(for odd +r) +(2r—1)[2Ho(m+1, n+1, a) —S(m, n+1, a) —S(n, m+1, a). 

[(27+1)/r(r+1) JH,’ (m, n, a) =(r+1)H,41(m, n, aw) —(27+1)H,(m+1, n+1, aw) +7H,-1(m, n, a). 


for r>1. 





"'N. Rosen, Phys. Rev. 38, 255 (1931). 
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[(27+3)(2r+1)(27—1)/(r+2)(7+1)7(7—1) JH," (m, n, a) 
= 7(27—1)[(27 +3) /(7 +1) (7 +2) JH’ -41(m, n, a) 
—(2r+3)(2r—1)[(27+1)/7(7 +1) JH,’ (m+1,n+1, a) +(7+1)(27+3)[(27—1)/ (7-1) 7 J’ (m, n, a). 


(The factors in brackets [ ] before 7’ and H” are left in explicitly because for purposes of compu- 
tation they may conveniently be considered part of the functions.) 

When X” has been tabulated for p=0 and p= —1, it can easily be found for positive values of p 
by means of the relation (with the usual abbreviation) 


X”(00002) = X(20000) +.X*(02000) +.X*(00200) +.X"(00020) — 2X»(00000) 
—2X°(11110) —2X*+4(00000). 


A corresponding relation holds also for Z”. 


Solution of secular equation, and determination of constants 

The complexity of the secular equation of course necessitates determination of the roots by trial 
and error. The following process seems to be the simplest available, and puts the work in a particu- 
larly satisfactory form for the determination of the constants. Choosing a trial value of \ in the 
proper neighborhood (known from experiment, or experience with a similar function), we may write 
the s equations which must be satisfied by the coefficients C in the form 


~4:-C, =0, 4=1,2,--: 5, dip = Hi, — Siz. 


t=1 


This set of equations is equivalent to another set, obtained from it by progressive elimination of 
the variables C,, in which there is one equation involving s of the C’s, one in (s—1)C’s, and so on. 
Let these equations be 


> D:-C, =0, 4=1,2,---s. 
t=1 


The D’s are conveniently found by the recurrence relation 


n—1 


Das =dn-— | 
v=1 


starting with D,,=d,,. Because the d’s form a symmetrical matrix, it will be found that D,, vanishes 
automatically for 7<m. A convenient numerical check is provided by the fact that the quantities 
~,D,, should be obtainable from the quantities 2,d,, by means of the same recurrence relation. 

If \ has been properly chosen, the last factor D,, will vanish, and the last equation D,,C,=0 will 
permit C, to be set arbitrarily equal to 1 (or any other number). Since the secular determinant is 
easily seen to be just the product of the diagonal factors D;;, the value of \ so determined is evi- 
dently that which makes the determinant vanish. The remaining C’s may be immediately determined 
by successively substituting C, in the (s—1)st equation, C, and C,_; in the (s—2)nd, and so on. 
The provisional coefficients are then used to compute { fdVidV2¥’; upon division by the square 
root of this quantity, a set of normalized coefficients is obtained. 

Should D,,, vanish for n<s, the meaning is that the value of \ chosen happens to be an exact 
solution of the secular equation containing only the first m terms. In this case (and also in practise 
if Dn» is very small), the procedure for subsequent D’s breaks down. If, nevertheless, it is desired 
to find the value of the whole determinant, the order of the original terms must be changed, placing 
at least one of the originally first » terms in a position later than the nth. 
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In general, of course, the chosen value of \ will prove incorrect, the last factor, and therefore the 
whole determinant, failing to vanish. In order to find a value of \ making D,, negligibly small, the 
above process is carried out for two assumed values of \ lying rather near the root, and preferably 
enclosing it, as shown by the last factor having opposite signs in the two cases. The root can then 
be quite accurately determined by considering the value of the determinant (not the last factor) as 
linear in X, and interpolating. If the coefficients are desired, they must be obtained by a new compu- 
tation using the interpolated 4, as it is not safe to try to interpolate the individual D’s which deter- 
mine the coefficients. 

It may be mentioned that by this method an 11-row determinant with six significant figures can 
be evaluated and checked by an experienced computer in a little over two hours. 


Order of magnitude of (H—E)y 


Suppose our approximate wave function y expanded in terms of the correct wave functions y, 
of proper symmetry, with coefficients a,, such that a; is nearly unity, and 2,a,?=1. Let E be the 
energy calculated, and £, the correct energy for the nth function. Then 


E= [vine v~ f (Ea4.) (Sa..v.)a V=)aZE, =aF,+ >07E,; 
T o T T=2 


=ark, +EDa? =aE,+E(1 —a,’), 


T=2 


where E is some value between E2 and E,,. Similarly we can show that 
2 = [Cur-mype V =a;(E,—E)?+ (E-E)*(1 —a;"), E,<E<E,,. 
Eliminating 1—a,’ gives &=(£,—E)[E,—E+} (E—E)?— (E,—E)*} /(E,—E)]. 


Assuming that (E—E)/(E—E,) is of the order of unity, and neglecting higher powers of e=E—F, 
we have 6?~(E,—E)(E,—E). This is a lower limit, since it can easily be shown that E>E. 
Note added in proof: 


A minute error has been found, affecting the energy by 0.003 e.v. Corrected coefficients will 
be published as soon as possible. 
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Dielectric Constant Studies. I. An Improved Voltage Tuning Resonance Method and 
Its Application to Aqueous Potassium Chloride Solutions 
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(Received August 30, 1933) 


An experimental study of the behavior of various types 
of apparatus used in the measurement of dielectric con- 
stants of aqueous solutions of electrolytes. Although simple 
theory indicates that the voltage tuning method of detec- 
tion should give results which are uninfluenced by the 
conductivity of the solutions and should, hence, be par- 
ticularly adapted to measurements on electrolytic solu- 


tions, it appears, however, that other disturbing factors are 
present. This is shown by the divergent results obtained 
when different types of detection are employed and leads 
to the conclusion that the change in dielectric constant of 
such solutions as observed by this method depends more 
or less upon the characteristics of the particular generating 
and receiving circuits used. 





HE experimental evidence for the effect of 
ions on the dielectric constant of liquids is 
very conflicting. Debye,! after a consideration of 
the results, obtained prior to 1929, states that, 
“no definite conclusion with regard to any 
theory can be drawn from measurements in their 
present condition.’’ Since 1928, numerous papers 
have appeared and one might expect considerable 
improvement; however, the pronounced lack of 
agreement in published data forces the conclusion 
that this is not the case. The seriousness of this 
situation was emphasized to the authors in the 
course of some work on the dielectric constant of 
gelatin solutions, the results of which will be 
reported in a later paper. In connection with 
that work, it became evident that in some 
solutions conductivity had an important bearing 
on the results. These circumstances led the 
authors to undertake a series of investigations 
in an attempt to clarify the situation in respect 
to aqueous solutions of electrolytes. This aim 
has necessitated the use of several methods of 
measurement in order to contrast directly their 
advantages and limitations. 
It is the purpose of this paper to discuss the 
so-called voltage tuning resonance method used 
in various forms and with varying success by 


*This paper constitutes a portion of the dissertation 
submitted by J. Gilbert Malone to the faculty of the 
Graduate School, University of Michigan, in partial fulfill- 
ment of the requirements for the degree of Sc.D. 

1 Debye, Polar Molecules, p. 124, Chem. Cat. Co., N. Y., 
1929. 


Jezewski,’? Kniepkamp,’ Lattey and co-workers,‘ 
and Astin.® Since the presentation of this work 
in the form of a dissertation, another paper by 
Jezewski® has appeared which describes work 
similar in its purpose to that presented here. 
Since the theory has been given by Jezewski, 
and can also be found in the papers of Lattey 
already referred to, it need be described but 
briefly here. The voltage tuning principle derives 
its significance from the fact that analysis of the 
tunable portion of the circuit of Fig. 1 shows that, 
if the value of C is adjusted to give a maximum 
voltage across the condenser, this maximum will 
be independent of the resistance x in parallel 
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Fic. 1. Diagrammatic resonating circuit. 








2 Jezewski, Zeits. f. Physik 48, 123 (1928). 

3 Kniepkamp, Zeits. f. Physik 51, 95 (1928). 

* Lattey (et. al.), Phil. Mag. 7, 985 (1929); 11, 997 (1931); 
12, 1111 (1931); 13, 444 (1932). 

5 Astin, Phys. Rev. 34, 300 (1929). 

6 Jezewski, Phys. Zeits. 34, 88 (1933). 
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with it. This analysis assumes an inductance 
without capacitance, and a condenser without 
inductance. Moreover, there must be no reaction 
of the resonance circuit on the oscillator. 


THE OSCILLATOR 


The oscillations were produced by the use of 
a UX210 tube in the circuit shown in Fig. 2. 


C=250 ap 


+ ~400Volts _ 


L— 


C 
0.025 pF 














~5 Volts 











Fic. 2. Oscillator. 


The filament was lighted and the plate supplied 
with the proper voltages by a transformer unit 
especially designed for this purpose. The in- 
ductances of the plate circuit and the grid 
circuit were wound on the same Bakelite core 
and fitted with plug connections facilitating their 
interchange. UX bases are convenient for this 
purpose as they will not permit the coil so 
mounted to be placed wrongly in the circuit. 
This oscillator functioned well for wave-lengths 
from 20 to 100 meters. 


THE RESONATING CIRCUIT 


When the voltage tuning resonance method is 
employed there are several factors to be con- 
sidered which do not appear from the simple 
theory. The simple theory takes no account of 
the inductance of the cell leads. It is obvious, 
however, from consideration of the circuit shown 
in Fig. 1 that the inductance will be changed 
because of the inductance in the leads represented 
by a and a’ if the total capacity is transferred 
from the cell K to the precision condenser C. 
It is evident also that some sort of a meter must 


be used to indicate the maximum voltage across 
the condenser C, and whatever sort is used, it 
presents complications by introducing capaci- 
tance, inductance, resistance or some combina- 
tion of these. Lattey and his co-workers employ 
a thermocouple galvanometer connected across 
the condenser terminals through a capacitive 
impedance. They believe this arrangement to be 
superior to all others for the purpose. 

Two distinct resonating circuits have been 
used in this work. The first followed that of 
Lattey and Davies’ as closely as possible, and is 
shown in Fig. 3. This circuit as constructed and 


L 




















KI] SKU _|K*0.0005mF 








Fic. 3. First receiver. 


used by the present authors, gave results in 
qualitative agreement with those obtained by 
Lattey and Davies. Since the second circuit is 
believed to represent an advance, and since it 
possesses far greater sensitivity and does not 
sacrifice or compromise any principle embodied 
in the first, no further discussion of the first 
circuit is considered necessary. 

The second circuit is shown in Fig. 4. In 
principle it is an alternating voltage indicating 
device. In this arrangement it is not necessary 
to balance out an initial plate current as is 
usually done; the plate current is initially zero 
because of the d.c. grid bias. When the resonance 
condition is approached an increase in the plate 
current of the tube results and can be measured 
by a milliammeter in the plate circuit. When the 
voltage across the condenser is a maximum, the 
plate current is likewise a maximum, and both 
of these conditions occur, as has been pointed 
out by Jezewski, when w*=1/LC. 

It is desirable that the parts of the resonator 








7 Lattey and Davies, Phil. Mag. 12, 1111 (1931). 
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Fic. 4. Second receiver. 


be arranged as compactly and as symmetrically 
as possible. The inductances were made from 
8X2 mm copper strip and were mounted on 
short extensions of the condenser terminals. The 
tube connections were placed at about the same 
point, and provision was made to allow the cell 
to be inserted directly across the terminals of 
the condenser and the tube connections. It was 
possible to bring all of these connections together 
by inclining the condenser to the extent that a 
radio tube could be hung just below its terminals. 
The extension of these terminals provided a 
means for securing the inductance and the cell. 

The condenser used was General Radio type 
222-L. A UX222 tube was employed. All leads 
and the tube were shielded and grounded. The 
dielectric cells used are illustrated by the one 
shown in Fig. 5. These were made of glass and 
platinum with very short leads. It has been 
concluded by Lattey that in such a case the 
apparent dielectric capacity will not be affected 
by the conductivity of the solution. 

The oscillator was normally kept at a distance 
of about 90 cm from the resonator. Tests showed 
that effects of coupling were observable at 30 cm. 
It has been observed also that unless one works 
at a considerable distance from conductors such 
as water pipes, etc., effects due to these may be 
attributed to coupling. On the basis of such 
observations the authors believe that the chang- 
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Fic. 5. Dielectric constant cell. 


ing of the relative positions of oscillator and 
resonator is undesirable. 

The setting of the condenser, giving maximum 
current, which is the condition of resonance, is 
difficult to determine with great precision and is 
often obtained, in the case of wave meters, by 
the insertion of a small capacitance in parallel 
with the variable condenser. Any setting of the 
condenser such that insertion of this capacitance 
makes no change in the reading of the indicating 
device requires that the capacitance has spanned 
the resonance point to give a corresponding 
point on the other side. This procedure, while 
admirable for wave meters, was not very success- 
ful here and an alternative was sought. 

It is generally known that the curve obtained 
when current is plotted against condenser read- 
ings is symmetrical if the calibration curve of the 
condenser is linear. To what extent this is true 
in the present case may be seen from the figures 
in Table I. 


TABLE I. Data on symmetry of resonance curve. 








Condenser Values 





Microamperes Above Below Mean 
400 5390 5260 5325 
300 5431 5230 5330 
200 5470 5172 5329 
150 5520 5122 5321 
100 5744 4864 5304 
Cell filled with solution 
110 5104 5502 5303 
140 5278 5322 5300 








The first column gives the current as shown 
by the ammeter A in microamperes, the next 
two columns give the two condenser settings on 
opposite sides of the resonance point correspond- 
ing to the currents in the first column. The last 
column shows the mean of these condenser 
values. If the above noted symmetrical relation- 
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ship holds, the values in the last column should 
remain constant. The extent of the agreement 
indicates that the symmetrical relationship is 
closely approximated until the curve is about 
400 condenser divisions in width. The values 
obtained with an electrolyte show that the curve 
is flattened as a result of electrolytic conduction. 
Only a limited range may, therefore, be followed 
and the accuracy is correspondingly decreased. 

If the calibration curve of the condenser were 
linear, a number of readings might be made at 
corresponding current values on either side of 
the resonance point, and the resonance point 
obtained as above by taking the mean. In 
practice it was found necessary to use the lower 
settings of the condenser where the _ linear 
relationship does not hold and in such cases 
three procedures are possible. In the first method, 
the condenser may be calibrated and the readings 
converted in terms of the calibration. This 
method was not employed in the present case, 
since experience indicated two other methods 
which were more convenient and of equal accu- 
racy. In both the average of the two condenser 
readings obtained with the same current but on 
opposite branches of the resonance curve was 
plotted against the difference between these 
readings, i.e., against the width of the resonance 
curve. The two methods differed, however, in 
the fact that in method II the plate current 
was always adjusted to the same value while in 
method III it was allowed to vary. 

The procedure used in taking measurements 
according to method II was as follows: Two 
condenser readings were taken, one on each side 
of the resonance point for the same plate current. 
The condenser was then set to a new value and, 
after the plate current had been restored to the 
original value by suitable variation of the grid 
bias, another condenser reading was taken on 
the opposite side of the resonance point. This 
procedure was then repeated for as many points 
as desired. 

In method III, on the other hand, the plate 
Current was allowed to vary with the closeness 
of approach to the resonance point, although 
again for each PAIR of readings the plate current 
was of course the same. 

These two methods gave satisfactory agree- 


ment over the node width of 200 to 600 divisions, 
which range included practically all of the 
measurements in this work. In practice, however, 
method III was preferred in that it gave some- 
what more consistent results. 

In applying the latter method to the measure- 
ment of the dielectric constant of an electrolytic 
solution, curves were determined for water and 
for the solution in the manner just described and 
both plotted on the same graph as illustrated in 
Fig. 6. The curves obtained approximate straight 
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Fic. 6. Plots of typical observational series for the 
dielectric constant of KCI solutions by the voltage tuning 
method. 


lines and the curve for the solution always falls 
below the one for water. The mean ordinate 
distance between these curves is the increase in 
the capacity of the cell when filled with solution 
over its capacity when filled with water. The di- 
electric constant of the solution is calculated from 
this value. Following this method of procedure, 
the troublesome conversion of condenser readings 
into comparable units is avoided and at the same 
time the effect of an asymmetric wave is mini- 
mized. 

The apparatus and cells described were de- 
signed for the measurement of the dielectric 
constants of aqueous solutions of electrolytes. 
It was assumed that the relationship between 
cell capacity and dielectric constant is linear, 
and that benzene and water can be employed as 
calibrating materials to determine this line. 
Several pure materials, covering a wide range of 
dielectric constants, were measured. The values 
obtained agreed well with those recorded in the 
literature. 






















































840 MALONE, 
THE DIELECTRIC CONSTANT OF AQUEOUS 
SOLUTIONS OF KCL 


A review of this subject has been given by 
Lattey and Davies.’ Many observers have meas- 
ured a large number of solutions.* The conflicting 
nature of the results obtained emphasizes the 
futility of making extensive measurements on 
solutions of a large number of carefully purified 
inorganic substances at the present time. It has 
been the purpose of the present authors to confine 
measurements to a few substances and to study 
the possibilities of several methods. The Drude 
method in several modified forms has been used 
by us, and is the subject of a separate communi- 
cation.* 

Lattey and his co-workers have used the 
voltage tuning resonance method most exten- 
sively. Some of his methods, on first considera- 
tion, appear different from and more elegant 
than those employed by us. Careful analysis, 
however, reveals that this is probably not the 
case. Let us consider the problem of finding the 
true resonance point. The impressed voltage 
under proper conditions, as has been indicated, 
gives a potential across the condenser which is 
symmetrical with respect to capacity and falls 
off with capacity values either less or greater 
than this. In order to obtain Ko, the resonance 
capacity, Lattey considers that two values for 
the capacity are obtained, one on either side of 
the resonance capacity Ko, which give the corre- 
sponding currents J; and J,’ and a current J at 
the resonance point. He is thus enabled to write: 


(Ko—Ky)/(0?/I?—1)'=¢ 
and 
(Ki! —K;)/(Iet/!*—1)'= 6. 


Eliminating Ko: 
o = (Ky' — Ki) /(Ue?/I? — 1)§+ Ie? /h? — 1)! 


¢ is the tuning factor. 

Lattey and Davies determine current values 
with the condenser set at definite values and 
calculate Ko by means of the preceding equations. 
It is obvious from the derivation of the equation 


8 Compare for instance the results of Milicka and Slama, 
Ann. d. Physik (5) 8, 663 (1931) with those of Lattey and 
Davies, Phil. Mag. 12, 1111 (1931); 13, 444 (1932). 

* This communication is really number II in this series 
and follows the present article. 
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that, if J’;=J;, the values of Ko will be the 
arithmetical mean of K; and K;’. The method 
followed in our work consists in reading the 
values of capacitance, first with the current at a 
definite value and then after the current has 
been brought to the same value on the other side 
of the resonance point by a proper adjustment 
of the capacity. Such a procedure is more 
convenient in observation and calculation than 
that used by Lattey and Davies and is more 
fundamentally sound in that the accuracy of the 
indicating device does not introduce a source of 
error. This is of special significance when a 
thermionic tube is employed as the voltmeter to 
measure the potential across the condenser. 

The measurements of Lattey and Davies show 
that the precision was not high. Thus, with a 
solution of approximately 0.001 N, they show a 
variation of from 1.5 to 3.0 percent and with a 
0.12 N solution the variation is about 10 percent 
in the extreme cases. Nevertheless, the trend of 
the values reported by Lattey is definite; and 
the first consideration in this work is to obtain 
results which, even though qualitative in nature, 
are fundamentally sound. 

Astin’s results’ with KCl by the resonance 
method are represented by three curves, the first 
showing a decrease to a minimum at a concen- 
tration of about 0.001 N followed by a rapid 
increase to values above that of water, while 
the other two, with closer coupling, show a 
continued decrease, more rapid with a cylindrical 
plate condenser than with a parallel plate 
condenser. 

The data obtained by the present authors do 
not confirm those of Astin but show instead a 
continued rise in the dielectric constant with 
concentration. The observed data are recorded 
in Table II. These values were taken at 25°C at 
a wave-length of 31.7 meters. Cell A had a 
capacity of 606 uwuf and cell B a capacity of 
28.2 uuf when filled with water. 

The value for pure water has been accepted as 
the reference value. Only low concentrations 
could be measured with the larger cell due to 
conductivity. The apparent dielectric constant 
change due to solute is small yet definite in 
direction, fairly reproducible and increases with 
increasing concentration. 

A concentration of 0.0025 N solution gave 
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TABLE II.* Summary of dielectric constants (uncorrected) of 
aqueous KCl solutions. 








Cell A Cell B 
Conc. 1 2 1 2 

0.0000 | 78.58 78.58 | 78.58 78.58 78.58 78.58 
0.0010 | 78.13 79.13 | 79.45 79.42 79.42 79.17 


3 4 








0.0020 79.69 | 79.93 80.95 80.27 80.44 
0.0025 81.21 | 80.27 81.97 81.12 
0.0050 81.63 85.71 














* Note: Each value in each column is the average of 
series of measurements made by refilling the cell with the 
same solution. The results in the different columns were 
made with freshly prepared solutions, at different times 
and under somewhat different experimental conditions, in 
order to determine whether there were any uncontrolled 
factors which might be having an effect. 


dielectric constants of 81.57, 81.35 and 81.42 at 
wave-lengths of 24.5, 46.5 and 82.4 meters, 
respectively. These values seem to indicate that 
the results are independent of frequency. 

The values here recorded for KCl are in 
agreement essentially with those obtained by 
Lattey and Davies. They are in contradiction 
to the values generaliy obtained with the Drude 
method. 

It has been pointed out that the theoretical 
treatment assumes a pure inductance. Every 
conductor, however, has a certain amount of 
capacitance associated with it. The authors 
attempted to change the distributed capacitance 
effect by immersing the inductance coil of the 
measuring apparatus in benzene. Results ob- 
tained on salt solutions with the coil immersed 
gave, within the experimental error, no change 
of dielectric constant with concentration of the 
salt. The dielectric constant of all solutions, 


within the experimental error, under these con- 
ditions, was that of pure water. The benzene 
was not replaced with materials of higher di- 
electric constant for the reason that conductivity 
would be a factor in the use of materials such 
as the alcohols and the amount of material 
required and the experimental conditions made 
the use of volatile or expensive materials im- 
practical. Even if the curve of apparent dielectric 
constants against the dielectric constant of the 
milieu surrounding the inductance is established 
the extrapolation to zero dielectric constant, i.e., 
to zero distributed capacity, is probably not 
justifiable. It was not practical to surround all 
of the inductance with benzene, and the problem 
may be one of reflection effects at the surface of 
the benzene. This circumstance, while it renders 
the problem even more complicated, does not 
change the basic fact that the distributed 
capacity of the coil is a factor of importance in 
using the voltage tuning resonance method and, 
as yet, has not been satisfactorily taken into 
account. 

The results given are consistent with the 
values obtained by Jezewski as well as qualita- 
tively in agreement with those of Wien.’ In 
contradiction to Jezewski the present authors 
are convinced, however, that, even though their 
data recorded in Table II indicate a fairly 
definite behavior for aqueous KCI solutions, yet, 
when all things are considered, the voltage 
tuning resonance method cannot be depended 
upon to distinguish between the effects of in- 
creased conductivity and dielectric constant 
change in aqueous salt solutions. 


® Wien, Ann. d. Physik 11, 429 (1931). 
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Dielectric Constant Studies. II. The Drude Method Applied to Aqueous Solutions of 
Potassium Chloride 


J. GILBERT MALONE,* L. O. CAsE AND A. L. FERGUSON, Chemistry Laboratory, University of Michigan 
(Received August 30, 1933) 


This paper is an experimental comparison of several 
adaptations of the Drude method as applied to conducting 
solutions. With three distinct types of detecting device, 
results for the dielectric constant of water and of water- 
alcohol mixtures agreed within experimental error. When, 
however, the same apparatus was applied to measurements 


with KCI solutions, different results were obtained with 
each. This indicates that any apparent change in the dielec- 
tric constant of electrolytic solutions, as observed by the 
Drude method, must depend in whole or in part upon the 
characteristics of the apparatus. 





ITH the intention of measuring the dielec- 

tric constants of some electrolytic solu- 

tions, the authors constructed several modifica- 

tions of the Drude apparatus. Comparison of the 

results obtained with these modifications, by 

using the same electrolyte, indicated a pro- 

nounced ambiguity, which, judging from our own 

experience and from work reported in the litera- 

ture, seems to be characteristic, also, of the results 
of other methods. 

Of the methods which have been employed for 
the measurement of dielectric constant of leaky 
dielectrics, the Drude method has much in its 
favor. It has, however, three obvious limitations: 
first, it requires a large amount of material; sec- 
ond, the material used must be in the liquid 
state; and third, because of the length of the tube 
required it can be operated only at very high 
frequencies. On the other hand, it appears to 
possess the best theoretical foundation, since the 
dielectric constant of a leaky dielectric is defined 
directly by the Maxwell equation: 


(eu/c?)0?E /d?+ (4rpp/c)dE/dt=V°E, 


where p is the conductivity of the medium, « its 
dielectric constant, w its magnetic permeability, 
c the velocity of light and E the electric intensity. 
The method has been applied to electrolytic solu- 
tions by several investigators, among whom may 


* Rewritten from a thesis presented by J. Gilbert Malone 
in partial fulfillment of the requirements for the degree of 
Doctor of Science at the University of Michigan. 
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be mentioned Drude,! Deubner,? Devoto,* Drake, 
Pierce and Dow.‘ In general the results with this 
method were interpreted as indicating practically 
no change of dielectric constant with concentra- 
tion of the electrolyte. (An exception might be 
noted in the result of Deubner? for CuSO,; he 
found a decrease in the dielectric constant with 
an increase in concentration, a result in disagree- 
ment with the earlier observations of Drude on 
the same material.) This, it will be noted, is at 
variance with the results of other methods, some 
of which’ seem to indicate an increase and others*® 
a decrease in the dielectric constant with concen- 
tration of electrolyte. 

The essentials of the Drude apparatus are a 
source of high frequency oscillations, a Lecher 
wire system and a means of detecting the nodal or 
antinodal points on the wires. 

The high frequency oscillations used by 
Drude! and other earlier workers consisted of 
damped waves. Deubner,? and Drake, Pierce and 
Dow‘ employed oscillations from the usual re- 


1 Drude, Zeits. f. physik. Chemie 23, 267 (1897). 

2 Deubner, Ann. d. Physik 84, 429 (1927). 

3’ Devoto, Gazz. chim. ital. 60, 208 (1930). 

4 Drake, Pierce and Dow, Phys. Rev. 35, 613 (1930). 

5 Wien, Ann. d. Physik (5) 11, 429 (1931). Lattey and 
co-workers, Phil. Mag. 7, 985 (1929); 11, 997 (1931); 12, 
1111 (1931); 13, 441 (1932). Jezewski, Phys. Zeits. 34, 88 
(1933). 

® Walden, Ulich and Werner, Zeits. f. physik. Chemie 
129, 389 (1927). Shanke and Schreiner, Phys. Zeits. 28, 597 
(1927). Sack, Phys. Zeits. 28, 199 (1927). Bramley, J. 
Frank. Inst. 205, 649 (1928). 
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generative thermionic oscillators, while Devoto* 
used oscillations from a Barkhausen-Kurtz 
generator. 

The earlier form of the Lecher wire system used 
by Drude and by Deubner consisted of the con- 
ventional arrangement of two parallel wires sur- 
rounded by the medium to be investigated. 
Devoto*® modified the Drude method by the use 
of a single wire instead of two. Drake, Pierce and 
Dow? substituted a concentric cylinder system 
for the parallel wires. 

The detecting device used by Drude was a neon 
bulb, Devoto and Deubner employed a milliam- 
meter of the thermocouple type, while Drake, 
Pierce and Dow used the reaction of the Lecher 
wire system on the oscillator to give the nodal 
points. 

The present authors have employed two types 
of Lecher wire systems, three distinct types of 
detection and a source of undamped waves. 


EXPERIMENTAL 


A diagram of the apparatus is shown in Fig. 1. 
For measurements of the wave-length in air, 
two wires, (see Fig. la) W and W’, 1.5 inches 
apart, were mounted vertically between the floor 
and ceiling of the room, the ends being secured 
by Bakelite supports. The lower ends were con- 
nected to the terminals of a Weston milliam- 
meter while a shorting bridge B, constructed of 
metal and Bakelite and provided with spring 
contacts, was movable on the wire system. The 
oscillator was placed on the floor within a foot or 
two of the wire system, and the positions of the 
nodal points of the waves induced on the wires 
were determined by moving the bridge and not- 
ing the points of maximum current as shown by 
the milliammeter. 

For the measurement of the wave-length in 
liquid media a similar arrangement (Fig. 1c) was 
employed except that in this case the wires were 
» inch apart and were surrounded for the greater 
part of their length by a glass tube 5 inches in 
diameter and about 7 feet long. The lower end of 
the tube was closed by a rubber stopper, through 
which the wires protruded to the milliammeter. 
The connections of the wires at the upper end of 
the tube are shown in the figure. In order to 
maintain a constant temperature, the liquid was 
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(a) (b) ¢) 
Lecher Pierce type Lecher System 
System System in glass tube 
in’ air 


Fic. 1. Diagram of apparatus. 


circulated by means of an electrically driven 
pump connected with glass and rubber tubing to 
the top and bottom of the large glass tube. Heat- 
ing units were located around the tube and the 
whole was heat insulated. A thermometer in one 
of the glass connections from the pump gave a 
convenient means of reading the temperature. 
The shorting bridge B providing reflection of the 
wave within the tube was moved from the top by 
means of a calibrated Bakelite rod (C), which, by 
means of a vernier, was read to a precision of 
0.01 cm. The ammeter at the bottom was read 
from the top by means of a telescope. 

In addition to the glass tube apparatus just 
described the system of Drake, Pierce and Dow 
was also employed. This apparatus (Fig. 1b) 
consisted of a brass tube 7.5 feet long and 2.5 
inches in diameter, closed at the lower end by a 
rubber stopper and containing an axially located 
phosphorbronze wire D. The circulation of the 
liquid and the detecting arrangement used here 
were of the same type as described above for the 
glass tube. 






































All wires, the interior of the metal tube, and 
the exposed parts of the pump were heavily silver 
plated, and where parts were exposed to con- 
siderable wear solid silver was used. In the first 
detecting arrangement employed in the present 
work, the ammeter was fixed at the generator end 
of the wire system and the shorting bridge moved 
to locate the maxima. This arrangement has been 
used by Laville.’ 

The undamped waves employed were gener- 
ated by a vacuum tube oscillator. The measure- 
ments were carried out in the following manner. 
The wave-length in air was first determined with 
the auxiliary wire system. In the early part of the 
work these measurements were checked on each 
of the other systems, but it soon became evident 
that this was not necessary. Without changing 
the generator the wave-length in the liquid was 
then determined with either the glass or brass 
tube arrangement. The wave-length in air was 
remeasured at the end of each experiment to 
make sure that it had not changed. The ratio of 
the half wave-lengths thus obtained was taken as 
the index of refraction, ~. For nonconductive 
media n?=e. 

With the object of comparing the two tubes 
and of testing the accuracy of the apparatus, 
some preliminary measurements were made with 
water and with water-alcohol mixtures. The re- 
sults with water are shown in Table I. It will be 
noted that the glass tube gave slightly lower 
values for e than the concentric cylinder system. 
This result was attributed to the fact that a 
sensible portion of the field lay outside the water 
in the tube. Support was given to this explana- 
tion by the still lower value of 72.90 obtained 
when the Lecher wires were further apart. 

The results obtained with alcohol-water mix- 
tures covering the entire concentration range 
agreed satisfactorily with those of Wyman* the 
average deviation being 0.75 percent. 

It seemed that the above determinations es- 
tablished the reliability of the apparatus for 
work with non-electrolytes and accordingly some 
work on electrolytes was attempted. For these 
measurements, KCI and CuSO,, representing two 
types of ions, were selected, solutions of these 

7 Laville, Les Ondes Electriques Courtes, p. 126, R. Mesny, 


Paris (1927). 
8’ Wyman, J. Am. Chem. Soc. 53, 3292 (1931). 
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TABLE I. Comparison of parallel wire and concentric tube 











systems. 
/2 d/2 € € 
Temp. air in liquid at T at 25°C 
20.4 176.2 g 19.74 79.68 78.03 
m 19.67 81.25 78.59 
16.6 176.2 g 19.55 81.23 78.61 
water changed m 19.54 81.31 78.69 
21.2 176.2 g19.74 © 79.68 78.36 
m 19.74 79.68 78.36 
24.2 176.2 g 19.86 78.71 78.36 
m 19.85 78.79 78.54 
29.0 176.2 g 20.06 77.15 78.59 
water changed m 20.04 77.31 78.75 
23.3 176.1 g 19.85 78.70 78.09 
m 19,77 79.34 78.73 
24.0 176.2 g 19.92 78.24 77.88 
m 19.87 78.64 78.28 
23.2 176.2 g 19.86 78.71 78.11 
m 19.81 79.11 78.51 
30.4 176.2 g 20.15 76.46 78.40 
m 20.13 76.61 78.55 
water changed 
20.1 176.2 g 19.69 80.08 78.22 
m 19.67 80.24 78.48 
22.6 176.2 g 19.78 79.35 78.49 
m 19.77 79.43 78.57 
26.5 176.2 g 19.96 77.93 78.47 
m 19.95 78.00 78.68 
water changed 
23.6 176.2 g 19.86 78.71 78.20 
m 19.80 79.19 78.69 
26.0 176.2 g 19.94 78.01 78.43 
m 19,93 78.16 78.52 
25.0 176.2 g 19.92 78.24 78.24 
m 19.86 78.71 78.71 








Average as measured by glass tube 78.36. 

Average as measured by metal tube 78.58 at 24°C. 
m is value measured by metal tube system. 

g refers to the glass container and the Lecher wire system. 


salts having been frequently investigated by 
others. The results obtained are given in Fig. 2. 
The values given are for ?; since in this work we 
are concerned only with relative values, no at- 
tempt has been made to apply a correction for 
conductivity, which would be uncertain at best. 
It will be noted that the values show a decided 
decrease with increase in concentration, whereas 
others, working with the Drude method, as was 
pointed out earlier, have found practically no 
effect on the dielectric constant due to electro- 
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lytes. These discrepancies may be due at least in 
part to differences in the equipment used. For 
instance, the various detecting devices no doubt 
introduce complicating factors in the application 
of simple theory, factors which may produce re- 
sults differing both in magnitude and direction 
from the correct nm? value. No quantitative 
theoretical relationship taking these factors into 
account has so far been developed. 

An experimental study of this question was 
undertaken by the use of two additional types 
of detecting device. The first was a voltmeter 
constructed from a UX222 tube and which had 
been found satisfactory in connection with some 
resonance work being carried on in this labora- 
tory. The potentials across the input end of the 
apparatus could be followed by the changes in 
the plate current of the tube. The grid was con- 
nected to the center wire, forming the inner 
cylinder of the concentric cylinder system, while 
the filament was connected to the tube by a 
short heavy wire and was also grounded. This 
apparatus functioned very well but it was more 
sensitive to changes in the surroundings than was 
the case with the previous detector. With this 
device, the values obtained for water and for one 
concentration of ethyl alcohol, 34.5 percent, were 
consistent with those previously obtained. While 
in these cases the position of the nodes was al- 
tered by the substitution of the new detector, the 
distance between nodes was the same as before. 
However, the precision was much less than in 
the case already described. With electrolytes, 
however, greatly different results were obtained. 





Although the results were not very reproducible, 
a large and fairly uniform increase of n® with in- 
creasing concentration was observed. 

A third detecting device was used, depending 
upon the reaction of the measuring system upon 
an additional resonant circuit, the changes in 
which were observed by means of the vacuum 
tube voltmeter which was inserted across the in- 
put terminals. With this arrangement the nodes 
in the measuring circuit could be observed by 
noting the points of maximum reaction on the 
added circuit as shown by the plate current of 
the vacuum tube in it. Here again the values ob- 
served for alcohol and water were consistent with 
those obtained in the preceding cases, but the 
values observed for KCl solutions were quite 
different. The n? values calculated from the ob- 
servations show again an increase with increas- 
ing concentration of electrolyte, but less than in 
the preceding case. The limited extent of this 
increase indicates that there is little change in 
dielectric constant when it is calculated from n? 
by means of the Maxwell relationship for a leaky 
dielectric. Such a result approximates that ob- 
tained by Drake, Pierce and Dow using a similar 
method. 

The last two procedures were not carried to 
high refinement by the present authors and al- 
though they gave fairly consistent results they 
were particularly sensitive to changes of capaci- 
tance in the surroundings. 

It is generally realized that the problem of the 
dielectric constant of electrolytic solutions is 
highly complex. This was the conclusion arrived 
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at by Debye® in 1929 as a result of a critical 
study of the literature up to that time. As stated 
in the first paper of this series, the authors under- 
took this work in an attempt to clarify this situa- 
tion. The results show, however, that the prob- 
lem is even more complex than has been realized. 
The measurements reported in this and the pre- 
ceding paper indicate that while self-consistent 
results can be obtained even for electrolytes by 
each of the various methods in the hands of differ- 
ent observers, there is no conformity in the re- 
sults of different methods, even when employed 


9 Debye, Polar Molecules, Chapter VI, Chemical Catalog 
Co., New York (1929). 
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by the same investigator.'’” Since, as has been 
pointed out, these discrepancies do not appear in 
the case of non-electrolytic solutions, it is evident 
that conductivity is responsible for the anomalous 
results obtained, and that its effect is different 
with different methods. The few attempts which 
have been made to treat this effect theoretically 
are inadequate, and until an adequate treatment 
has been made the changes observed by any 
method cannot be regarded as representing true 
changes in dielectric constant. 


10 cf, Lattey, Phil. Mag. 41, 829 (1921); 12, 1111 (1931); 
13, 444 (1932). 
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The Heat of Formation of Binary Liquid Solutions from Their Liquid Components 


J. J. COLEMAN AND FRANK E. E. GERMANN, Depariment of Chemistry, University of Colorado 
(Received October 9, 1933) 


It is pointed out that an equation giving the relation of 
the heat of formation per mol of liquid solutions from their 
liquid components, to temperature and composition can 
be used as the foundation equation for their thermodynam- 
ical theory. The empirical equation 

AH = N,(1—N))(a+pe™), 
in which AH is the heat of formation per mol of the solu- 
tion, N; is the mol fraction of one of the components of the 


binary solution, and a, B and y are constants or functions 
of temperature, is presented. There are two degenerate 


forms of this equation. (Eqs. (7) and (8).) There are listed 
values of a, 8B and y (or the values of the corresponding 
quantities in Eqs. (7) and (8)) for 59 solutions for one tem- 
perature and sometimes for several temperatures. It is 
noted that according to the data assembled in the Jnterna- 
tional Critical Tables, some solutions composed of ‘‘unas- 
sociated” liquids present a behavior more irregular than 
those composed of ‘‘associated”’ liquids. It is suggested that 
this unexpected behavior may mean that some of these 
data are not as precise as they are indicated to be. 





|. USE OF THE HEAT OF FORMATION EQUATION AS 
THE FOUNDATION EQUATION FOR THE COM- 
PLETE THERMODYNAMICAL THEORY OF 
THE SOLUTION 


T has long been known that many of the prob- 
lems of solutions can be reduced to that of 
finding the dependence of their thermodynamical 
properties upon temperature and composition. 
This has been done for dilute solutions, for 
“ideal’”’ solutions, for van der Waals solutions? 
(i.e., those in which the van der Waals equation 
can be used as the equation of state of the com- 
ponents of the solutions in both the liquid and 
gaseous state) and for one or two other types of 
solutions. Only the theory of van der Waals solu- 
tions has any general application and this yields 
only qualitative results. 

We cannot get all of these relations from 
thermodynamics alone. We must have at least 
one of them to start with. In the case of ‘‘ideal’”’ 
solutions we can take as the fundamental extra- 
thermodynamical relation 


fi=Nift (1) 


where f; is the fugacity of any component in the 





‘Lewis and Randall, Thermodynamics, Chapters XIX 
and XX, 

*van der Waals and Kohnstamm, Lehrbuch der Thermo- 
Statik, Part II. 
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solution, N, is its mol fraction, and f;° is its 
fugacity in a pure phase of the same kind (solid, 
liquid or gaseous) as that of the solution. 

Now by the methods of K. Watanabe* we 
can derive the equation, valid for any component 
in any solution 


H,-H,° IT 


In fi=ln Nifi?— { RT? ’ (2) 


(P, Ni) 


in which i; is the partial molal heat content of 
the component in the solution, //,° is its molal 
heat content in the pure phase, 7 is the tem- 
perature and R is the gas constant. The integra- 
tion is indicated to be at constant pressure and 
composition. This integral is evidently the cor- 
rection term necessary to make (1) (written in 
the logarithmic form) applicable to non-ideal 
solutions. And the Hi—H,° appearing in the 
integrand can be related to the heat of formation 
of the solution by differentiation of the obvious 
equation 


TH, —nylTy® — nell2® = AH (ny + ne) 


(in which JJ, is the heat content of a solution 
containing ”, mols of component one and m2 mols 
of component two and A// is the heat of forma- 
tion per mol of the solution) with respect to 


3 Watanabe, On the Vapour Pressure of Liquid, Tokyo 
Science Reports, Sect. A, No. 6 (1931). 
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n(T, p and nz constant). This gives 








H,—H, 

AAT 
> On, a T, p 
rd ATTq ON, 

= (+m) +4H (3) 
LON, 4n2, T, p On, ne, 1’, p 
TOA 

=(1-N,)|——| +4H. 
5 ON, i 7. 2 








(There is a similar equation for solutions of more 
than two components but we are concerned with 
binary solutions only in this article.) 

It is evident then that Eq. (2) enables us to 
use any relation giving the dependence of AJ/ 
on N, and T as the fundamental extra-thermo- 
dynamical equation for non-ideal solutions. But 
before turning to such equations we will give a 
derivation of (2). 

For the change 


dn, (pure liquid) 
+(,+7n2) (solution) ->(:+2+dm) (solution) 


we have a change of free energy and of heat 


content 
dn(Fi— F,), dn,(H,—H,°). 


These two quantities are related by an equation 
from the general theory of thermodynamics 


Ni, P T? 








ih~ F,°)dny (H, —H,")dn 
(aa). Sm 


Here dn, the amount dissolved, is, of course, 
independent of 7. Cancelling this out and inte- 
grating we get (after dividing by R) 

Fi-FY fi 

AN tl 


RT fi 


in which g(M,, p) is an arbitrary function. 
Consider now the solution under any pressure 
p and any temperature T. If we first increase the 


Ih-HY 
—— +g(Mi, p), (4) 


 P) 
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pressure until it exceeds the critical pressure of 
the solution and then increase the temperature 
until it exceeds the critical temperature, we will 
have passed by a series of continuous states from 
a liquid solution to a gaseous solution. It seems 
safe to assume that the form of the function 
g(Ni, p) is not altered in the first step of this 
process and that it is not changed at all by the 
second step; and since the gaseous solution is 
ideal or can be made so by a further increase in 
temperature we can use the relations for ideal 
solutions 


H,—H,=0, fi=Nify’, 


which give us 
g=In M. 


Substitution into (4) then gives the desired 
Eq. (2). 


Il. THE Baup-HEITLER EQUATION 
In 1915 Baud‘ brought forth the equation 
AH=bN,(1-— WN), (5) 


in which 0} is a constant, and showed that his 
equation served for a number of solutions. Later 
Heitler® gave a derivation of this equation based 
upon statistical mechanics and gave some of the 
thermodynamical relations for solutions for 
which this equation holds. Hildebrand® also used 
this equation in developing the thermodynamics 
of a new class of solutions which he has called 
“regular solutions.”’ It is probably the most serv- 
iceable equation for the heat of formation of 
binary solutions that has been previously 
developed. 


Ill. A New Heat OF FORMATION EQUATION 


It is evident then that the study of the func- 
tion AH(T, N,) offers a promising approach to 
the problem of extending the thermodynamical 
theory of solutions. We have sought first to find 
an equation for the relation between AH and .\; 
which would serve for solutions in which the 


4 Baud, Analyse Thermique des Melanges Binaires, Bull. 
Soc. Chim. (4) 17, 329-45 (1915). 

5 Heitler, Ann. d. Physik (4) 80, 630 (1926). 

6 Hildebrand, J. Am. Chem. Soc. 51, 66 (1929). 
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Baud-Heitler one fails. A purely empirical analy- 
sis of the results assembled in the International 
Critical Tables resulted in the following equation: 


AH = N,(1—N;)(a+-Be™!), (6) 


where a, 8 and y are either constants or functions 
of temperature. This formula was found to have 
a wide application, but before turning to this we 
will consider some interesting special cases. 

If a and B are of opposite sign then it may hap- 
pen that AH changes in sign as we proceed over 
the composition range. This is a frequent phe- 
nomenon. If y is small enough to enable us to 
use the approximation e*=1+x then Eq. (6) 
reduces to 


AH = N,(1—N:)(ai1 +6811), (7) 
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in which a; and £; are new constants or functions 
of temperature. If either 8 or y is zero then the 
equation reduces to the Baud-Heitler form 


AH = Ni(1—Mi)az. (8) 


Finally if a, 8 and y are all zero we have, of 
course, the case of the “‘ideal’’ solution. 

In Table I we have listed values for a, 8 and y 
(or a, and #1, or a2 as the case may be) for a 
number of solutions. In the first four columns we 
give these and the temperature (or interval with- 
in which the temperature lies). In columns five to 
eight we give some data which are intended to 
give some notion of how well the equation serves. 
In column five is listed the maximum value of a 
deviation index which is the numerical value of 


AH given by the equation —AH given by experiment 





x 100. 


Average of experimental values (all considered positive) 


We did not use the percentage deviation for if we 
had sought to make this a minimum, it would 
have meant giving the greater weight to observa- 
tions of smaller numerical value, where the rela- 
tive precision is lower. We assumed that, unless 
otherwise indicated, the absolute precision was 
the same for all measurements. Hence we sought 
to keep the absolute deviations as small as pos- 
sible and our index indicates how well we were 
able to do this. In some cases all of the data ex- 
cept one could be represented nicely by our 
equation. In these cases we ignored this one result 
in the determination of a, 8 and y and the index 
of deviation. Hence we have given in columns six 
and seven the number of observations ignored 
and the total number of them. (In this last we in- 
cluded the two zero values of AH at N,=1 and 
N,=0.) In some cases while what data there were 
could be easily represented by our equation, their 
distribution over the composition range was so 
poor as to make uncertain the applicability of 
the equation or, assuming this, the values for a, 
8 and y given. Hence in column eight we have 
commented on the distribution of the data. 

In Figs. 1 and 2 we give two illustrations. In 
Fig. 1 we have plotted the heat of formation of 
CS: - CioHig (pinene) solutions against the mol 
fraction of CS». The solid line is the plot of the 








Fic. 1. Variation of heat of formation of CS.—C,Hi¢ solu- 
tions with mol fraction of CSs. 
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Fic. 2. Variation of heat of formation of C,H,Cl.—C;Hs 
solutions with mol fraction of C2H4Cle. 
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TABLE I. 
(The heat unit used is the kilojoule.) 
» Q z 
& 3 a % S 2 3 c | 8 
“2 S65 |25| § 2 5Si5 |Se!| 8 
g . BSlcs~is5) S | s 2 : Baleezlss| S | s 
S sie ESEE\=8| 2 |e 3 s |a| |geSbees! 2 | 
F >|] s BSlez else) 2 | é > | &§ | 5 | |BS/S8el82| 3 | 3 
5 £ a e ~leeizerlee] A 1S fe) q y elem moe") & 15 
H20 (Component one) C2HiCle Ethylene chloride (Component one) 
CH,O Methyl C7Hs Toluene | 15-20 |— 690] 111 } —] 4] 1 | 9 |Good| 5 
> bel 0 — 2.95 |— .0525/ 5.00) 3 0 21 | Good! 1 
4 ethy . r . 
coalgghel | 19.69 |~ 2.85 |— .108 4.00] 9 | 0 | 21 | Good | 4 CEE hnkine® Compe ane 
4 ethy Ae on al abe om ’ 
cileohol, [4287 | 1.88 |— 284 |8.00 3 | 0 | 21 |Good} 1 CHOEgicther [2% |-2 | — 1-1 31 © | 8 Heed] 3 
aeabal 0 |— 1.75 |— .038 |6.06} 5 | 0 | 21|Good| 1 C2H.O: Acetic acid (Component one) 
C:HeO2Glycol [32 |— 1.40 |— .232 |3.00| 7 | 1 | 11|Good| 2 
C2H¢02 Glycol 55 |— 1.27 |— .213 | 3.00} 2 1 | 11]|Good| 2 |}CcéH::Cyclohexane 15-20; 6.77 | — |—] 4| 1 9 |Good| 5 
C2Hs02 Glycol 76 |— 1.27 |— .191 |3.00} 9 | 0 | 111] Good| 2 ||C7Hs Toluene ¢/ 15-20; 1409 | — |—]| 3] 0 8 | Fair | 5 
H2O Water ? 2.25 |—.376| 2.00} 9 0 11 | Fair 9 
CCl, (Component one) 
CoHBee Ethy! C2Hc0 Ethyl aleohol (Component one) 
2HaBre Ethylene 
Bromide 25 1.92 31 —|5 0 9 | Good| 3 |/CsHeO Acet 25 4.63 i buat © 0 10 | Good! 3 
pe ered 25 44 — —|7 0 8 _— _ — mf | | | | | 
7N Aniline 25 4.67 01 | 6.20) 3 0 13 | Goo 4 ; 
C;Hs Toluene 17 1061 :106 | —|7 | 0 | 11|Good| 6 CARD Anstons (Compuneat one 
: CHO Methylalcohol | 25 | 2.32 82|—| 2| 0 10 | Good| 3 
CS (Component one) CiHs02 Ethyl acetate |25 | .544| — | oe E | 0 | 9|Good| : 
CCl4 Carb a 
otrachloride 25 1.14 0239 | 3.00] 2 0 11 |Good| 3 Cs3HeO2 Propionic acid (Component one) 
3 Chloroform | 25 1.93 48 —j| 2 0 11 | Good| 3 
a mid H20 Water } 8 | 120 | 7.90 |—| 4] O | 13|Good| 9 
romide 25 2.52 .222 | 2.00) 1 0 8 00 3 
CsHs02 Ethyl C3HsO2 Methyl acetate (Component one) 
acetate 25 3.95 .159 | 2.70) 1 0 9|Good| 3 
CsHi00 Ethyl ether | 25 — 1.52 |— .63 —|3 0 9|Good| 3 ||CiHsO2Ethylacetate |]16 | .279| — |—I 3] 1 | 11 | Good, 6 
— ame 4 1.77 .0531 = 3 0 11 vee 4 
6He Benzene 14.5 1.73 .055 | 3.00) 2 0 10 | Goo 6 » » Ethy 
CoH Bensene, 95 1193 67 oe Fe 0 8 |Good| 3 CiHs02 Ethyl acetate (Component one) 
'6H1203 Paralde- . - ° : ; 
CioHic Pinene 25 1.10 .0790 | 3.00) 5 0 9 | Good 3 i sheokel sate 25 5.88 .0106| 6.00 4 ‘ » Good 3 
thy] ether 25 1.04 — _ 700 3 
CHCl; (Component one) CHO Isobutyl 
y 
CCl, Carbon alcohol 25 7.43 — |-—| 4 0 9 | Good 
cli Ethyl 25 eee ae sat fs 0 8 | Good} 3 CsH100 Ethyl ether (Component one) 
] — 6. , P 
C:HsO: Ethyl és ee . : ome e ae ate - 25 57 1248 |—]| 4 0 9 | Fair 3 
acetate 25 — 5.77 |—4.55 —|5 0 10 | Fair 3 ||VatisU n-rropy “ 
CsHioOEthylether}25  |-105 | — |—| 4] 0 | 10|Good| 3 |}, alcohol 25 1.94 | .189/ 3.00} 4) 0 | 10 | Good 
CiH100 Ethylether|14 |— 9.96} — |—]|7 | 1 | 11| Good] 6 ||C+Hi00 Isobutyl 4 
CcHsC! Chloro- alcohol + 25 2.27 | .417/2.00} 1) 0 | 10|Good| 3 
benzene — et P, CsH120 Isoamy] alcoho) |25 1.95 .397 | 2.00} 6 0 10 | Good} 3 
15-20 598 .210 1 0 6 ‘oor 5 
CeéHe Benzene 15-20|— 1.00 |— .75 ane 2 0 8 |Good| 5 ||CeH:203Paraldehyde [25 1.03 45 |—| 5 0 9 | Good| 3 
CeHe Benzene 25 — 1.40 |— .70 —|8 0 8 |Good| 3 
CeHioCyclohexene | 15-20)}— 443) — —| 3 0 5 | Poor | 5 CeH;Cl Chlorobenzene (Component one) 
CeHi2Cyclohexane | 15-20 2.63 — —| 3 0 7 | Fair 5 . 
CeH120s Paralde- CsHi0 p-Xylene [15-20/— 48 | — | —] 5] 0 { 7|Goodj 5 
hyde 25 — 6.05 |—4.65 —| 5 0 9|Good| 3 
CsHi0 p-Xylene 15-20 |— 3.23 |— .80 —j1 0 8|Good| 5 CsHs (Component one) 
CHA Methyl sleobol (Component one) CH.O Methyl alcohol | 15 1.10 | .240/3.00, 10] 0 | 11| Good] 6 
C2H:O Ethyl C2Hs02 Acetic acid ? 1.73 47 |—]|] 9] O 15 | Good | 12 
alcohol 3 053} ~— nae 0 5|Fair | 8 ome encom 2° = —|- . : 7 ost 5 
J , : His n-Hexane 5- 86 —|— air | 5 
coe 3.00 | 0.146 {3.001 3] 0 | 5| Fair | 3 ||C7Hs Toluene 15-20) 208) — |—| 7] 0 | 8|Good) 3 
Calfso0 Ethyl ether 25 57 | 2.48 eae) 0 9|Good| 3 C;H: Toluene 16 298} — | —| 10 0 11 | Good| 6 
m 
pm te " 25 22 | 1.08 ee 0 5 | Fair 3 CéH1:Cyclohexane (Component one) 
C2HiBrz Ethylene bromide (Component one) C7Hs Toluene | 15-20] 243 | 47 |—|10 0 | 9|Good| 5 
C2HiO2 Acetic acid ? 2.20 0722 | 3.00} 2 0 10 | Good | 12 CsHio p-Xylene (Component one) 
CeHs Benzene 15-20 1.16 — —|3 0 9|Good} 5 
CcH12Cyclohexane | 15-20 5.56 -- —|3 0 14 | Good} 5 |/CeHzN Aniline 10 | —3.80 |—.150 |3.00} 3] 0 {| 8| Poor! 11 
OBSERVERS 
1 Bose, Zeits. f. physik. Chemie 58, 585 (1907). 8 Bose, Ann. de Chim. et de Phys. 27, 89 (1912). 
2 Schwers, Rec. trav. Chim. 28, 42 (1909). ® Faucon, Ann. de Chim. et de Phys. 19, 70 (1910). 
’ Hirobe, J. Faculty of Science (Tokyo) 1, 155 (1926). 10 Drucker and Moles, Zeits. f. physik. Chemie 75, 405 
‘ Hartung, Trans. Faraday Soc. 12, 66 (1917). (1911). 
5 Baud, Bull. Soc. Chim. de France 17, 329 (1915). 11 Clark, Phys. Zeits. 6, 154 (1905). 
* Schmidt, Zeits. f. physik. Chemie 121, 211 (1926). 12 Baud, Ann. de Chim. et de Phys. 27, 89 (1912). 
7 Winkelmann, Ann. d. Physik 30, 529 (1873). 
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equation 
AH = N,(1 — Nj) (1.10+0.0790e**:), 
the dotted line is the plot of the equation 


AH = N,(1— N,)1.64 


and the circles represent the experimental data. 
The dotted line indicates that in this case a 
Baud-Heitler equation will serve as a rough ap- 
proximation. The maximum value of the devia- 
tion index (for the first equation) ignoring none 
of the experimental results is five. In Fig. 2 we 
give asimilar graph for CeH,Clz—C;Hs (toluene) 
solutions. The equation is 


AH = Ni(1—Ni)(—0.690+1.11). 


The maximum deviation index, ignoring one of 
the experimental results, is four. The result ig- 
nored (AH =0.661, N,=0.385) cannot be shown 
on the graph. The Baud-Heitler equation cannot, 
of course, be applied to solutions of this type. 
There are a large number of cases in which the 
equation cannot be made to represent the data 
and it is quite evident that it is not perfectly 
general. However, solutions composed of such 





notoriously bad behavers as water and the alco- 
hols gave points on a plot of AH/N,(1—™i) 
against N, which could be represented by a 
simple smooth curve and so it was expected that 
solutions of simpler (i.e., ‘‘unassociated’’) liquids 
would do the same. But these in many cases gave 
very irregular plots. This unexpected behavior 
suggests that it might be well to check some of 
these data. 

There are not sufficient data available to en- 
able us to attack the remaining part of the prob- 
lem, the determination of the dependence of a, 
B and y on 7, in this direct empirical fashion 
without some preliminary experimental work. 
We could proceed at once to deduce thermo- 
dynamical relations for solutions for which 
Eq. (6) serves, leaving the a, 8 and ¥ in as un- 
determined functions of temperature. But the 
resulting equations are so clumsy that their ap- 
plication would be extremely difficult. Hence it 
seems best to delay this until we have further in- 
formation concerning a, 8 and y. One looks to 
statistical mechanics for an interpretation of this 
equation and such investigations when successful 
will, no doubt, throw some light on these 
quantities. 
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Contact Potentials and the Effects of Unimolecular Films on Surface Potentials. I. 
Films of Acids and Alcohols* 


Wi1Lu1AM D. HARKINs AND E. K. FiscHEer,t Department of Chemistry, University of Chicago 
(Received August 19, 1933) 


Simultaneous measurements of film pressure and surface 
potential have been carried out by the use of an apparatus 
designed in such a way that the potential may be deter- 
mined for any location on the surface of the film. At film 
pressures above that of the gaseous films, organic sub- 
stances with homo-heteropolar molecules give a single 
smooth curve for the relation between surface potential 
and molecular area. At areas sufficiently great to reduce the 
pressure to that of the gaseous film the surface potential 
becomes variable and remains variable until the area be- 
comes so great that the continents and islands of con- 
densed film evaporate in the two-dimensional system to 
give a gaseous film alone. For example with films of myristic 
acid at 17° the surface potential is represented by a single 
curve below a molecular area of about 50 sq. A, and by any 
value below 170 mv at higher areas, at which islands in the 
film persist. The areas above which the surface potentials 
become variable, due to the effects of islands or continents 
of film, are found to be about as follows: stearic acid, 28; 
palmitic acid, 28; pentadecylic acid, 41; lauric acid (on 
NaCl-HCl solution), 36; oleic acid, 52; cetyl alcohol, 23.7; 
myristyl alcohol, 27. The surface potential rises rapidly 
with the number of carbon atoms in the hydrocarbon chain, 
and is about 400 mv for condensed films of stearic acid. 
It is higher for alcohols than for acids with the same num- 
ber of carbon atoms. The increment is about that given by 
the addition of two carbon atoms to the acid. The change 
of one of the single bonds of stearic acid to the double bond 


of oleic acid as the film forming substance lowers the sur- 
face potential by from 37 to 65 percent, and makes this 
potential much more greatly affected by the molecular area. 
The following numbers give the number of carbon atoms in 
several normal acids and the corresponding potentials at 
20 sq. A at 17°: (12), 276; (14), 336; (15), 375; (16) 383; 
and (18), 410. The surface potential increases with decrease 
of temperature, with an increase of concentration in the 
film, and for the substances investigated, by a replacement 
of a basic by an acid solution. A collapse of the film lowers 
the surface potential since it changes the molecular orienta- 
tion in the film. The films exhibit in the condensed state 
considerable hysteresis, both with respect to film pressure 
and surface potential. This feature has not been sufficiently 
noted in other work. The relation between surface potential 
and the molecular moments of the molecules which consti- 
tute the film is discussed. The dipole moment as calculated 
from surface potentials by the equation of Helmholtz are 
only about one-sixth as large as those found by ordinary 
methods for independent molecules in nonpolar solvents. 
The discrepancy is attributed partly to a lowering of poten- 
tial by an orientation of the dipoles with respect to each 
other, and partly to the too great simplicity of the equation. 
The apparatus used for this work has been found to give 
good determinations of the contact potentials of metals. In 
its design use was made of ideas suggested principally by 
Lord Kelvin, and by Guyot and by Frumkin. 





I. INTRODUCTION 


T has been supposed that the surface of a pure 
liquid or of a solution is the seat of several 
different types of potentials. Of these the greatest 
prominence in the theory of the stability of col- 
loids, as developed by Hardy and later by others, 
has been given to the electrokinetic or zeta (¢) 
potential, the original theory of which was de- 
veloped by Helmholtz. It has not been found 


* The work reported in this article was carried out with 
the aid of a fellowship from the Julius Stieglitz Fund for 
Research in Chemistry Applied to Medicine, established 
at the University of Chicago by the Chemical Foundation. 
The objective of the present work is a study of membrane 
formation. 

t Julius Stieglitz Fellow, University of Chicago. 


possible to measure the quantity designated by ¢ 
directly as a potential, and it is necessary to in- 
troduce relative motion of the phases along the 
interface if an electromotive force (£) is to be 
developed, and this appears as a potential differ- 
ence along and not through the interface. For this 
reason ¢ is sometimes considered as a tangential 
potential, although it is E, which may have an 
entirely different order of magnitude, which is 
measure! tar sentially. 

Wha e considered as the contact poten- 
tial of th: surface is less fictitious, since it may be 
determined by a direct measurement. If the sur- 
faces of two metals are separated by a gas, and if 
the gas is sufficiently ionized, a potential, known 
as the contact potential (V) is developed between 
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the two metals. This particular method for the 
determination of contact potentials was de- 
veloped by Kelvin.! 

What may be designated as the surface poten- 
tial (AV) of a solution was determined in 1883 by 
Bichat and Blondlot? and later by other workers.* 
The method used by Kenrick was to allow the 
reference liquid to flow down over the inner wall 
of a moderately large glass tube placed with its 
axis vertical, and to allow the second liquid to 
flow from a fine glass tip in a narrow jet down this 
axis. The two solutions were connected through 
calomel electrodes to a quadrant electrometer. 
It is claimed by both Kenrick and Frumkin that 
this method gives satisfactory results for aqueous 
solutions of inorganic electrolytes and for not too 
dilute solutions of derivatives of the hydrocar- 
bons. The errors inherent in the method of the 
“dropping electrode’”’ are discussed by Williams 
and Vigfusson.* 

Unfortunately the values of the surface poten- 
tials as determined for similar systems by differ- 
ent observers did not exhibit a good agreement. 
In an attempt to obtain better values Guyot and 
Frumkin independently’ developed a method in 
which the potential between a metal, coated with 
a radioactive substance, and a solution, is deter- 
mined. The method of Guyot and Frumkin was 
used by Frumkin and Williams*® and by Schulman 
and Rideal.’ 

The development of this method did not, how- 
ever, remove the difficulty that the values ob- 
tained by different workers, or even by the same 
worker, on seemingly identical systems, were 
extremely variable. The researches on surface 


1Sir Wm. Thomson, Lord Kelvin, Mathematical and 
Physical Papers, Vol. VI, Contact Electricity of Metals, 
pp. 110 to 145, Cambridge Press, 1911. See also Righi, J. de 
Physique [2] 7, 153 (1888). 

* Bichat and Blondlot, J. de Physique [2] 11, 548 (1883). 

* Kenrick, Zeits. f. physik. Chemie 19, 625 (1896); Guyot, 
Ann. de Physique (10), 2, 506 (1924); Frumkin, Zeits. 
f. physik. Chemie 111, 190 (1924); 116, 485 (1925); Garri- 
son, J. Phys. Chem. 29, 1517 (1925); Biihl, Ann. d. Physik 
[4] 84, 211 (1927); [4] 87, 877 (1928). 

* Williams and Vigfusson, J. Phys. Chem. 35, 348 (1931). 

* Guyot, Ann. de Physique [10] 2, 506 (1924). Frumkin, 
Zeits. f. physik. Chemie 116, 485 (1925). 

*Frumkin and Williams, Proc. Nat. Acad. Sci. 5, 400 
(1929), 

* Schulman and Rideal, Proc. Roy. Soc. A130, 259, 270, 
284 (1930). 


films which had been carried on for many years in 
this laboratory had shown that many films of in- 
soluble organic substances of the polar-nonpolar 
type are much less homogeneous than had usually 
been supposed in the literature of the subject. 
From this point of view it was suggested to 
Frumkin by one of us (Harkins) in 1928 that the 
discrepancies would disappear if simultaneous 
measurements of potential and film pressure were 
to be taken, the latter by the use of a film balance 
of the type used by Adam. This suggestion, re- 
peated by Frumkin to Rideal, was the origin of 
changes in methods abroad, and was taken as the 
basis of work in the present paper, which agrees 
remarkably well with that published recently by 
Adam and Harding.* The agreement is as good at 
low values for the area of the organic molecule as 
can be expected, especially when the great diffi- 
culty of obtaining really pure specimens of 
normal long chain acids is taken into account. At 
higher areas all of the difficulties disappear as 
will be shown later, if the potential is determined 
for all of the different regions of the film. 

Although the relation between the Kelvin- 
Righi method and that of Guyot and Frumkin 
seems to have been obscured in the literature of 
the subject? it is obvious that the latter method is 
the same as the former, with the substitution in 
the latter of the surface of a solution for that of 
one of the metals in the former. Thus the poten- 
tial determined by the method of Guyot and 
Frumkin is a contact potential. For obvious 
reasons in the determination of the effect of sur- 
face films it is more satisfactory to take the 
potential of the clean surface as an arbitrary 
zero, and to consider the change of potential 
caused by the presence of a film as the surface 
potential. 

A new method for the determination of the 
contact potential of metals and of the related 
surface potential of a solution has been devised by 
W. A. Zisman, " = who causes the metal elec- 


§ Adam and Harding, Proc. Roy. Soc. A138, 411 (1932). 
See also Schulman and Hughes, ibid. 430 (1932), and Fos- 
binder and Lessig, J. Frank. Inst. 215, 425 (1933). 

*See, however, a paper which has appeared since this 
paper was first written; Whalley and Rideal, Proc. Roy. Soc. 
A140, 484, 489, 497 (1933). 

10 Zisman, Rev. Sci. Inst. 3, 7 (1932). 

11 Zisman and Yamins, Physics 4, 7 (1933). 

22 'Yamins and Zisman, J. Chem. Phys. 1, 656 (1933). 
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trode to vibrate 250 to 500 times per second in- 
stead of coating it with polonium. This method 
of the vibrating condenser seems to give the 
same results as that of Guyot and Frumkin, 
but thus far it has been tested only with respect 
to the general order of magnitude of the values for 
the condensed film. 


II. APPARATUS AND PREPARATION OF FILMS 


The film balance, of the type designed by 
Adam" and modified by Freud and Harkins" 





Fic. 1. Apparatus for the simultaneous measurement of 
film pressure and surface potential, and for the determina- 
tion of the location and contour of continents and islands 
of condensed film. This view is taken inside the large 
thermally insulated metal box which surrounds the appa- 
ratus. The Compton electrometer is inside the smaller 
lead box in the left-rear. 



































Fic. 2. Diagram of apparatus shown in Fig. 1. 


13 Adam, Physics and Chemistry of Surfaces, Clarendon 
Press, 1930. 

4 Freud, Doctoral Dissertation, University of Chicago, 
1927. 
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(Figs. 1 and 2), is supported on four blocks of 
paraffin, P, inside a large grounded metal box 
supported on a heavy concrete pier. Heat insula- 
tion of the box is provided by a thick layer of 
wool felt on the outside. Inside this box is a 
smaller box L, the inner walls of which are heavy 
sheets of lead. This box holds the sens‘tive Comp- 
ton electrometer used to determine the potential. 
The electrical circuit is shown in Fig. 3. This is so 
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Fic. 3. Diagram of electrical connections for determina- 
tion of film potential. 7, trough of paraffin or brass; P, 
potentiometer; E4, one to six calomel cells; E, Compton 
electrometer; £;, battery for potentiometer; E2, standard 
cell, E3, battery for charging electrometer; S2, double- 
throw double-pole switch: upper position for setting R, 
with P at zero; lower position for measuring E, on poten- 
tiometer with zero reading for E. 


arranged that the potential may be either given 
by a potentiometer, or read directly from the 
electrometer deflection, which amounts to 680 
mm per volt at a scale distance of about 1.5 
meters. A much greater sensitivity is easily ob- 
tained. 

The front of the large box is occupied by two 
glass doors with metal frames, but the inside of 
the glass is covered with earthed copper gauze. 

As has been customary in this laboratory for 
many years, the entire operation of the film 
balance system, except for filling the trough and 
spreading the film, is carried out by means of 
handles outside the box. This procedure was also 
used with the electrical system. Thus the handle 
on the ebonite rod R controls the film balance, 
while that connected with the screw H controls 
the barrier B by means of two forked supports, 
one at each end of B. A handle on the rod R’ 
serves to shift the position of the silver electrode, 
coated with polonium, from side to side with re- 
spect to the surface of the solution, while the car- 
riage D is moved lengthwise of the film by another 
rod and handle not shown. Lengths along the 
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trough are measured by the scale S and those 
across it by S’. The electrode E is insulated by 
the amber supports A and its height is adjusted 
by the rack and pinion F’. 

The float F of thin sheet metal is attached to 
the end supports of the balance by vertical sheets 
of thin gold shaped as the letter S. The graduated 
circle G is read by a vernier V below it. The 
trough J is 25 cm wide and 60 cm long. The sur- 
face back of the float F is kept clean by sweeping 
with various barriers which are not shown. 
Troughs of glass, paraffin, and of brass carefully 
coated with paraffin, were found to give results 
which agree remarkably well. 


Preparation of films 


The organic acids used for the preparation of 
the films were the purest that could be obtained. 
Each sample was washed with a large quantity 
of pure hot water, and was then recrystallized 
twice from pure alcohol, and finally dried care- 
fully. 

A dilute solution of the film forming substance 
in a low boiling petroleum ether was prepared. A 
portion of this solution was put in a special weight 
pipette designed by Harkins and Freud in 1926. 
The carefully ground caps and stopper of these 
pipettes must be kept dry, and if this is done they 
prevent the evaporation of the solution during the 
process of weighing. Thus the error of 10 percent 
in the work of Schulman and Rideal caused by the 
determination of the amount of solution by 
counting the drops was avoided. Adam and Hard- 
ing and Schulman and Hughes have greatly re- 
duced the error below that of Schulman and 
Rideal by the use of the volumetric pipettes sug- 
gested by Fosbinder. 

It is obvious that with dilute films the distribu- 
tion of the islands is highly dependent upon the 
way in which the drops of the dilute solution are 
placed upon the surface of the liquid in the 
trough. An extensive study was made of the ef- 
fects on the surface potentials of dilute and con- 
densed films as affected by various methods of 
putting the dilute solution upon the surface. 


Electrodes 


Many experiments were carried out to deter- 
mine the best electrode to use with a radioactive 
substance. A platinum gauze electrode in the form 





of a small basket with from six to twenty radon 
tubes laid directly on the metal was found to 
give potentials which varied (1) with the number 
of radon tubes, (2) with the distance of the eiec- 
trode above the surface of the solution, and (3) 
with the previous treatment of the platinum. For 
example, the potential in millivolts against the 
surface of 0.1 molal hydrochloric acid was 95 after 
immersion of the platinum in cold cleaning solu- 
tion, — 200 after cathodic electrolysis, 312 after 
anodic electrolysis, and 520 after heating to red- 
ness in a Bunsen flame. After standing for several 
hours the potentials of the electrodes listed above 
became more nearly the same, but constant re- 
sults could not be obtained. 

Electrodes of various metals covered with 
polonium were much more satisfactory. Wires 
and disks of platinum, cadmium, zinc, tin, gold, 
copper and silver were tried. Good values were 
obtained with copper but even more constant 
values with silver, so this metal was used in all of 
the experiments reported. 

Connection with the solution was made by 
from one to six calomel half cells, which seemed 
to give better results than a large silver-silver 
chloride electrode. 


III. PoTENTIAL AND ForRCE-AREA VALUES FOR 
FitMs OF Myristic AcID 


Interesting results were obtained with mono- 
molecular films of myristic acid on 0.01 N hydro- 
chloric acid (Fig. 4), since the curves which give 
the variation of potential with the mean molecu- 
lar area of the myristic acid are much more 
simple than those of Schulman and Rideal. 

At molecular areas from 15 to 52 sq. A the 
values of the surface potential are represented by 
a moderately smooth curve. If any inflection at 
all appears at 32 sq. A, which is the location of a 
marked inflection in the pressure-area (FA) 
curve, it is a very minor one. At areas above 52 
sq. A the pressure falls to that of the gaseous film. 
The four sets of measurements at areas higher 
than this are of considerable value in connection 
with an interpretation of the characteristics of 
the film. In experiment 1, represented by circles 
with the lower half inked in, the film was at the 
beginning very dilute. As the film was compressed 
from high areas down to 57 sq. A the potential 
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Fic. 4. Lower left=pressure-area relations. Upper right=potential-molecular area 
relations. Films of myristic acid at 17°. The surface potentials are variable at molecular 
areas greater than 52A. The lower curve corresponds to the potential of a gaseous film, and 


the higher curve to that for an island. 


values increased along a moderately smooth 
curve. At 115 sq. A the surface potential was 
19 mv, and at 65 sq. A about 60 mv. 

In experiment 2 (open circles), the potentials 
for areas below 65 sq. A agree moderately well 
with those for experiment 1, but at higher areas 
they vary between 30 and 90 mv. In experiment 
3 (black circles) the potentials are much higher, 
and between 150 and 170 mv. In experiment 4 
(left half of circle open) the range of values was 
more limited, but of the same order as for 
experiment 4. 

In experiments 3 and 4 the silver-polonium 
electrode was directly over a continent or island 
of film in what has been called the ‘‘expanded”’ 
state, which represents a moderately closely 
packed film. In experiment 1 at 115 sq. A molec- 
ular area the electrode was largely over an area 
of gaseous film, but the potential was raised by 
adjacent continents or islands. As the area was 
decreased the area of gaseous film under the 
electrode was at first decreased only very slowly, 
but after the molecular area had decreased to 
less than 75 sq. A the continents of concentrated 
film occupied such a great fraction of the area 
that further decrease of area caused the gaseous 
film region to disappear rapidly, with a resultant 


rapid increase of potential. However, when the 
whole film became of the ‘“‘expanded”’ rather than 
of the gaseous type, the rate of increase became 
much more slow, since the gaseous phase was no 
longer present, and thus could not disappear. 

At areas below 52 sq. A the pressure was above 
that of the gaseous film, and the continents of 
concentrated film were pressed together, so be- 
low this area a single curve is sufficient to repre- 
sent the potential values. However, all irregulari- 
ties do not always disappear, since the effect of 
the contour of the continents is not entirely ef- 
faced. Between 53 and 22.5 sq. A for the molec- 
ular area the results given in the figure are those 
of experiments 2 and 3, and are represented by a 
single smooth curve. Curves of this type were 
obtained in a considerable number of similar ex- 
periments between 10° and 30°. The potential 
exhibited a general decrease as the temperature 
was increased. In a few experiments a minor in- 
flection occurred at the same area as that at 
which the marked inflection occurs in the force- 
area curve. 

A geographic plot of the variation of the sur- 
face potential for films of myristic acid at 23° is 
given in Fig. 5. The variation in potential at the 
different areas was as follows: 
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Molecular area in sq. A 56.5 37.1 25 
Minimum potential in mv 0.0 155 238 
Maximum potential in mv 176 175 281 


In general the variation of the potential de- 
creased as the area decreased toward 20 sq. A. At 
12° Adam and Harding find a variation of 166 mv 
at a molecular area of 70.6 sq. A, with a minimum 
of 1 mv and a maximum of 167 mv. 


Palmitic Acid 





























Area 
48 A2 68 217 315 330 
Barrier—> 35 59 141 200 |<—Float 

233 130 50 61 

282.2 A? 345 315 324 315 
297 310 311 312 
319 314 308 305 
331 313 313 331 

Pentadecylic Acid 

48.5 A? 155 163 170 175 
160 20 165 170 
165 24 101 170 

39.7 A2 179 183 183 180 
179 180 181 180 
175 179 180 179 
175 179 181 186 











Myristic Acid 





56.5 A? 10 22 31 42 90 165 158 
0 15 30 38 164 160 156 
5 35 51 165 167 165 165 

48 98 166 174 171 176 164 

145 148 170 169 170 103 199 











37.1 A2 175 170 170 160 157 
172 165 168 160 155 
169 166 167 164 158 
166 169 188 165 155 
163 167 190 195 — 











25.0 A? 255 248 255 281 
255 248 255 245 
255 252 258 243 
254 277 261 238 
253 254 261 256 











Fic. 5. Distribution of surface potential values according 
to the location of the electrode over the film. The higher 
values correspond to the more concentrated portions of 
the film. Barrier on the left, and float attached to film 
balance on the right. Palmitic acid film at 23°C. Other 
films at 17°C, 


IV. IsLANDS AND OTHER NON-HOMOGENEOUS 
AREAS IN FILMS 

One of the first steps in the work reported in 
the present paper was to study the effect on 
the surface potentials at different positions of 
various methods of dropping the dilute solution 
of the film forming substance on the surface. 
For this purpose the silver-polonium electrode 
was mounted in such a way as to be easily moved 
in either direction from barrier to float. It was 
found that when the potential was plotted against 
the distance a high peak or a wide plateau ap- 
peared near the position where a drop had been 
placed, provided the measurements were made 
soon after the drops had been deposited. Often, 
however the islands move quickly from their 
original positions. The best results are obtained 
with large islands anchored by some solid object. 

The preceding section illustrates the existence 
of such islands in gaseous films of myristic acid. 
These appear as floating independent islands at 
molecular areas higher than those given in Fig. 5, 
where the islands have largely merged into con- 
tinents. At extremely great molecular areas the 
islands disappear or become less prominent. 

With a film of palmitic acid at 23° and a 
molecular area of 48 sq. A, the surface potential 
varied from 35 to 330 mv, while at 28.2 sq. A the 
potential had become more constant and between 
297 and 345 mv. That the potential is not more 
variable than this is surprising, as the film pres- 
sure is still very low. At lower areas and higher 
pressures the potential becomes still more con- 
stant. 

With pentadecylic acid the film pressure rises 
at much higher molecular areas, so that, while at 
48.5 sq. A the potential varies from 20 to 175 mv, 
at the somewhat high area of 39.7 sq. A most of 
the variation has disappeared and the potential 
varies only from 175 to 180 mv, except in one 
corner of the film where one value was obtained 
at 185 mv. 

At lower areas than those given in Fig. 5 the 
potentials in general become more constant and 
except for some peculiar positions, such as a 
corner of the film, do not vary more than 3 to 5 
mv. However, at molecular areas much below 
20 sq. A the potential again becomes more 
variable, and in general decreases, since most 
films crumple at these lower areas, and regions 
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are produced in which the molecules of the film 
do not preserve their regular orientation with 
respect to the surface of the liquid. 


V. GENERAL POTENTIAL RELATIONS FOR FILMS 
oF Fatty AcIps 


The surface potential of a fatty acid film in- 
creases rapidly with the length of the hydro- 
carbon chain. Thus at 20 sq. A the potentials for 
stearic (18), palmitic (16), pentadecylic (15), 
myristic (14), and lauric (12) acids, are respec- 
tively 400, 383, 375, 336 and 276 mv. 

Each curve in Fig. 6 represents a single experi- 
ment. At areas greater than that at which, ac- 
cording to the force-area curves, the pressure ap- 
parently falls to zero, other experiments gave 
other results, but at smaller areas there was good 
agreement. 

The presence of the double bond of oleic acid 
produces a lowering of about 200 mv, a sur- 
prisingly great lowering of the surface potential. 
The presence of the single double bond lowers 
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Fic. 6. Surface potentials (P) and film pressure values 
(F) for films of organic acids on 0.01 molal hydrochloric 
acid. The points plotted correspond to individual experi- 
ments. The potentials for areas greater than that at which 
the special force-area curve apparently merges with the 
x-axis, are uncertain. Temperature, 17°. 
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the potential to about the same extent as the 
removal of six carbon atoms. 

The values for myristic acid (14 C) given in 
Fig. 4 lie below those for pentadecylic acid (Fig. 
6). The flat portion of the curve lies at about 150 
mv for the former, and 170 mv for the latter. 


VI. ALCOHOL FILMs 


The two higher alcohols investigated (Fig. 7) 
show the same relation as the acids in that the 
potential is higher for the alcohol with the longer 
chain. The surface potential is, however, higher 
for the alcohols than for the corresponding acids. 


VII. PARAFFIN LAYERS 


Nonpolar hydrocarbons on a clean water sur- 
face exhibit different potential relations from a 
fatty acid. Thus when paraffin was spread from a 
solution in ligroin, the potential of the interface 
became 60 mv more positive with respect to the 
clean surface when the area was decreased to one- 
fourth the initial area. At the same time the pres- 
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Fic. 7. Surface potentials (above) and surface pressures 
(below) for films of normal alcohols. As with acids the 
potential increases with the length of the carbon chain. 
Temperature, 17°. 
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sure exerted on the float was of the order of one 
dyne per cm. The film of paraffin was visible to 
the eye. 


VIII. SurFACE POTENTIALS, DIPOLE MOMENTs, 
AND THE ORIENTATION OF MOLE- 
CULES IN SURFACES 


The surface potential of a film, defined as the 
contact potential of the surface covered by the 
film minus that of the clean surface of the under- 
lying phase, is undoubtedly related to the dipole 
moments of the film producing substance, to the 
extent of the orientation of the molecules, and to 
other factors, some of which are geometric. 

The orientation of the molecules in a surface 
was supposed by Hardy" to be due to the un- 
symmetrical field of the molecules, and Harkins, 
Davies and Clark'* took a somewhat similar view 
when they assumed that the potential differences 
at surfaces, as well as the other energy relations 
may be explained on the basis of the polar setting 
of molecules in surfaces. An addition to this elec- 
trical theory of surfaces was made by Harkins 
and King!’ who consider especially the relations 
between surface energy, cohesion, and the 
polarity of the molecules which constitute the 
legends or liquids involved. 

In a monomolecular film on water, on any 
other liquid, or on a solid, the molecules in the 
film, and to some extent the adjacent molecules 
of the surface region, orient themselves in such a 
way as to decrease the free surface energy. 

In this discussion it will be assumed that just 
as a stone in a position above the surface of the 
earth possesses an energy of position, usually 
designated as its potential energy, a molecule in 
the surface possesses an energy of position which 
will also be termed its potential energy. The po- 
tential energy of the weight may be referred to an 
arbitrary zero position at the level of the earth’s 
surface, while that of the molecule is referred to a 
position in the liquid phase underneath the sur- 
face region. The work done in lifting the weight 
is done against the force of gravitation, while that 
of lifting the molecule is opposed by the cohe- 


* Hardy, Proc. Roy. Soc. A86, 634 (1911-12). 

* Harkins, Davies and Clark, J. Am. Chem. Soc. 39, 541 
(1917), 

" Harkins and King, J. Am. Chem. Soc. 41, 970 (1919), 


sional forces which may also be considered as the 
resultant of the electrical forces. 

The potential energy of the stone with refer- 
ence to its zero position, may be determined by 
lifting the stone, initially at rest, to the specified 
position. The potential energy is said to be equal 
to the work done. If the stone, when at the level 
of the earth’s surface, had a velocity in the up- 
ward direction, then the work done would be 
less, but the final energy of position at a specified 
height would not be changed. 

A molecule inside the liquid can in general 
move into the surface region only when its energy 
of position, which we will here call its potential 
energy, is supplied in the form of work, heat, or of 
some other class of energy. 

The total energy necessary to form one sq. cm 
of surface is supplied partly as work vy and partly 
from the kinetic energy of the molecules them- 
selves. The molecular kinetic energy is the source 
of almost the whole latent heat / of the surface. 
Thus a surface has a lower temperature than the 
liquid when it is first formed. 

The value of the heat function / for unit area 
of a surface is given by the relation 


h=y+l. (1) 


In an isothermal change h is the total amount of 
energy which must be added to the system when 
a unit area of surface is created, and is usually 
designated as the total surface energy. 

It will be assumed, on the basis of our general 
knowledge of surface energy relations that 7, the 
sum of the potential energies of all of the mole- 
cules in 1 sq. cm of the surface region is very 
nearly equal to / the energy supplied to form 
unit area of the surface, or 


a=h. (2) 


The orientation of the molecules causes a de- 
crease in the free energy (vy) of the surface, while 
on the whole the latent heat is not greatly 
changed, so the total surface energy h is decreased 
by the orientation. Thus the sum (7) of the po- 
tential energies of the surface molecules is also 
decreased (Eq. (2)) by the orientation below 
that for a random distribution. Since the poten- 
tial energy of the molecules in the surface is as- 
sumed to be due to electrical forces, it may be 
called the total electrical energy. 
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The mutual energy due to the interaction of 
the dipole moments of the molecules constitutes 
only a part of this electrical energy. In the older 
theory the remainder of this energy is attributed 
to quadrapoles or multipoles but more recently 
the theory has been modified by the introduction 
of some of the principles of quantum mechanics, 
as, for example in papers by Lennard-Jones.'® 

While thermodynamics predicts only the de- 
crease of free surface energy as associated with 
the molecular orientation, the above discussion 
seems to indicate that this is in general accom- 
panied by a decrease in the total electrical sur- 
face energy. While it is not from these relations 
certain that the electrical energy of the dipoles 
must decrease with the decrease of total elec- 
trical energy, the dipoles make an important 
contribution to this energy, so it seems probable 
that such a decrease occurs. Also from the known 
behavior of large scale dipoles which are free to 
change their orientation and position it may be 
deduced that the molecular dipoles will shift and 
orient themselves in such a way as to decrease 
their mutual electrostatic energy. However, the 
decrease may not be to a minimum since it is 
possible that there may be in some cases a re- 
sultant increase of the energy of the multipoles. 
The dipoles involved include those of the liquid 
(usually water) as well as those of the organic 
molecules of the film. 

The potential difference (AV) between two 
parallel plane metal plates of an ordinary parallel 
plate condenser, with a unit dielectric constant, is 


AV =A4nes, (3) 


in which oa is the surface density of the charge and 
s is the distance between the plates. 

Since o=ne, where u is the number of elec- 
tronic charges per sq. cm and e¢ is the electronic 
charge, 

V =4rn(es) =4rnm, (4) 


in which m is the product of the value of the 
electronic charge by the distance between the 
plates. 

Schulman and Rideal!® consider that this equa- 


18 Lennard-Jones, Proc. Phys. Soc. (London) 43, 461 
(1931). 

19Schulman and Rideal, Proc. Roy. Soc. 
(1931). 
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tion, which they attribute to Helmholtz, may be 
applied to calculate the dipole moment of the 
molecules of the film on the basis that u,=m if n 
is considered as the number of molecules in one 
sq. cm of film. Thus 


us =AV/4rn. (5) 


They find, however, that the value thus obtained 
is only about one-sixth, or even less, of that de- 
termined by the standard methods. Their ex- 
planation of the discrepancy may be found in 
their papers. 

Our point of view is different in that we con- 
sider that the equation of the parallel plate con- 
denser is not valid for the case of a film. 

(1) The equation assumes that the charges are 
on metal plates, which is far from true, since they 
are assumed to be associated with the molecules. 

The correct equation for parallel plates of an 
insulating material is 


AV=2rnm’. (6) 


If now it is assumed that u,=m’, then the value 
of u, thus obtained is, obviously, much closer to 
the standard value of the dipole moment. How- 
ever, this does not show that Eq. (6) is suited to 
the calculation of the dipole moment. 

(2) The Eqs. (5) or (6) were developed for a 
large scale condenser. There is no evidence that 
they are valid for a condenser in which the 
charges are attached to molecules, and are a 
distance apart which is of the atomic scale. 

(4) The charges of this atomic scale condenser 
are only imaginary or apparent charges, though 
the dipole moment is real. 

(5) Even if these charges are assumed to be real 
and to have a definite position, the charges would 
move with the heat vibrations of the molecules in 
such a way that any instant both positive and 
negative charges would probably be on both 
sides of a plane which represents the height of the 
mean position of all of the imaginary charges, 
both positive and negative. Even if this is not the 
case the charges do not lie in a plane. 

(6) Both the molecules and the positions of the 
imaginary charges are in rapid motion, which is 
not taken into account in the development of the 
equation. 

Thus yu, does not seem to be a dipole moment, 
but is a quantity which may be found useful, 
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since it is equal to the surface potential divided 
by the number of molecules in unit area of the 
film and this quotient multiplied by the con- 
stant 1/(47). It seems preferable to make the 
constant unity. This gives 


m=AV/n. (7) 


Either m or yu, is thus related in a simple way to 
the surface potential and the concentration of 
the molecules in the film, and in a complicated 
and unknown way to the standard dipole mo- 
ment. The values of m are given in Table I. 


TABLE I. Values of m=V/n for films of fatty acids and 
alcohols, as calculated from surface potentials. 


Substrate=0.01 N HCl (for the lauric acid films this 
was nearly saturated with NaCl to reduce the solubility of 
the lauric acid). 








m X10" volts molecule“! cm=~? 





at area a 

Compound Do & DD © 8 SA 
Lauric acid iS 29 20 22 
Myristic acid ni 20 23 28 th 
Pentadecylic acid 24 26 2.5% 234 
Palmitic acid 2.54 2.94 
Stearic acid 24. 3a 
Oleic acid 16, 20 2.i¢ 2.2 2225 
Myristyl alcohol 2.66 2.7 
Cetyl alcohol 2.9 








In the determination of the factors which are 
of the most importance in their effect on the 
surface potential, the variation of this potential 
with temperature is of prime importance. In the 
experimental work the surface potential was 
found to usually exhibit a decrease with rising 
temperature, but a rise in temperature often 
changed the state of the film. The effect of in- 
crease of temperature on the surface potential of 
condensed films will be studied more carefully 
later. 


IX. GENERAL DISCUSSION 


There are certain features of the potential-area 
curves which are of importance. There is no indi- 
cation in these curves of the numerous ‘‘surface 
phases,”’ such as the “‘homalic”’ and ‘“‘prehomalic”’ 
phases of Schulman and Rideal. Since these 
“phases” are not indicated by any known force- 
area determinations it is evident that there is no 
evidence that such phases have any existence, 
especially since both Adam and Schulman and 


Hughes have also failed to find them. Doubt is 
also thrown upon the value of the consideration 
of condensed and expanded films, as defined by 
Adam, as distinct phases. If films in general, after 
being compressed to a molecular area of about 
20 sq. A, are subjected to a constantly reducing 
pressure, it is found that there is considerable 
hysteresis, both with respect to pressure and po- 
tential. With increasing pressure there is in gen- 
eral no special variation in the curvature of the 
potential-area curve at the point where the force- 
area curve indicates what has been considered to 
be the transition point between expanded and 
condensed phases. This fact seems to indicate 
that there is no sharp change in orientation or 
structure in the film at this point. The above con- 
siderations do not show that the terms ‘“‘con- 
densed”’ and “‘expanded”’ with reference to films 
should be abandoned, but they indicate that 
there is a need for a more detailed study of the 
reversibility of the supposed transition between 
the two states. 

As has been stated earlier the curves for the 
potential-area variations as given in this paper, 
represent individual sets of determinations. 
These have been checked by numerous sets which 
have not been reported. The curves presented are 
not the smoothest which were obtained, but 
they were taken from several sets in good agree- 
ment with each other. The curves exhibited by 
Adam and Harding, and by Schulman and 
Hughes, in recent papers are in general either 
averaged or idealized, so that no specific data are 
presented. In some cases the same substances 
were used by Adam and Harding and in the 
present investigation, with no apparent disagree- 
ment in the results. However, since our tempera- 
tures were higher the potentials obtained are 
lower. 

The specific experiments on films of organic 
acids, as presented in Fig. 6, exhibit several 
features of general interest. For example the 
surface potential for films of stearic acid at areas 
greater than 27.5 sq. A may be as high as those 
given by the curve, or lower values may be found. 
However, in this individual set of determinations 
the potential-area curve gives no indication of 
any transition whatever. The same is true of 
pentadecylic films above 40 sq. A. Of the two 
condensed films that of stearic acid is much more 
























































coherent than that of pentadecylic acid. With 
such highly coherent films as that of stearic acid 
the surface potential is lowered only slightly as 
the pressure on the film is reduced. Thus it ap- 
pears that the molecular orientation in such a 
portion of a solid film is not materially affected 
by the release of pressure on the film to almost a 
zero value. If time is given for a two dimensional 
evaporation of the film substance, at extremely 
low pressure, the surface potential is reduced to 
a few millivolts. 
The surface potential is found to vary with 


1. the length and shape of the hydrocarbon part of the 
organic molecule, 

. the nature and number and position of the polar 
groups in the molecule, 

. the hydrogen ion concentration of the solution, 

. the temperature, 

. the concentration of the film forming substance in 
the film, 

. the nature of any ions of salts present in the solution, 

. time, including phenomena of hysteresis. 


bt 
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With changes of acid concentration between 
0.01 and 0.1 N the surface potentials for films of 
acids and alcohols was found to be only slightly 
affected. 

The effect of hydrogen ion concentration may 
be illustrated by experiments on polymers of 
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hydroxydecanoic acid with molecular weights 
from 208 to 16,900 and general formula 
HO[(CHe)sCOO ],H. A change from 0.01 normal 
hydrochloric acid to the same concentration of 
sodium hydroxide increased the molecular area in 
the condensed film by about 30 percent and de- 
creased the surface potential by ten to twenty 
percent. Proteins were even more largely affected 
by the hydrogen ion concentration, as corre- 
sponds with the work of Gorter. The details of the 
effects of hydrogen ion concentration upon the 
areas and surface potentials of organic substances 
of high molecular weight will be left to another 
paper. 

In order to test the assumption presented in 
Section I of this paper, that the method used in 
this work is suitable for the determination of 
either the surface potentials of solutions or the 
contact potentials of metals, the carefully cleaned 
surfaces of several metals were substituted for the 
liquid surface. The differences in these values 
were found to give the contact potentials of the 
metals as accurately as they are known for like 
conditions. It is obvious that the calomel cell is 
eliminated from the circuit for such a measure- 
ment. 

The writers wish to thank Dr. Julius Stieglitz, 
and the Chemical Foundation, whose aid has 
made this work possible. 
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Acoustical Studies. II.! The Behavior of a Gas with Several Independent Internal 
Energy States 


WitiiaM T. Ricuarps,? Princeton University 
(Received August 24, 1933) 


The method founded by Einstein for the description of 
the acoustical behavior of a dissociating gas has been 
extended to cover a nondissociating gas in which five 
groups of internal energy states have different relaxation 
times. Since the resulting expressions are unwieldly, 
approximations based on them, which permit the rough 
description of experiment, have been given. The study of 
the variation of the velocity of sound with frequency 
demonstrates differences in the relaxation times of the 
various states. The variation of the velocity of sound 


with pressure shows in addition effects due to three-body 
collisions or to radiation from optically active states. It 
has been suggested to treat the temperature coefficient of 
the kinetic relaxation time by the introduction of an 
empirical quantity which may be called ‘‘the activation 
energy of collision.’”” By this means a rough kinetic analysis 
of the transition probability is possible. In gaseous mix- 
tures, the relative spatial and energetic effectiveness of 
various types of collisions in exciting internal energy may 
be compared. 





INTRODUCTION 


REAT progress has been made during the 
past few years in the measurement of the 
velocity of sound in gases. Because of the 
invention of the sonic interferometer by Pierce® 
it is now possible to determine this quantity in 
a gas of any desired purity with an accuracy of 
0.1 percent from audible frequencies to fre- 
quencies somewhat above a million cycles per 
second. There is every reason to suppose that, 
with sufficient incentive, even this accuracy may 
be considerably improved. Similarly, interfero- 
metric absorption measurements, particularly 
owing to the investigations of Hubbard,‘ are 
rapidly becoming capable of exact interpretation. 
It is worth while, therefore, to discuss in detail 
the information which may be expected to be- 
come available from extremely accurate acousti- 
cal measurements. 
Two methods have been used to describe the 
propagation of sound in acoustically dispersive 
gases. The first of these evaluates the absorption 





‘It is considered that the communications of Richards 
and Reid, J. Chem. Phys. 1, 114,737 (1933), are Parts I 
and Ia of this series. 

*Fellow of the John Simon Guggenheim Memorial 
Foundation. 

* Pierce, Proc. Am. Acad. Sci. 63, 1 (1928). 

* Hubbard, Phys. Rev. 41, 833 (1932), and earlier papers 
there cited, 


863 


and dispersion by calculating as a function of 
frequency the number of molecules making the 
transition between two given energy states under 
the influence of the adiabatic acoustical cycle. 
It was founded by Einstein,’ who considered 
the case of a dissociating gas, and has recently 
been extended by Kneser® and by Rutgers’ to in- 
clude nondissociating gases. The second method 
uses the conception of the fictional temperature 
of an energy state or group of energy states, 
and calculates as a function of frequency the 
relation between this fictional temperature and 
the temperature of the translational terms in 
the acoustical cycle. It was first published by 
Herzfeld and Rice® and was generalized immedi- 
ately thereafter by Bourgin® to cover gaseous 
mixtures and effects due to radiation. The two 
methods are equivalent descriptions of the same 
process, and with similar postulates must yield 
similar conclusions. They have converged upon 


5 Einstein, Sitz. Ber. Akad. 380 (1920). 

6 Kneser, Ann. d. Physik 11, 761 (1931); 16, 337 (1933). 

7 Rutgers, Ann. d. Physik 16, 350 (1933). 

8 Herzfeld and Rice, Phys. Rev. 31, 691 (1928). Among 
the numerous and important conclusions reached in this 
paper is included a demonstration that absorption due to 
viscosity and heat conduction plays a minor réle if a heat 
capacity effect is also present. This is formal justification 
for disregarding this quantity in using Einstein’s method. 

® Bourgin, Nature 122, 133 (1928); Phil. Mag. 7, 821 
(1929); Phys. Rev. 34, 521 (1929) and 42, 721 (1932). 
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expressions for the velocity and absorption of 
sound which are practically identical, and are 
suitable for the description of experiment if the 
dispersion is produced by the gradual disappear- 
ance from the heat capacity at high frequencies 
of the transition from the zero point energy state 
of a vibration to its first excited state. 

The purpose of the following discussion is to 
amplify these expressions to include as many 
transitions as may practically be necessary. The 
criterion of practical necessity is the calculation 
of the velocity of sound at very low frequency 
within limit of experimental error from the same 
relationship which is used to describe the disper- 
sion. The method of Einstein has been used in 
preference to that of Herzfeld-Rice-Bourgin 
because of its greater simplicity. It must again 
be emphasized that either method is suitable for 
the purpose, and that the description of cases of 
any desired complexity is merely a matter of 
industry. Sections 1, 2, 3 and 4 of the following 
study are devoted to the development of ex- 
pressions for the velocity and absorption of sound 
suitable for use with gases of increasing numbers 
of frequency dependent internal energy terms. 
Section 5 considers the bearing of the earlier 
sections on experiment, and shows how suitable 
measurements may lead to an increased knowl- 
edge concerning collisions of the first and second 
kinds. An effort has been made to make Section 
5 intelligible to a reader who does not wish to 
follow the mathematical detail of the earlier 
Sections. Finally, in Section 6, gaseous mixtures 
and more complex transitions in pure gases are 
briefly considered. 


1. THE METHOD OF EINSTEIN 


If we consider the propagation of sound in a 
pure nondissociating gas which has a molecular 
heat capacity at constant volume greater than 
(3,2)R, we may expect with increasing frequency 
one of two alternative types of behavior. Either 
the rotational and vibrational energy will be 
brought into equilibrium by each collision, in 
which case the gas will not become acoustically 
opaque until frequencies of the order of magni- 
tude of the collision frequency have been reached, 
or the adjustment of rotational and vibrational 
energy may be hindered by some spatial or 
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energetic property of the colliding molecules, in 
which case it will lag behind the translational 
adjustments and the gas will become partially 
opaque at lower frequencies. It is with the second 
of these possibilities which we are here concerned. 
If the internal heat capacity is due to changes 
between many states, each of these changes 
may be characterized by a separate rate of 
transition. The method of Einstein may be 
adapted to supply a quantitative description of 
the rates of these transition processes. 

If w is the excess pressure in an element of 
volume of the gas at a distance x from the source, 
and A is similarly the density increment of this 
element, one may write for the time variations 
of these quantities 


m= moeiteivtt—az) (1.1) 
and 


A =Apeie(t-22), (1.2) 


where w is the frequency of sound in radians, 
and ¢ refers to time. The velocity and absorption 
of sound are related to a by the expression 


a=(1/V—18/w) 


if V is the phase velocity (27 V=d) and 6 the 
absorption coefficient for the amplitude ex- 
pressed as a reciprocal distance. The propagation 
of sound in an absorbing medium is therefore in 
close analogy to the propagation of light in a 
medium with a complex refractive index. The 
angle ¢ expresses the phase lag between pressure 
and density increments” and characterizes both 
absorption and dispersion. From the real and 
imaginary parts of (1.1) and (1.2) and the equa- 
tion of motion of a plane compressional wave it 
is found that 


cos g=(mo/Ao)(1/ V?— B?/w") 
sin ¢ = (10/Ao)(2B/wV). (1.5 


Einstein has demonstrated that the quantity 
x/A corresponds to the square of a complex 
velocity of sound, and from this it follows that 


(1.3) 


(1.4 
and 


r.p. t/A=(0P/dp)adiav. = V2u (1.6 


1° Kneser, second reference under footnote 6, has given 
an illuminating specific picture of this lag and the irre- 
versible energy change it implies. 








IS 


1S 


1) 


nm 


lex 
iat 





and that 
ip. A/r = —128/wV = —th/2xrV?, (1.7) 


when hk, in (1.7), designates the absorption 
coefficient for the intensity expressed as a 
reciprocal number of wave-lengths. The factor 
win (1.6) is 


u=(1—h2/1672)/(1+h2/1672)2 (1.8) 


and may be taken as unity except in the neigh- 
borhood of the absorption maximum. 

In general there will be absorption only when 
the dispersion d V /dw is appreciable.® It is difficult, 
however, to formulate the relation between the 
two, since it depends upon the characteristics of 
7/A. It is more convenient, therefore, to compare 
absorption and velocity, which are more simply 
related as in (1.4). This may be expressed by 
combining (1.4) with (1.6) to give 


(10/Ao) (1 +h?/162?) = V?. (1.9) 


The complex velocity of sound and the prop- 
erties of the gas are related by the characteristic 
acoustical equation. This may be written 


a/A=(P/p)(1—APv/PAv) (1.10) 


if P and p refer to pressure and density re- 
spectively, and APv and PAv are complex energy 
increments which vary with time about their 
respective amplitudes by a factor e*‘, and which 
are out of phase in the dispersive region of 
frequency. The problem of calculating the 
acoustical behavior of a gas is therefore the 
problem of evaluating APv and PAv in terms of 
pressure, frequency, temperature, and quantities 
descriptive of rate of energy adjustment and 
heat capacity. This is conveniently effected by 
means of the equation of state and the equation 
for adiabatic change in the gas, the equation for 
adiabatic change being related to frequency by 
means of a series of reaction equations. 

In finding APv the gas will be considered ideal, 
and effects due to departures from ideal behavior 
later discussed. For an ideal gas the displaced 
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equation of state is 


APv=RAT, (1.11) 


where R and v are considered gram molecular 
quantities, and AT refers to a temperature 
displacement of the translational terms. A rela- 
tion between AT and the displacement of the 
internal energy is given by the equation for 
adiabatic change, which may be written 


— PAv =(3/2)RAT+E, © kAm, 
k 
+E. ¥ lAnj+-++. (1.12) 
l 


The first term on the right refers to changes in 
translational energy, and the subsequent terms 
express energy changes in sets of rotational and 
vibrational states corresponding to the fre- 
quencies 4.=£,/hNa, ve=E2/hNa, etc., h being 
Planck’s constant and N4 Avogadro’s number. 
The symbols m;, m), etc., refer therefore to 
fractions of gram molecular weights in the kth 
and /th, etc., states, and are subject to the 
conditions }/,m,=1, etc., and }>,An,=0, etc., 
owing to the law of the conservation of matter. 
The actual number of molecules per cubic 
centimeter in the kth state is Nn,, if N is the 
total number of molecules per cubic centimeter. 
If C, is the equilibrium heat capacity due to the 
first set of states, and AT, its fictional tempera- 
ture displacement in the sound wave, £,>°.kAn, 
= C,AT; where AT=AT,=0. The calculation of 
the velocity and absorption of sound in a pure 
gas of given properties follows from an appro- 
priate development of (1.12) on the basis of a 
reaction equation, and the subsequent substitu- 
tion of (1.11) and (1.12) in the characteristic 
acoustical equation (1.10). 

In Sections 2, 3, and 4 of this communication 
the characteristic acoustical equation of the gas 
will have the form 
x P R 
—=—/14 , (1.13) 
A p C.+Rx?no[ (a+b) /(c+id) ] 


where C,, is that part of the heat capacity which is 
independent of frequency in the range considered, 
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x=E/RT=hv/kT, and the constants mm, a, b, c, and d vary with the complexity of the internal 
energy of the gas. According to (1.6), therefore, the velocity is given by 


r 





Vu=—{1+R 
p 


(C,,+Rx?noa/c) + eee | P 
a es 
(C..+ Rx?noa/c)?+(C,,d/c+Rx?nob/c)? 


(1.14) 
p 


where ., may be called the ‘‘apparent heat capacity ratio” at the frequency w, and the dimensionless 
absorption coefficient, from (1.7) to the first order of small quantities is 





(C+ Rx?noa /c)(C,d/c+Rd/c+Rx*nob/c) —(Cd/c+Rx?nob/c)(C,+R+Rx?noa/c) 
™~ (C,.+Rx?noa/c)(C,+R+Rx?noa/c) +(C,d/c+Rx?nob/c)(C,.d/c+Rd/c+Rx?nob/c) 


}. (1.15) 


Suitable constants for substitution in (1.14) and (1.15) will now be evaluated. 


2. A GAS WITH THREE SETS OF ENERGY STATES 


The simplest case of practical interest is that 
in which the dispersive region results from a lag 
in the adjustment of a single quantum transition, 
as, for example, the transition from the zero 
point state to the first excited state of.a defor- 
mation vibration. To describe this various parts 
of the energy may be considered as contained in 
three types of energy states. The states EF; are 
the translational states. The states £,; contain 
in addition all rotational and vibrational states 
in their equilibrium distribution with the excep- 
tion of a single mode of vibration, which is 
considered to have its zero point energy. Eyry 
differs from E,; in that the single mode of 
vibration is considered to be in its first excited 
state. It is then the gradual dropping out of the 
transition Ey@E;;; from the effective heat 
capacity at high frequencies which produces 
dispersion and absorption. For the equation for 
adiabatic change (1.12) may then be written 


— PAv = (3/2)RAT 


+ (0(Ey — Ey) /0T),AT+(Emr—En)An, (2.1) 
or, in more abbreviated form, 
— PAv=C,AT+EAn. (2.2) 


This case will be discussed with pedantic atten- 
tion to detail in order to illustrate the method 
of calculation. 

In going from (2.1) to (2.2) a step has been 
taken which is not permissible without further 
justification. By including in C,,, which repre- 
sents the effective heat capacity of the gas at 
frequencies so high that Am,=0, all the heat 








capacity except that due to the E,;,y,@Ey, tran- 
sition, two assumptions have been made. The 
first, which supposes that a terminal velocity 
corresponding to 


V..= {(P/p)(1+R/C,,)}3 


may actually be found can, at present, be 
justified only by experiment. The second, which 
supposes that, in case all energy states except 
one are adjusted at a rate much greater than 
the frequency, their influence on the acoustical 
behavior of the single slowly adjusted state is 
negligible, is justified by the more complicated 
models which are developed in Sections 3 and 4.** 

At mitall frequencies of sound the elasticity 
of the gas has equilibrium value, and (2.2) it 
becomes 


— PAv=C,AT+E(an,/aT),AT=CAT. (2.3) 


From Boltzmann’s principle m/m=e-* and m 
=(1+e-*)—'. Hence, in the three-state gas, 


Co— C,, = Rx*e—*/(1+e7*)?. (2.4) 


It is now necessary to relate Am and An, to 
the acoustical frequency. This relation is given 
by the reaction equation, which may be written 


—dno/dt =foino—fiom1 (2.5) 


if ¢ denotes time and m, m the fractions of 
molecules having the energy of the first excited 
state and the zero point state respectively. The 


% Bourgin (reference 9, 1932 paper) has demonstrated 
this explicity in his three-state model by carrying all 
rapidly adjusted states through the entire calculation 
with specific finite transition probabilities. 
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symbols fo: and fio are used to denote the proba- 
bilities of the transitions 0-1 and 1-0 due to 
any cause; they are taken, for reasons which 
will later become apparent, as rates of excitation 
and de-excitation per molecule per second. Thus 
fi is the mean number of times per second 
which a molecule in state 1 leaves state 1 for 
the state 0; it may be expected to depend on 
pressure and temperature. At equilibrium fiom 
= fom, and therefore fo: = fie”. 

Eq. (2.5) is now displaced from equilibrium 
by small, finite increments of the independent 
variables, these increments being written Amp, 
AT, and AN. WN as before is the numerical 
density and AT refers, as always, to the temper- 
ature of the translational terms. The three 
partial differential coefficients resulting from the 
displacement of the left-hand side will not be 
explicitly written. Since fio, the equilibrium value 
of m, has become 7ip+ Any it is sufficient to write 


d(No+Ano) /dt=dAno/dt=iwAno, (2.6a) 


the evaluation of the time derivative following 
from the fact that Amp varies about its amplitude 
according to a factor e“*'. The displacement of 
the right-hand side involves three sets of partial 
coefficients. Displacement with respect to mp 
gives, if fo: and fi are independent of mo, 


(—- 


Ono 





) Ano= — (fiotfor)Am, (2.6b) 


because in this case Amp+An,=0. The displace- 
ment with respect to T gives, in close analogy 
to the van’t Hoff isochore, 


(O(fo1%0—f 1071) /8T) non AT = (fiomx/T)AT. (2.6c) 


The result of the change N=N+AN may be 
represented in functional form, or some relation 
between fo, fio and pressure may be assumed. 
It is more convenient to take the latter course, 
and to assume that 


foi/N =kor and fio/N =k, (2.7) 


Where ko: and ky correspond to the bimolecular 
rate constants of the reaction No +N@N,+N. 
It is supposed, therefore, that the excitation 
of vibrational energy takes place solely because 
of a double collision of the first kind, and that 
its de-excitation is the result solely of a double 


collision of the second kind. It is further sup- 
posed that any molecule is equally capable of 
causing these transitions whether or not it is 
itself excited. This assumption can hardly be in 
strict accord with fact, and a divergence between 
the predicted and the experimental variation of 
the acoustical properties of the gas with pressure 
must therefore be expected. The interpretation 
of this discrepancy is touched on in Section 5. 
On the assumption (2.7) the pressure displace- 
ment yields 


(O(for1Ano—fioAn 1)/ON) tno AN 
=(AN/N)(fono—fiom1) =0. (2.6d) 


The displaced equation follows from (2.6a, b, c 
and d) and gives the result 
AT 
) (2.8) 


1+e-7+iw/fi0 


and the constants for substitution in (1.14) and 
(1.15) are evidently a=e-*, b=6, c=(1+e7*), 
d=w/fi, m= (1+e-*)—. The velocity of sound 
of small amplitude as a function of frequency 
will therefore be given by 


P Co(fiotfor)? +C,,w? 
Viu=—|14R 
Co?(fro+fo1)?+ C0 


p 

and this expression will also hold for the pressure 
dependence of the velocity insofar as (2.7) is 
valid. It is not yet suitable to describe the effect 
of temperature, since the relation between fio 
and temperature has not been examined. The 
absorption of sound is correspondingly 


h=2ny[(Q?—1)/(1+Q%y*)] (2.10) 


to the first order of small quantities if Q= V./ Vo 
and y = C,,w/Co(fio+fo). These are essentially the 
expressions given by Kneser, Bourgin, and several 
others. Their properties and implications have 
been widely discussed." 





EAn, = atm 





(2.9) 


3. A Four-STaTE GAs 


The majority of gases which have been 
investigated do not conform to (2.9) within 0.1 


11 Tn addition to the references already given the papers 
of Luck, Phys. Rev. 40, 440 (1932), and of Henry, Proc. 
Camb. Phil. Soc. 28, 249 (1932), should be cited, since 
they deal with aspects of acoustical absorption. 
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percent. For carbon dioxide, the most popular 
gas, Cy at 30°C is 3.435 R cal. deg.—', whereas 
C..+2Rx?no[e-*/(1+e-*) ] is 3.291 R® and Vo cate. 
is greater than Vo obs. by about 0.5 percent. In 
these cases it is, of course, possible to use (2.9) as 
an approximation, but the meaning of fo+/fio is 
then obscured. More detailed analysis of 7/A is 
needed to make available additional information. 
This may be exemplified by a model containing 
four internal energy states, E,;, Ey, and Ejy 
being as before, and Eyy being Ey plus the 
energy necessary to raise one gram molecule of 
the gas from the first to the second excited 
state of the slowly adjusted vibration. 
The equation for adiabatic change is then 


— PAv=C,,AT+ E(Ani+2Anz) (3.1) 


if the oscillator which represents the slowly 
adjusted vibration is harmonic (i.e., Eyjy— Ey, 
= Ey—Ey,=£). Two reaction equations are 
now necessary 


—dny/dt = (for +fo2)n0 — fiom —foone, 


(3.2) 
—dn,/dt = (fiotfie)n1—foino— fone. 


Two assumptions will be made in displacing 
these from equilibrium. The first, which supposes 
fu=Nkou, fiz= Ni, etc., has been discussed at 
length in Section 2 under (2.7). The second 
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supposes foe and fo to be much less than the 
other transition probabilities, and terms con- 
taining them consequently to be negligible. In 
the interpretation of experiment it must be 
remembered that these assumptions may not 
conform to fact. The displaced equations, to- 
gether with the simultaneous equation Am+ An, 
+ Anz,=0, give solutions for Am; and An, in terms 
of the exponential factor x and the various 
transition probabilities. When Am and Am, are 
substituted in (3.1) the appropriate form of the 
characteristic acoustical equation follows at once 
from (1.10). This has the form of (1.13) above 
where m= (1+e-*+e-**)—! and 


a=e—*+4e**+¢-*, 
b=w(e**/fiot+e-*/fa1), 
c=1+e-*+e-2*—w?/fiofer, 
d=w(1+e-*)(1/fiot1/for). | 


L (3.3) 





The velocity and absorption of sound follow 
from appropriate substitution in (1.14) and 
(1.15); the resulting equations are so complicated 
that a glance can hardly suffice to reveal their 
significance, and it appears therefore purposeless 
to print them in full. A few of the reductions of 
the expression for the velocity may, however, 
be examined with profit. For example, when 
wXfio, for 











P R 
Vet=—{1+ (3.4) 
p Co+Rx*nol (e-*+4e—-** ++ €-**) /(1+-e-*+e-**) ] 
which represents the terminal velocity at low ys _P R - 
frequency. When w>f;0, foi =" an C4 Rene (1++e-2)) ” 


V..® = (P/p) {1+R/C,,} (3.5) 


gives the terminal velocity at the highest 
frequency for which the model is adapted. This 
is of course not infinite; the terms in C., are 
still considered in equilibrium in all parts of the 
wave. When fa<w<fio an intermediate constant 
velocity 





12 Both these values have been calculated on the basis 
of fundamental frequencies 667.5, 1285, and 2326 cm™, 
as reported by Villars, Chem. Revs. 11, 408 (1932). 
C., is taken as 2.5085R cal. deg. in the second. 


is reached, and when fu<w< foi 
P R 

Vino —| f- ! (3.7) 
p C.+Rx*no(e~**/(1+e-*)) 





the alternate intermediate constant velocity, is 
obtained. A third intermediate constant velocity, 
due to the complete failure of transitions to and 
from the 1st state, may of course be imagined, 
but is considered so unlikely to be realized in 
practice that it has not been provided for in the 
four-state model above. It may simply be 
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evaluated by taking 
On2/dT = — O(noe~**) /OT 


which gives 





V trex Sis20 
P R 
ll PO }. (3.8) 
p C.,+4Rxono(e-27/(1+e-27)) 


The practical realization of constant velocities 
corresponding to (3.6), (3.7) and (3.8) does not 
appear very probable, but the influence of a 
decided difference between fio and fo: may be 
detected by experiment. 
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Fic. 1. Diagrammatic representation of the frequency- 
dependent part of the heat capacity for the five-state 
model. Transitions involving energy changes greater 
than +h have been considered of a low order of proba- 
bility, and harmonicity of the oscillator has been postulated 
by setting Eo: = Ej2= E23; = hyp. Where these approximations 
are permissible the model may be applied to any mode of 
vibration, whether transverse or longitudinal. 





4. A Five-STATE GAS 


Although the four-state gas model developed in the previous section is adequate to describe the 
behavior of carbon dioxide at all ordinary temperatures, it cannot be used for gases in which the 
fundamental vibrational frequency is much lower. It is exactly in these gases, however, that differ- 
ences in the transition probabilities of the various states may most easily be investigated. It is 
advisable, therefore, to carry the development of 7/A one step further, and examine a five-state 
model. As before the energy states E,;;— FE, will be considered in equilibrium with the adiabatic 
cycle at all frequencies. The frequency-variable transitions have been illustrated in Fig. 1, which, 
with the explanations in the previous sections, gives the model sufficient definition. The development 
has been carried out exactly as in the previous section with the help of the assumptions that the 
probabilities are directly proportional to pressure, that all transitions involving an energy change 
greater than +/Ay are so rare as to be negligible, and that the oscillator is harmonic. Three appro- 
priate reaction equations are displaced from equilibrium, and yield, with the equation Am+Am 
+ An.+ Anz; =0 

Anye*—xnAT/T ) 


twe*/fio+1 
An = (iw/ for +1) Anze* + (iw/ fer) Anze* — xnoe~*(AT/T) 
Ane => (tw/fse2+ 1) Anse* —xnoe~?*(AT/T) 





Any= 


b. (4.1) 





These are solved simultaneously and substituted in the equation for adiabatic change 
— PAv =C,AT+E(Ani+2An,+3Ans3) (4.2) 
to evaluate the constants for substitution in (1.14) and, (1.15). These are m= (1+e-*+e-*+e-*)-! 
a=e-*+ 4e-**+ 10e~87 + 4e-47 + 697 — w?(e3*/frofer te **/firofse te */farfse) 
b=wl_ (e-*+4e-* + e- 82) (e-*/f 9+ 1/fa2) + (e-*# +672 + 6-3 + €-4*) / fos | 
c=1+e-*+ e727 + e-82 — w?(1 +e-*) (1/fiofert1/firofse+1/forfse) 
d=w(1+e-*+e-**) (1/fro+1/fert+1/fs2) +we-*/foi—w?/froferfse 


Only two of the reductions of the resulting expression for the velocity need be examined. When 


-. (4.3) 
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i 
ve=—| 1 


: | 
C+ Rx? (e-*+4e-** + 10€-#* + 4e-47 + €-7) /(1+e-*+e-2* + e-3)2] 


(4.4) 


and the low frequency terminal velocity is obtained. The high frequency terminal velocity is identical 


in form with (3.5). If fox<w< fio, for 





a 
V? 55020 = —| 1 + 
p 


R 
C..+Rx?nol (e-* + 4e-27 + e- 37) /(1 pore ? 


(4.5) 


which represents a wholly different situation from (3.4) although the two expressions bear a strong 
superficial resemblance to each other. The other intermediate constant velocities are easily obtained 
by suitable modifications of the constants given in (4.3). They will not be printed. 


5. THE BEHAVIOR OF A GAS WITH MANY INDE- 
PENDENTLY ADJUSTED ENERGY STATES 


Before proceeding further it will be well to 
discuss the experimental bearing of what has 
already been written. This will be undertaken 
from two viewpoints: first, to discover what 
fresh information has become available and 
second, to find reasonably satisfactory approxi- 
mations which obviate, at any rate during the 
course of an experimental research, the necessity 
of calculating from the cumbersome equations 
resulting from (1.14) and (1.15) with appropriate 
constants. Since a relation between absorption 
and velocity is already at hand in (1.19), the 
velocity of sound will alone be given attention 
below. 

Figs. 2 and 3 illustrate the pressure and 
frequency dependence of the models developed 
in Sections 3 and 4, respectively. In each case 
the ordinates are expressed as the “apparent 
heat capacity ratio’? which is an abbreviation 
for the apparent ratio of heat capacity at 
constant pressure to that at constant volume at 
a particular value of w. The abscissae scale 
corresponding values of w/fo: which have, how- 
ever, been expressed as w/ No: since the foregoing 
exposition is committed to the formal equality 
of these quantities. This method of plotting 
describes, therefore, with equal suitability, the 
behavior of the models at constant frequency 
and variable pressure, or at constant pressure 
and variable frequency. Fig. 1 has been obtained 
by calculating numerically, from the constants 
(3.3) the cases fo = for, 10f 10 = for, and 100fi0 = fo 
in a gas for which the exponent hy»/kT is 3.00. 





Fig. 3 gives numerical values for the effective 
heat capacity ratio according to (4.3) in a gas 
where hy/kT is 2.00 for the cases fio=for=f0, 
100f 10 = 10fo1 = fse and 10*fi9 = 107fo1 = fz. To avoid 
ambiguity the slowly adjusted vibration has 
been considered to have an a priori probability 
of unity, although this will not be the case when 
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Fic. 2. The acoustical behavior of a model in which the 
frequency-dependent part of the heat capacity is adequately 
described by the 0—1 and 1—2 transitions. The ordinates 
represent the apparent values of C,/C, and consequently 
indicate also the variations of V?. The abscissae give cor- 
responding values for log w/foi, and serve therefore to 
describe the behavior of the model at constant pressure 
and variable frequency or, conversely, at constant fre- 
quency and variable pressure if the assumption fo = Nhu, 
which has been discussed in the text, is valid. The value of 
hvo/kT has been chosen as 3.0 for the purpose of plotting, 
and it has been assumed that the frequency-dependent 
part of the heat capacity is due to a single valence vibra- 
tion, the contribution of all other vibrations being neg- 
ligible. 
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Fic. 3. The acoustical behavior of a model in which the 
frequency-dependent part of the heat capacity is ade- 
quately described by the 0—1, 1—2, and 2=3 transitions 
of a valence vibration. In the case plotted hyo/kT has been 
given a value of 2. (This is suitable for comparison with 
carbon bisulfide at ordinary temperatures if the value of 
Co—C,, is doubled throughout.) The graph indicates that 
differences in the probabilities fo, and fiz of one order of 
magnitude should be detectable by measurements at 
constant frequency and variable pressure or at constant 
pressure and variable frequency. 


the dispersive region is caused by the failure at 
high frequencies of a degenerate vibration. C,, 
has been uniformly assigned the arbitrary value 
2.5R. 


The dependence of the velocity of sound on 
frequency at constant pressure and temper- 
ature 

Several conclusions follow from the consider- 
ation of Figs. 2 and 3. ; 

I. If the transitions 10, 21, 3-2 --+ have 
the same probabilities, the square of the meas- 
ured velocity of sound plotted as a function of 
the logarithm of the frequency will conform 
within experimental accuracy to a symmetrical 
S-shaped curve. An arbitrary number of energy 
states cannot of course be treated by the method 
used above, but there appears to be no danger 
in making this statement general. 

II. Any detectable deviation from a sym- 
metrical S-shaped curve may be attributed to a 
difference in the transition probabilities fio, fer 
‘++. Such a difference has been predicted by 


Zener,!* and its direct experimental demonstra- 
tion should not be difficult by this means. If the 
higher energy states are more rapidly adjusted 
than the lower, the high frequency part of the 
S-curve is flattened; if they are less rapidly 
adjusted than the lower, the curve is flattened 
in its low frequency region. 

III. If the difference between the rates of 
adjustment of any two states is greater than 
about two orders of magnitude, the deformed 
S-curve becomes irregular, and two or more 
maxima appear in the variation of the absorption 
coefficient with frequency. The appearance and 
separation of these maxima are conditioned by 
the relative importance of the two energy states 
in the heat capacity. 

IV. If the rate of adjustment of any state is 
greater than the others by about three orders of 
magnitude, its contribution to the heat capacity 
may, from the point of view of experiment, 
legitimately be incorporated in C,, and not 
assigned a specific rate of adjustment. This is 
the justification for giving 0(£;,—£;,)/dT its 
equilibrium value throughout the dispersive 
region in the proceeding sections. 

V. If the transitions 20, 3-0, 3-51, etc., 
involving energy changes of two or more quanta 
are not less frequent than those 1-0, 2-1, etc., 
by about two orders of magnitude the low 
frequency portion of the S-curve may be de- 
formed, and the only legitimate remedy is to 
describe the situation by a model containing a 
greater number of transitions. Since, however, 
the maximum possible contribution to the heat 
capacity of these transitions is small, their 
disturbing effect will probably escape experi- 
mental detection. The uncongenial task of still 
further complicating the calculations in this 
way will be postponed in the fervent hope that 
experiment may not demand it. 


The dependence of the velocity of sound on 
pressure at constant frequency and temper- 
ature 

A purely formal consideration of the relation- 
ships given in Sections 2, 3 and 4 leads to the 
conclusion that the increase of the frequency by 

a certain factor should be exactly balanced by 


13 Zener, Phys. Rev. 37, 556 (1931); 38, 277 (1931). 
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the decrease of the pressure by the same factor. 
All the conclusions just stated for the dependence 
of the velocity of sound on frequency therefore 
appear to hold in the inverse sense for its 
dependence on pressure. The correspondence 
between frequency and pressure is due, however, 
to the simplifying assumptions and approxi- 
mations which have been used. In practice the 
conclusions reached from considering frequency 
may be applied to pressure only with the 
following important qualifications. 

I. In a real gas the heat capacity at constant 
volume is not independent of pressure, and 
pressure and numerical density are not exactly 
proportional. Although it is always possible to 
eliminate P/p by dividing V, by Vo at any given 
pressure or frequency, errors due to these causes 
may become experimentally manifest at high 
pressures. 

II. The assumption that fio, for, etc., vary 
directly in proportion to the pressure may be 
invalid on several grounds. If, for some reason 
which is not at present clear, triple collisions are 
more effective thai double collisions in adjusting 
the internal energy fio~ Nkiw. The mathematical 
treatment will not seriously be disturbed by this 
since (2.6d) will vanish for any reasonable 
functional relationship between fio, etc., and N. 
A diminution of the pressure will, however, 
have more effect on the velocity and absorption 
than a corresponding increase of frequency. If 
the advantage of triple collisions is marked, this 
effect should be very strikingly manifest. If only 
slight, it may escape experimental detection 
because of the relative infrequency of triple 
collisions. A correction, if necessary, may be 
applied to the high-pressure region by expressing 
each f as the sum of two probabilities, one of 
which is proportional to the number of double 
collisions, and one to the number of triple 
collisions. 

III. Again, it is possible that a molecule which 
is vibrating has an advantage over one which is 
not vibrating in bringing about energy transfor- 
mations. Since no example of this is yet known, 
we need not concern ourselves with a mechanism 
for it, but merely with its effect on the acoustical 
behavior of the gas. A more complete statement 
of the assumption fis=RiN reads fio= kiN 
+k'iNm+k''Nne+-++ where Rio=R' p=’ 19 
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-++. If this identity is for any reason invalid, 
Ofio/dn, no longer vanishes, and (2.6b) is in- 
correct. A variation in m may, in a perfect gas, 
be brought about only by changing the temper- 
ature, and this introduces a set of approximations 
which have not yet been considered. Both for 
this reason, and because 7, m2, etc., are usually 
small fractions, it seems unlikely that any 
inequality of ki, k’1, etc., may at present be 
demonstrated by experiment in pure gases. 

The presence of another mode of vibration 
which, although excited to an appreciable extent, 
remains in equilibrium with the sound wave 
throughout the dispersive region, may, of course, 
cause exactly similar effects. There is some 
reason to believe, however, that resonance be- 
tween the two vibrational modes must be close 
in order for collisions of this type to be effective. 

IV. When the dispersive region is produced by 
the failure of transitions between vibrational 
states which are optically active, the excited 
molecules may radiate their energy between 
collisions. The effect of this will be to cause the 
mean lifetime of the molecules to increase less 
rapidly with decreasing pressure than is indicated 
by the purely kinetic assumption fi9 = Ni. This 
possibility was first envisaged by Bourgin,"* who 
has pointed out that the effect becomes impor- 
tant at low pressures. It has recently been 
stressed, especially in connection with acoustical 
absorption, by Kneser.'® It is not intended to 
extend the scope of this study to cover radia- 
tion effects. It will therefore merely be pointed 
out that if the relaxation time is written as 
8=1/[N(Rkiot+ko)+1/7] where 7 is a tempera- 
ture and pressure-independent quantity express- 
ing the radiation lifetime, a rough estimate 
of + may be obtained from measurements of 
the velocity of sound as a function of pressure 
in the dispersive region. Since the temperature 
and pressure derivatives of 7 vanish, its intro- 
duction in this way does not disturb the exposi- 
tion given in Sections 2, 3, and 4. The radiation 
lifetimes of molecules are probably so great that 
it is not to be expected that radiation effects will 
become appreciable at pressures which are at 
present suitable for the measurement of the 


44 Bourgin, reference 9, 1928 papers. 
15 Kneser, reference 6, especially 1933 paper. 
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velocity of sound. It should be noted, however, 
that acoustical measurements at very low pres- 
sures may ultimately become an effective means 
of investigating radiation lifetimes. 


The dependence of the velocity of sound on 
temperature at constant frequency and 
numerical density 


The numerical density of a gas increases with 
decreasing temperature at constant pressure, and 
it will therefore be convenient to discuss the 
variation of the velocity of sound with temper- 
ature at constant numerical density rather than 
at constant pressure. Some assumption must 
now be made concerning the temperature de- 
pendence of C,, and the f probabilities. 

C., will depend on temperature if some of the 
non-lagging degrees of freedom are not fully 
excited.'® In this case C,,=$R+,Rx;°e-1/(1 
—e1)?4---+ if p, etc., are the a priori proba- 
bilities of the degrees of freedom, and x, etc., 
abbreviations for hv;/kT, etc. 

In discussing the temperature coefficient of the 
f probabilities only those of the type fio, fo, etc., 
will be considered, since the relation between 
these and the fo, fiz, etc., probabilities is already 
known. The simplest assumption is the C,, and 
fio, etc., are independent of temperature. In this 
case the three, four, and five-energy state 
models behave similarly: the S-curve is merely 
compressed as the temperature is lowered. This 
is roughly described by stating that (V,?— V.,”) 
/(Vo?— V2) is independent of temperature for 
a given value of w. Fig. 4 diagrammatically 
illustrates this situation. 

Such measurements as exist on the tempera- 
ture dependence” of the velocity of sound in the 
dispersive region indicate that the transition 
probabilities in several gases are by no means 
independent of temperature. 

To account for this a quantity called the 


© Kneser and Ziihlke, Zeits. f. Physik 77, 649 (1932), 
have shown in nitrous oxide that the deformation vibration 
fails at frequencies lower than the other vibrations. 
Similar evidence concerning carbon disulphide has recently 
been obtained in the writers’ laboratory in Princeton by 
Mr. J. A. Reid. 

‘7 Apparently the measurements of Richards and Reid, 
Nature 130, 739 (1932), are the only ones reported present 
which deal directly with this question. 
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Fic. 4. Diagrammatic representation of various possi- 
bilities for the temperature dependence of the velocity of 
sound at constant frequency and constant numerical 
density. In practice fio will depend on the number of 
collisions per unit time, and the straight line ‘‘fio inde- 
pendent of T”’ is therefore an idealization. Relatively small 
differences between the activation energies of collision of 
two states should be experimentally manifest by measuring 
the temperature coefficient of the velocity at the low 
frequency and high frequency limits of the dispersive 
region. 


activation energy of collision has already been 
suggested.'” This bears an analogy to the activa- 
tion energy of a chemical reaction in that 


Afio fivAro @for _for(A1o+ £) 
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etc. 


and represents the ‘‘threshold kinetic energy” 
necessary to effect the transitions 1-0, etc. 
Such an impediment to the transference of 
vibrational to kinetic energy might have several 
causes. For example Heil,'* following Oldenberg,'® 
has suggested that the dynamics of collisions of 
the second kind demand that in some cases only 
molecules of particularly high kinetic energy 
should be effective. On the other hand, a po- 
tential energy barrier of the same general type 
that found in chemical reactions might be 
imagined.'’: 2° This possibility will be discussed 
at slightly greater length in order to make clear 


18 Heil, Zeits. f. Physik 14, 31 (1932). 

19 Oldenberg, Phys. Rev. 37, 194 (1931). 

20Franck and Eucken, Zeits. f. physik. Chemie A20, 
460 (1933). 
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certain characteristics which the activation en- 
ergy of collision may conceivably possess. 

If two molecules approach one another with 
a given relative kinetic energy they influence 
increasingly each other’s chemical binding en- 
ergies as they come increasingly into each other’s 
repulsive fields. This process will, with a certain 
probability, induce a change in the vibrational 
states of the colliding molecules. A simple though 
approximate method of calculating the transition 
probability has been given by Kailmann and 
London.”! It may be assumed, in analogy to 
their mechanism, that the process of approach of 
the two molecules occupies an extremely short 
period of time, and that during this time the 
vibrational wave function remains unchanged. 
This means that, owing to the mutual “‘soften- 
ing” of the chemical bonds,”"* there is a certain 
probability that a change will occur in the vibra- 
tional quantum numbers of one or both of the 
colliding molecules. The process of separation is, 
conversely, considered to be very slow, so that 
no further change of the vibrational quantum 
numbers occurs; this must be assumed in order 
to have a nonvanishing transition probability, 
since in a rapid separation the wave function 
would again remain unchanged, and the vibra- 
tional quantum numbers would revert to their 
original values. 

During the second, or adiabatic, process the 
hardening of the bonds as the molecules separate 
will effect an increase of the vibrational energy. 
If vy represents the frequency of the original bond, 
and v’ that of the softened bond, the initial vibra- 
tional energy may be increased by a factor, v/v’. 
For a molecule originally in the mth state this 
excess of energy, which is supplied at the expense 
of the relative translational energy, will amount 
to (v/v’)(n+1/2)hy—(n+1/2)hv, and the prob- 
ability of an excitation will be roughly in propor- 
tion to this divided by hy, or to (v—v’)/v(n+1/2). 
This quantity increases rapidly with increasing 
values of m, and is three times larger for the 
second than for the first vibrational state. Thus 


21 Kallmann and London, Zeits. f. physik. Chemie B2, 
207 (1929). 

218 The “‘softening’’ of the chemical bonds as a forerunner 
of chemical decomposition has been studied, in connection 
with Raman spectra, by Cremer and Polanyi, Zeits. f. 
physik. Chemie, Bodenstein-Festband, 771 (1931). 


when n is.small (v—v’)/y must be relatively 
larger, and a greater softening of the bond is 
necessary to effect an increase of vibrational 
energy. This can be brought about only by closer 
approach of the molecules, and consequently by 
collisions involving greater relative kinetic en- 
ergies. The difference will be experimentally 
manifest by larger activation energies of collision 
for the excitation of the lower than of the higher 
quantum states. 

While the mechanism described above is in- 
correct in that the rapidity of approach and 
separation must, in actuality, be of the same 
order of magnitude, it is believed to represent a 
plausible picture for the expected relationship 
Ay>An>Az..... It is hoped that a more 
satisfactory description of energy transitions 
during collision of the first and second kinds may 
be obtained from a quite different method of 
approach, which cannot here be included. 

It is not proposed, however, to introduce the 
activation energy of collision on deductive 
grounds, but to use it as a purely empirical 
quantity calculable from the temperature co- 
efficients of fio, fo1, etc., according to the expres- 
sions derived from (1.14) and (1.15) in appro- 
priate cases. Any reasoning from analogy with 
the chemical activation energy will be sup- 
pressed, since there is at present no convincing 
reason for believing the activation energy of 
collision to be independent of temperature over 
any considerable range. Furthermore, it is not 
at present possible to predict its order of magni- 
tude. It has been introduced only to facilitate 
the description of the temperature coefficient of 
the velocity in terms of experimentally deter- 
mined quantities, and to permit a rough kinetic 
analysis of the probabilities fio, foi, etc. 

Before proceeding further it must be pointed 
out that the existence of an activation energy of 
collision introduces a new approximation into 
the treatment of the previous sections. If the 
molecules on the high energy side of the Maxwell- 
Boltzmann distribution lose translational energy 
to vibrational energy on every collision, the 
effective pressure of the gas must also vary 
somewhat with frequency in the dispersive 
region. It does not appear possible at present to 
compute the extent of this approximation, but 
it is believed to be of small importance. 
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The effect of an activation energy of collision 
on the dispersive behavior of the gas is clear. 
If fio, for, etc., all have equal activation energies, 
lowering the temperature will correspond to 
increasing the frequency, and the gas may be 
made to pass through its dispersive region by 
means of temperature changes alone. This case 
has been diagrammatically represented in Fig. 4; 
the S-curve is slightly deformed by the variation 
of Cy with temperature. The greater the activa- 
tion energy the more rapidly will the gas pass 
through its dispersive region as the temperature 
is lowered. If the activation energies are different 
for fio, for, etc., the S-curve will be more greatly 
deformed. Since it seems probable, on the 
argument given above, that A1> Ado, etc., this 
case has also been diagrammatically represented 
in Fig. 4. The extreme upper and lower temper- 
ature of the dispersive region will in this case 
show different activation energies of collision, 
which may be interpreted by means of a gas 
model containing a suitable number of energy 
states. 


The kinetic analysis of transition probabilities 


All modern writers on the dispersion and 
absorption of sound have used a kinetic mech- 
anism as the basis for their interpretation.” 
If experiment has established the magnitude of 
the activation energy of collision and the 
essential importance of double collisions in the 
excitation process, a little further information 
may be obtained by subjecting the transition 
probabilities to a rough kinetic analysis.'7? The 
results so obtained must not, however, be 
considered significant in any exact quantitative 
sense. 

The number of double collisions suffered by 
the molecules in the 0 state per cubic centimeter 
per second is approximately 


Z=4N2no?(eRT/M)}, (5.1) 


where N is the total number cc~' and m the 
fraction of these in the 0 state as before, M is 
the molecular weight, and o the distance of 
nearest approach of the molecules in a double 





*? Apparently Jeans (Dynamical Theory of Gases, Cam- 
bridge (1904)), was the first to take this standpoint. 
Einstein, Herzfeld and Rice, Bourgin, and others have 
Suggested or stressed a kinetic mechanism. 
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collision. On the assumption that o is unaltered 
by the changes of state 0-1, 0-2, etc., similar 
expressions for Z;, Z2, etc., follow. Carrying 
through the previous assumption that either 
steric or energetic properties of the molecules 
may hinder the transformation of kinetic into 
vibrational energy, the number of collisions cc 
sec.~! which are effective in removing a molecule 
from the 0 to the 1 vibrational state may be 
written. 


20 =4Nonokoi0?(t RT /M)ie—40/®7, (5.2) 


where 27é 0? is a steric factor not strongly 
dependent on temperature, which may be termed 
the ‘‘effective cross section’’ of the molecule for 
the transition 0-1. Since the f probabilities are 
expressed per molecule per second it is apparent 
that to find fo: 2 must be divided by Nm, the 
number of molecules in the 0 state cc~'. Thus 


for — (4No?(xRT/M) *) Eye~ (Art E)/RT (5.3) 
and 
fio = (4No2(eRT/M)})E e427, (5.4) 


from which the equality of the effective cross 
sections in the two directions is immediately 
deducible from fo:/fio=e-*. Similar relations 
between any pair of transition probabilities 
Fn, n+tfn4i, n May similarly be written. There is 
no reason at present for supposing the various 
€ values to be equal for different states. 

From the comparison of 2707& with the ordi- 
nary kinetic cross section deduced from viscosity, 
effusion, and the like, it is possible that inter- 
esting information concerning the spatial char- 
acteristics of the excitation process may be 
deduced. For example if a triatomic colinear 
molecule is visualized by the naif model of three 
balls connected by springs, a ‘‘head on”’ collision 
of another molecule with either of the end 
atoms must be far less effective in exciting the 
deformation vibration than a ‘‘broadside’’ col- 
lision of another molecule with the center atom. 
The effective solid angle for collisions may in 
this case be small, and the ordinary cross 
section far greater than that indicated by 
acoustical measurements.”* 

It should be noted that if the exponential 


23 Heil, reference 18, has shown this from Kneser’s 
measurements on carbon dioxide. 
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factor in (5.4) is near to unity the effective cross 
section must be expressed with greater circum- 
spection than has here been considered necessary, 
since in this case it may dominate the tempera- 
ture coefficient of the dispersion. The proba- 
bilities fio, etc., must also be expected to di- 
minish slightly with temperature owing to the 
decreased frequency of collisions. 


Approximations suitable for orientation in ex- 
perimental investigations 

Since calculation from the four and five-state 
modifications of (1.14) and (1.15) is a laborious 
and uncongenial process, two approximations 
are here given which it is hoped will be useful 
in orienting the experimenter during the course 
of his measurements. They are in certain cases 
also suitable for the final expression of measure- 
ments. 

To describe the frequency and pressure de- 
pendence of the velocity the expression 

P ROo+RC,,w* dn? 

Vp =—t i+ 
p Cor+C,,20dy? 
where dy = (e7 +1) /(e7—1)*fo1, is very nearly iden- 
tical in its behavior to that of the multistate 
models of Sections 3 and 4 when fip =f21 =f32. The 
correspondence with the four state model is so 
close that no difference is apparent in Fig. 2. 
The correspondence with the five state model is 
within all reasonable experimental limits, as the 
dotted line in Fig. 3 shows. Any experimentally 
detectable deviation from (5.5) may therefore 
be attributed to inequality of fio, fo, etc., if the 
frequency dependence of the velocity is in 
question. Conversely if fio, fo1, etc., have been 
proved identical by varying the frequency at 
constant pressure and temperature, deviation 
from (5.5) may be attributed to radiation or 
triple collision effects if the pressure dependence 
of the velocity is receiving attention. 

It has been proposed to study the temperature 
coefficient of the velocity in the dispersive 
region merely to determine the activation energy 
of collision for each of the various excitation 
processes. Since this is given by fio/f’10, etc., 
corresponding to JT and 7”, an error in the 
absolute magnitude of fio, etc., is of no im- 
portance. It is somewhat better for the determi- 
nation of the activation energy of collision to 


(5.5) 
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use the approximation 
ROot+RC,,w?d 7? 
Coer+C,,202d 7? 





r (5.6) 


where #7 = 1/fio(1-+e-7) rather than the approxi- 
mation (5.5). The values of fio so obtained may 
be noticeably in error, but this will not affect 
the calculated activation energy. If Ao is differ- 
ent from Ags, etc., the low and high frequency 
parts of the S-curve will show different temper- 
ature dependence, as illustrated by Fig. 4. The 
exact significance of the difference can be demon- 
strated only by a more profound analysis, but 
its existence may conveniently be demonstrated 
by (5.6). 

Both (5.5) and (5.6) assume a knowledge of 
the energy levels of the molecule, since they 
require the value of e*. Cases may occur where 
these are not obtainable. A rough analysis may 
then be made by supposing Cy)— C,,”™ to be due to 
the failure of a single vibrational transition 
involving a suitable quantum jump. This will, 
of course, lead to more or less erroneous f values, 
but will give a variation of 3; with temperature 
which is approximately correct. 


6. A Few FurTHER CASES. GASEOUS MIXTURES 


It is believed that the suggestions given in the 
previous sections are sufficient to interpret the 
acoustical behavior of all pure gases which have 
at present been examined experimentally. They 
by no means suffice, however, to cover all 
possible cases of interest. These are, indeed, so 
numerous and diverse that the labor and space 
required to examine them individually are pro- 
hibitive. Since the description of a system of any 
necessary complexity by the method of Einstein 
or of Bourgin is a matter of labor rather than of 
ingenuity, more complicated cases may reason- 
ably be postponed until required by experiment. 
The purpose of this section is only to give a 
superficial survey of further territory. 


The simultaneous failure of two modes of 
vibration 

Kneser’s contention that in gases in which the 

deformation vibration is excited it is this vibra- 


24 This must in any case be known in order to identify 
the energy transitions which are absent at high frequencies 
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tion alone which causes the first dispersive 
region®® appears to be justified in all cases which 
have at present been reported. It has neverthe- 
less been noted that chlorine and oxygen, which 
can have only longitudinal vibrations, show 
dispersion of sound in the region of frequency 
which is at present accessible to experiment, and 
the simultaneous failure of two vibrational 
modes in the acoustical cycle must therefore be 
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considered as possible. If the excitation and 
de-excitation of the two modes is due only to the 
relative kinetic energy of the two molecules on 
collision, the behavior of the gas will be a 
composite of two examples of the type already 
discussed. An approximate description is ob- 
tained by writing for the characteristic acoustical 
equation 


Tv P R(1+iwd 1) (1+iwd2) 





+ , 
A p C..(1 + tad) (1 +282) + Ci(1 +e) + Co(1 +7001) 


in which, of course, Co=C,.+Ci+C: and 3; and 
3, are the approximate relaxation times of the 
vibrational modes k and /, respectively. Service- 
able if somewhat inexact expressions for velocity 
and absorption follow from (6.1) by (1.6) and 
(1.7). The absorption will show maxima at two 
frequencies if #, and #, differ by about two 
orders of magnitude. 

A more interesting case is provided if vibra- 
tional resonance phenomena increase the proba- 
bility of the exchange of energy between the 
two modes of vibration. Such a possibility has 
already been touched on in Section 5, paragraph 
III of the discussion dealing with the variation 
of the velocity of sound with pressure. Each 
transition may then be brought about by one of 
several means, and it will in general be the most 
effective of these which determines the position 
of the dispersive region. The probability fo: will, 
for example, be the sum of the number of times 
per second which a molecule leaves the 0 state 
for the first state due to each of the various 
types of collisions. It is clear that unless a 
considerable fraction of the molecules are excited 
the number of transitidns per second due to 
collisions with vibrating molecules will escape 
detection. It is also clear that, since the rate of 
transition due to each type of collision cannot 
be measured separately, a large number of 
extremely accurate measurements at different 
temperatures and pressures must be made in 
order to separate one from the other. It is 
therefore more likely that information of this 
type will be obtained from the study of gaseous 
mixtures than from pure gases. 





*° Kneser, reference 6, 1933 paper. 





(6.1) 





The behavior of a dissociating gas in which 
the heat capacity may vary with frequency has 
been treated at great length in Parts I and Ia 
of this series. It is perhaps worthy of note that 
if the conclusions stated above are justified this 
treatment is adequate for all present purposes. 


Binary gaseous mixtures 

The discussion of the previous paragraphs 
immediately suggests a parallel treatment for 
gaseous mixtures. This branch of acoustical 
theory has already been comprehensively de- 
veloped by Bourgin and nothing of importance 
need here be added to what he has said. The 
following remarks are primarily to be interpreted 
as a translation of Bourgin’s work into the 
language of Einstein’s theory. The analysis of 
the transition probabilities is, however, some- 
what more specific, and consequently more 
vulnerable than that which Bourgin has pro- 
vided, although it is believed to be fully in 
accord with his intentions. 

The simplest dispersive binary gaseous mixture 
consists of a dispersive gas A and a nondispersive 
gas B. The heat capacity in equilibrium of A is 
Coa mole deg.—!, and of this a part Coa—Caa 
is dependent on frequency. The heat capacity 
of B is Cog at all frequencies which are experi- 
mentally considered. The numerical density of 
A is mN and of B(i—m)N. The effective heat 
capacity at very low frequencies of sound will 
then be Co=Coam+Con(1—m) and at high fre- 
quencies will be C= C,am+Con(1—m). It will 
probably be sufficient to describe the acoustical 
behavior of A by a three-state approximation, 
particularly if measurements are compared at a 
single frequency in the dispersive region. In this 
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case the characteristic acoustical equation of the mixture is 


R(1i+iws) 6.2) 





; 
A pam+pp(1—m) 


and expressions for the velocity and absorption 
follow from (1.6) and (1.7). P in (6.2) is of 
course the sum of the partial pressures of A and 
B. The relaxation time #? is now, however, a 
more complicated quantity; collisions of an A 
molecule with another of the same species or 
with a B molecule may produce energy transi- 
tions, and it is the more effective of these 
encounters which will dominate the relaxation- 
time. This may be expressed by writing, for a 
pressure great enough to exclude radiation effects, 

3 =const./(faat+fas)(1+e*), (6.3) 
where faa and fas are the probabilities per 
molecule per second of removing vibrational 





fas =4N(1 —m)oasp*éan(tRT (M+ Mp)/MaMz) te-AA pitt 


the temperature variation of fag will permit a 
comparison of the steric and energetic effective- 
ness of collisions of the types AA and AB. 

The next simplest binary mixture consists of 
two gases A and B each of which has a dispersive 
region in the frequency range under observation. 
If direct transmission of energy from the vibra- 
tions of A to those of B is not possible, the 
characteristic acoustical equation of this mixture 
will be given by (6.1) if P=Pam+P2,(1—™m), 


Ci=(Coa—Cpa)m, C2=(Cos—C,n)(1—m), 
3; =const.4/(faat+fas)(ite* ) 


and #,=const.s/(fae+faa)(i+e-*8) are substi- 
tuted. In view of what has been said in con- 
nection with (6.1) and (6.2) this requires no 
further discussion. 

It should be pointed out, however, that this 
case provides the best means of investigating 
advantages due to vibrational resonance effects 
in collisions between unlike molecules. The 
probabilities faa and fzez may be separately 
determined for the two pure components, and 
the analysis of measurements on such a mixture 
will therefore be at a great advantage over 


C..(1 + tw?) + (Coa —Ca)m 





energy from an excited state by collisions of the 
type AA and AB, respectively. The constant in 
(6.3) may be taken as unity without losing the 
significance of the measurements; it will become 
larger than this if Coa—C.,4 is great. The values 
of faa are obtained from measurements on pure 
A at various temperatures and partial pressures. 
In this way the activation energy of collision 
and effective cross section for collisions of the 
type AA may be evaluated by means of (5.3). 
Measurement at a certain temperature, fre- 
quency, and composition determines the value 
of # in the binary mixture. The value of faa at 
a corresponding partial pressure of pure A is 
then substituted in (6.3) and faz is obtained. 
Since, roughly 


(6.4) 





similar measurements in pure gases. By suitably 
varying the number of molecules in the excited 
states of the two molecular species with temper- 
ature, there appears to be a reasonable chance 
that significant information concerning vibra- 
tional resonance may be obtainable. 


The behavior of a gas on illumination 

If a gas is illuminated by light of wave-length 
suitable to raise it to a state of electronic 
excitation, the velocity of sound may change. 
The vibrational levels in the excited state will 
be, of course, quite different from those in the 
normal state. On radiating the gas may be 
expected to drop back to those vibrational levels 
in the normal state which disturb the relative 
separation of the nuclei as little as possible, in 
accord with the Franck-Condon principle. These 
will in general be of higher energy than those 
excited by collision. An alteration of the heat 
capacity and a change in the velocity of sound 
may therefore be expected. The situation so 
created is, of course, exceedingly complicated, 
and does not encourage general treatment. The 
number of molecules raised to the electronic 
state will also be only a very small fraction of 
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the total since the absorption coefficient for light 
is customarily small in gases. Nevertheless it 
may be possible to find a system which shows 
the effect of illumination to a marked extent 
with strong light sources and, in addition, 
possesses a sufficiently simple energy structure 
to permit its quantitative study without pro- 
hibitive labor. The possibility is mentioned here 
since experiments which are open to this inter- 
pretation have already been reported in the 
nondispersive region of several gases.”° The study 
of a dispersive gas with a high optical absorption 
coefficient such as bromine might, in the dis- 
persive region, yield valuable results if the 
optical conditions were quantitatively controlled. 


28 See for example, Kupper, Ann. d. Physik 43, 905 
(1912); Trautz, Zeits. f. Elektrochemie 18, 513, 651 (1912). 
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* * * * * * * * * * * * kK Kk * 


It is hoped to make the next two parts of this 
series a detailed account of the experimental 
investigations on carbon dioxide, carbon di- 
sulphide, and binary mixtures of ethylene with 
other gases which have already been partially 
reported.!” 

It is a pleasure to acknowledge the assistance 
given throughout this study by Dr. Max Del- 
briick of the Kaiser Wilhelm-Institut fiir Radio- 
chemie. Not only has he helped me to understand 
what Bourgin has written by correcting the 
immense number of printer’s errors which occur 
in the work of this writer, but he has also 
frequently helped me to understand what I have 
written myself. 
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According to previous experiments liquid helium shows 
a ‘transition point” with discontinuity in the values of its 
specific heat and expansion coefficient. It has not been 
possible, however, to observe the coexistence of the two 
“‘phases.”’ This fact is here explained on a thermodynamic 


basis. Coexistence of two phases at a transition point can 
only be anticipated if the two phases exhibit a discontinuity 
in energy. An analogy is made with similar cases like that 
of the ‘transition point of superconductivity.”’ 





XPERIMENTS! have recently led to the 

surprising result that liquid helium shows 
a ‘transition point’ at extremely low tempera- 
tures. This transition point distinguishes itself 
from those with which we are usually acquainted 
by the fact that the two phases of helium exhibit 
no discontinuity as far as energy content or 
specific volume is concerned. They distinguish 
themselves from each other essentially by a jump 
in the values of the expansion coefficient (dv/dT) p 
and of the specific heat c,. It has been found 
experimentally that this transition point also 
depends upon the pressure so that we should 
speak more strictly of a transition curve in the 
p—T diagram. Furthermore, relations have been 
found to exist between the values for the dis- 
continuity of the specific heat and the expansion 
coefficient. These empirical relations have been 
derived theoretically as a consequence of the 
second law by Keesom! and Ehrenfest.? We 
reproduce briefly the line of reasoning developed 
by Ehrenfest, since it will help to disclose the 
main theme of this short note more easily. 

Let us denote by Z the thermodynamic 
potential of each helium, phase; Z is a function 
of » and 7. Ehrenfest considers the two modifi- 
cations of helium as two different phases which 
when in equilibrium are characterized by the 
condition Z;=Z»2. If we differentiate Z with 
respect to p and T, we obtain, 0Z/dp=v and 
0Z/dT = —p. Since there is no discontinuity in 
volume or energy (and therefore entropy) the 


1For references cf. W. H. Keesom, Proc. Amsterdam 
Acad. Sci., Sect. A 36, 147 (1933). 

2? Ehrenfest, Proc. Amsterdam Acad. Sci., Sect. A 36, 
153 (1933). 


first differential coefficients of the thermody- 
namic potential are continuous at the transition 
point. This does not hold true for the second 
derivatives. 


Z/dT*? = —0/dT = —c,/T, 
&Z/aT ap =dv/aT, 
3Z/dp? = dv/dp. 


Since as mentioned above the values for the 
entropies do not show a discontinuity at any 
point of the transition curve, we may take the 
total differential of both sides of the equation 
S;=S_ provided the differentials dp and dT are 
so chosen that we always remain on the transition 
curve. These special values of dp and dT shall 
be denoted by Dp and DT, following Ehrenfest. 
We thus obtain 


— (0 /dT),Dpt (¢m/T)DT 
= — (002/8T) »Dpt (¢y2/T)DT 


or 


Dp/DT = (Cy: —Cy2) /T(001/0T — 002/dT), 


a relation which has been established experi- 
mentally by Keesom. Several other similar 
relations can be derived which need not be 
discussed at present. 

The question which has been raised by Ehren- 
fest, and which we shall try to answer is the 
following: Though we have treated liquid helium 
on both sides of the transition point as two 
different phases in a thermodynamical sense, 
and though the results derived from the equi- 
librium condition for the two phases are correct 
and in agreement with experiment, the observa- 
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tion of the coexistence of the two phases at the 
transition point as is usually possible at all other 
transition points has not been accomplished. 
What is the nature of the distinction between 
the helium transition point and the commonly 
observed transition points? 

This question is readily answered if, for the 
sake of lucidity, we take as the independent 
variables of our thermodynamic system energy 
and volume instead of temperature and pressure. 
It then becomes obvious that a coexistence of 
two phases is only stable, i.e., observable, if these 
two phases show a discontinuity in energy (the 
first differential coefficient of the thermodynamic 
potential) and not only in specific heat (the second 
differential coefficient of the thermodynamic 
potential). 

This will be clear from the following consider- 
ations. Let us plot, for example, the melting 
curve in the energy-volume diagram. Each point 
in this curve represents a certain equilibrium 
state with a definite ratio between the respective 
amounts of molten and solid substance present. 
To an infinitesimal variation of the energy there 
corresponds an infinitesimal displacement of the 
point on this curve or an infinitesimal change in 
this ratio. Temperature does not change unless 


we have reached the terminals of the curve at 
which points the whole substance is either solid 
or molten. 

In the case of the helium transition point the 
situation is essentially different. Here all ratios 
of the composition of helium I and II correspond 
to one value of the energy. An infinitesimal 
change in energy would immediately eliminate 
one phase, and raise the temperature of the 
substance. This would hold true for every point 
in the transition curve, and not merely for the 
terminals, since there is no energy difference 
between the two phases. The coexistence of the 
two phases therefore never would be stable 
with respect to infinitesimal fluctuations in 
energy. This explains why the two phases could 
not be observed simultaneously. 

It has already been pointed out by Ehrenfest 
that the discontinuity in the specific heat of 
helium shows much similarity to the discon- 
tinuity which is observable at the temperature 
at which a metal exhibits superconductivity. 
Here too it has not been possible to observe 
the coexistence of the two ‘‘metal phases.’’ We 
can infer from the preceding line of reasoning 
that, at the transition point of superconductivity, 
no appreciable energy change occurs. 
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On Electrophoretic Mobilities and the Isoelectric Point of Protein-Coated Particles 


Otto HALPERN, Physics Department, New York University, University Heights 
(Received October 14, 1933) 


It is a fact known by experiment that quartz or glass 
particles of microscopic size suspended in a liquid and 
coated with a superficial layer of protein exhibit in almost 
every instance the same mobility in response to the applica- 
tion of an electric field. This indicates that their electro- 
phoretic potential is independent of the size and shape of 


the particles. Furthermore, the isoelectric point of these 
particles coincides with the isoelectric point of the 
proteinic ions in the solution. Some very general thermo- 
dynamic arguments are given in this paper indicating 
that we should expect a similar behavior as a consequence 
of simple model properties. 





HE difference between the electric potential 
as it enters into the formulae of thermo- 
dynamics, and the electrophoretic potential as 
determined by kinetic experiments (experiments 
of flow in general) is now clearly recognized.! 
While the first refers to the electrical work done 
in the transference of a charged particle from 
one phase to another, and therefore can be due 
to excess charges in one phase as well as to 
double layers rigidly connected with one phase, 
an electrokinetic potential, in the second place, 
can only be observed provided one of the phases, 
which are in motion relatively to each other show 
an excess charge. Keeping this in mind, it is 
not difficult to understand that in extreme cases 
the two potentials might change in an opposite 
sense, with the variation of the ionic concentra- 
tion, as has been found experimentally by 
Freundlich,? and interpreted later by Stern.* 
From the preceding remarks it should not be 
expected that it would be possible to draw any 
inferences relative to electrokinetic potentials 
from thermodynamics. This is generally true, 
although in the case at present under considera- 
tion it seems as if the combination of a simple 
assumption touching the model properties of the 
protein with general thermodynamic arguments 
permits of an understanding of experimentally 
discovered uniformities of behavior. 


1 For a treatment of the theory and experiments pertain- 
ing to electrokinetic phenomena, see H. A. Abramson, 
Electrokinetic Phenomena, Chem. Cat. Co., New York. 
(In Press.) 

2H. Freundlich, Kapillarchemie, Leipzig, Vol. I, 1932. 

30. Stern, Zeits. f. Elektrochemie 30, 508 (1924). 


The theory of electrokinetic phenomena as 
developed by Helmholtz, Smoluchowski, and 
others shows that the electrokinetic potential ¢ 
of a particle is connected with its mobility u 
under the influence of an electric field by the 
relation 


¢=(Cn/D)u (1) 


(n=coefficient of friction, D=dielectric con- 
stant). The numerical factor C is equal to 4z for 
particles the radius of which is large as compared 
to the thickness of the ionic atmosphere sur- 
rounding the particles. 

Proteins dissolved in a solution of known pH 
act as amphoteric electrolytes; they are able to 
ionize by parting with either H*+ or OH™ ions. 
The dissociation equilibria are given by the 
relations 


(CyCon-)/C= K,, (2) 
(C_Cy+)/C=K_, (3) 


where the notation is self-explanatory. 
Combining the relations (2) and (3) with the 
expression for the solubility product of water 


Cy+Con-=K (4) 
we obtain as a condition for C, = C_ the relation 
Cu+=[(K_/K4)K ]}. (5) 


The pH for which C,=CL- is called the iso- 
electric point. At this isoelectric point the 
protein particies do not show any mobility in 
the electric field. This can be understood by a 
slight reinterpretation of the thermodynamic 
formulae. C, = C_ does not mean that half of the 
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ions are positively, and the other half negatively 
charged, so that both types will move, only in 
opposite directions, but that every protein- 
particle is alternatingly positively and negatively 
charged, with the charge 0 in the time average. 
If the successive changes in the sign of the charge 
succeed each other fairly rapidly, the small 
motion in both directions becomes unobservable, 
the particles seeming to have the mobility 0. 

If small particles of quartz or glass of micro- 
scopic dimensions but otherwise of arbitrary 
shape and size be immersed in the liquid it has 
been found by Abramson‘ that these particles 
coat themselves with a surface layer of protein, 
and subsequently show identical mobility in an 
electric field. Furthermore, the isoelectric point 
of the microscopic particles, ie., the pH at 
which they show the mobility 0, coincides with 
the isoelectric point of molecularly dissolved 
protein. It is this general behavior which we are 
here trying to explain. 

Let us denote by y’‘, u.', w_*, respectively, the 
specific thermodynamic potential of the neutral, 
the positively, and the negatively ionized protein 
molecule in the surface of the particle denoted 
by z We are considering for the following 
purposes the surface layer of adsorbed proteins 
as a phase in the thermodynamic sense, which 
assumption is very nearly correct, since the 
surface layer contains a large number of mole- 
cules. Removal or addition of a protein particle 
would influence the thermodynamic properties 
but very slightly. We shall furthermore assume 
that the individual specific thermodynamic po- 
tentials of the neutral, positively, and negatively 
charged protein particles apart from the concen- 
tration terms are all changed by the same amount 
through the process of adsorption. This assump- 
tion seems justified in view of the fact that the 
protein particles are giant molecules containing 
thousands of atoms, so that the change in free 
energy produced by adsorption will probably 
affect very little, if at all, the ionogenic properties 
of the protein. Our second assumption would be 
untenable if we were dealing, for example, with 
comparatively small molecules or ions, with 
which the adsorption would strongly influence 
the ionization properties. 





*H. A. Abramson, J. Phys. Chem. 35, 289 (1931); 
J. Gen. Physiol. 15, 575 (1932). 


In the case of thermodynamic equilibrium 
the following relation, as one out of the series of 
many similar, must hold 


Bei peAi=py=py*+edy. (6) 


In Eq. (6), A stands for the electrical potential 
difference in the thermodynamic sense. Eq. (6) 
simply states that there must be equilibrium 
between any ion on the surface of any particle 
and the same ion in solution, and therefore 
between any two identical ions on the surfaces 
of two particles selected arbitrarily. According 
to our assumption as to the influence of adsorp- 
tion we are justified in equalizing uw,‘ and p,*, 
so that Eq. (6) reduces to 


Ai= Ax. (7) 


This relation was to be expected from the general 
behavior of ‘large electrodes” in electrochemis- 
try. 

Nevertheless, since A might still contain, as 
mentioned above, potential differences due to 
double layers rigidly connected to the particle, 
Eq. (7) does not express the equality of the 
electrophoretic potentials (i.e., the equality of 
the mobilities). To derive this second result we 
write the equilibrium condition Eq. (6) in the 
form 


bso’ +RT log Cyi+eAi=yyot+RT log Cy, (8) 
uot + RT log C_i—eA;=y-0t+RT log C_, (9) 


where we have assumed, as is certainly correct, 

that we are dealing with an ideal mixture. We 

subtract Eq. (9) from Eq. (8), and apply the 

resulting equation 

yo0'—w_o' + RT log (C,*/C_*) +2eA; 
=p+0—M-o+RT log (Cy/C_) (10) 


to the isoelectric point 
Ci=C_, Cy*a C*, 


At this isoelectric point there is no more excess 
charge present on the particle (mobility O), so 
that any difference in potential would be due to 
double layers. Rewriting Eq. (10) we obtain 


2eAi = w40—M+0'— (u-o— H-0'), (11) 


in which the right-hand side is 0 (or at all 
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events very small), according to our second 
assumption. This shows for our model the ab- 
sence of potential differences due to dipole layers, 
and therefore leads, with the help of (7), to the 
desired result: the potential difference due to 
excess charge and (therefore the mobility) is the 
same for all particles, independently of size and 
shape. 

It becomes quite obvious in our presentation 
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that the identity of the isoelectric point for 
free and adsorbed particles is by no means 
necessary, and in the case of proteins is due to 
the fact that adsorption apparently does not 
influence the ionogenic properties. In more 
general cases in which the dissociation is changed 
by adsorption the identity of the isoelectric 
point is not to be expected, and then our conclu- 
sions would become invalid. 
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The Thermal Dissociation of Cyanogen into Cyanide Radicals. A Correction 


(J. Chem. Phys. 1, 432, 1933) 


E wish to correct a number of numerical In the second column of the same page, on 
errors which are present in our article of line 3, 147 must be changed to 146 and on line 
the above title. 32, 161 changed to 164. 
On page 438, column 1, line 9, 567 should be The conclusions reached before remain sub- 
written instead of 563, and on the following line stantially unaltered by these corrections. 
73 instead of 77. In the same column, lines 5 and G. B. KisTIAKOWSKY 
4 from the bottom, 234 should be corrected to H. GERSHINOWITZ 
238; on the next line 22 changed to 27, and on Harvard University 


the last line 153 to 152. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


not later than the 15th of the month preceding that of the 
assue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


The Formation of Formaldehyde by the Action of Ultraviolet Light on Carbon Dioxide and Water: An Application of 
the Allison Magneto-Optic Apparatus 


Since the Allison' magneto-optic apparatus gives read- 
ings specific for each compound in solution, regardless of 
others that may be present, and since it will detect concen- 
trations as low as 3 or 4 parts in 10", it should prove to be 
a powerful research tool in dealing with many problems in 
chemistry, physics, and biology. We have recently used 
this apparatus for the detection and quantitative determi- 
nation of formaldehyde formed by the action of ultraviolet 
light on carbon dioxide and water. Formaldehyde gives 
minima at scale readings 21.83 and 21.92 Allison units. 
These minima are not found in pure water nor in carbo- 
nated water that has not been exposed to ultraviolet light. 
Upon placing the carbon dioxide saturated water in a 
quartz flask and irradiating the solution for several hours 
with the light from a mercury arc lamp, we obtained the 
light minima characteristic of formaldehyde. We have 
made ultraviolet light exposures up to 25 hours and find a 
uniform increase in the formaldehyde concentration with 
an increase in irradiation period. That this reaction is 
accompanied with a liberation of oxygen is indicated by 


the fact we have been able to oxidize (simultaneously with 
the formaldehyde formation) certain metallic ions from a 
lower to a higher valence. ‘‘Blanks’’ were run in all experi- 
ments to insure the absence of higher valence ions at the 
start. 
This work is being continued and a full account will be 
published elsewhere upon completion of the investigation. 
Joun H. Yor 
R. E. WINGARD 
Cobb Chemical Laboratory, 
University, Virginia 
and 
Department of Physics, 
Alabama Polytechnic Institute, 
Auburn, Alabama, 
October 28, 1933. 


1 Allison and Murphy, J. Am. Chem. Soc. 52, 3796 
(1930); Allison, Ind. Eng. Chem., Anal. Edition 4, 9 
(1932); Allison, J. Chem. Ed. 10, 71 (1933). 


Helium Hydride Ion 


The existence of helium hydride ion had been predicted 
on the basis of considerations regarding isoelectronic 
structures! and this ion has actually been found by the 
methods of mass-spectroscopy by Bainbridge.? Similar 
rare gas hydride ions have been found by Lukanow and 
Schiitze.* 

The structure of this ion has not yet been considered 
on the basis of wave mechanics and we thought the 
appropriate calculation should be made. The problem was 
considered in two ways: 

(1) We supposed that the structure resulted from the 
interaction of an excited helium atom 


v2 = (b+er2)ePr2 


with a proton. These ¥-functions for the two electrons of 
excited helium were found satisfactory when used to 


vy = on™; 
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calculate the total energy of the helium atom (24S, 19.77 
e.v.) with the values of the parameters a and 8, 2 and 0.7, 
respectively. However, it was found that at all distances 
of separation of the atom and the proton repulsive forces 
predominated. 

(2) We next assumed that the structure of helium 
hydride ion was made up of one helium ion and a hydrogen 
atom. Applying the method of Rosen‘ for the unsym- 


metrical system with a=2, B=1 and m=n=1 which 


1G. Glockler, D. L. Fuller and Charles P. Roe, J. Chem. 
Phys. 1, 703 (1933). 

2K. T. Bainbridge, Phys. Rev. 44, 59 (1933). 

3H. Lukanow and W. Schiitze, Zeits. f. Physik 82, 610 
(1933). 

4N. Rosen, Phys. Rev. 38, 255 (1931). 
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implies that the y-functions used are purely hydrogenic for 
the normal state, we repeated the calculations. We evalu- 
ated all the integrals with exception of Jz (page 272 of 
Rosen’s paper) and now we find that there is an energy 
minimum of 8.1 e.v. at a distance of 1.3 Bohr hydrogen 
radii. The effect of the omitted integral Jz will be to reduce 
the value of the energy. 

The calculation of the integral J; mentioned above 
appears very laborious and at the present time it is not 
possible for us to spend the necessary time in its evaluation. 
However, we believe from physical considerations that its 
inclusion in the calculation will not change the main result 
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of our work, i.e., that a minimum exists in the potential 
energy curve of helium hydride ion and that therefore wave 
mechanics agrees with experiment in making HeH* a 
stable structure. 

It is a pleasure to thank Professor E. L. Hill of the 
Department of Physics for the advice given and the 
interest he has shown in these calculations. 

GEORGE GLOCKLER 
D. L. FULLER 
University of Minnesota, 
Minneapolis, Minnesota, 
November 13, 1933. 
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